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SUMMARY 


This  is  the  third  report  in  a series  of  three  published  at  the  Seismic 
Data  Analysis  Center  in  1974,  which  describe  our  studies  of  and  programming 
experience  with  the  ILLIAC  IV  computer.  The  present  report  is  a computer 
listing  ol  the  ILLIAC  IV  version  of  a scientific  program  called  FKCOMB . The 
main  program,  FKCOMB,  and  two  data-editing  and  formatting  modules,  DEMI  and 
DEM2  were  written  in  Computational  Fluid  Dynamics  Code  (CFD) ; some  subroutines 
were  written  in  ASK  code. 

The  first  report  in  the  series,  SDAC-TR- 74-16,  presents  an  overview  of 
the  ILLIAC  IV  System,  describes  the  suitability  of  the  ILLIAC  IV  computer  as 
a processor  of  seismic  data,  and  contains  project  notes  on  programming 
techniques  and  languages.  The  second  report,  SDAC-TR-74-17,  is  a complete 
documentation  of  the  preliminary  version  of  the  FKCOMB  software  designed  for 
processing  long-period  seismic  data;  it  also  compares  results  from  the 
ILLIAC  software  with  those  of  the  original  version  of  FKCOMB. 


I  <\E  RR>CFU  .DEMI  ( 1 M()N  8-JIJL-74  8«46AV 


4AIN  DRIVER  FOR 


FNTERED  PY  GENE  WAGENBRFTH,  APRIL  24,  1974 


★ 9E  INTEGER  CNTRL(*,6) ,OUTBUE(*,64,6)  ,PINTI  ( * ) , INBUFC  *,  I 2B)  , 

I T I ME  ( ★ ) , OLDT I M ( * ) , 

1 SAVBCT, SAVPTN , OIJPAGE ( 6)  ,TGTEPS(6)  ,SCANS, 

2 OUPTWA ( 6 ) , OTIMHA(A) 

*CU  I NTEGER  ADBPUF(R)  , ARRAY,  I NPTR,  I NPTN ,S A V ADR , ADBOUT( 6 ) , OlJPTW , 

1 PYTS , WORDS, Tl  ,T2,T3,T4,T5,T6,  IT.PRTI AL.ADDRS, 

2 WORD,  BYTCNT(6),ADBWR0, I N3YT , OIIPYT,  ORGCOR,  PAGE  , 

3 DEBUG, 3CT,ADP, ENDADR 

★CU  LOGICAL  LADBBIK  8) ,LARRAY, LI NPTB ,L I NPTW , L5AVAD,LADPOU( 6) ,L<MJPTW, 

1 LBYTS ,LWORDS , LTI  , L72 , LT3, 1.T4 , LT5 , LT<A , LOUR YT , L I T , LPRT I A , 

2 LADDRS  ,LWORD,I.I  NBYT,LBYTCN(  6)  ,LADRlN'R,LoRGC(),LPAGF 

3 ,LDF  JUG, L3CT, LADB, LENDAD 
★EXTER JAL  RDPRM ,CETP YT, PUTRYT ,CNVTI M 

★COMMON/MA I N/CNTRI.  ,0'UTbUF  , INBDF  ,PI  NTI  ,TI  MR  , OLDTI M, SAVRCT,SAVPTW , 

I TSTFP5,  SCANS,  OUPTWA, OUPAOE,  OTI  ME  A 

★ FUJI  VALENCF(  1 , ADRBUF  ( I ) ,LADRHU(  I ))  , ( 9,  ARRAY  ,I.ARRAY ) , 

I ( 10, I NPTB, LI NPTR) , 

I (11  , 1 NPTW ,L INPTW)  ,(  12,SAVADR,LSAVAD), 

1 ( 1 3 , ADBOUT ( I ) ,LADR()IJ(  1 ) ) 

2 ,(  I 9,  OUPTW  ,L()'JPTW)  , ( 20,  PYTS , LBYIS  ) , ( 2 I , WORDS, LNORDS)  , 

3 ( 22 , T 1 ,LTI ),( 23, T2,LT2),(  24, T3,L73) ,(25, T4,LT4), 

4 (26,T5,LT5),(27,T6,LT6)  , <2B,OU8YT,LOUBYT)  ,(29,  IT,  LIT) , 

5 ( \),PRTIAL,LPHTIA) ,<31 ,ADDRS,LADDRS>, <32,WORD,LWORD>, 

6 ( ,3,IN3/T,LINBYT) ,<34,BYTCTT<  I ) , LBYTCN ( I )), 

A (40,ADBA'RD,LADB1"R) 

7 ,(43,()RGCOR,LORGCO)  , 

8 (44,  PAGE,  LPAGE),<  45, DEBUG,  LDEB'JG) , < 46, BCT,LBCT) , 

9 ( 47 , ADR , LADB ) , ( 4R , ENDADR , LENDAD ) 

★ DISK  AREA  0UPUT1  ( 20)  ,0UPUT2(  20)  ,0UPUT3( 20)  , 0UPUT4  ( 20)  ,()UPUT5  (20)  , • 

1 0UPUT6(20)  , I*JPUT(  1 00) 


1 ( 1 3 , ADBOUT ( I ) ,LADR()IJ(  1 ) ) 

2 , (19,  OUPTW  ,L()'JPTW ) , ( 20, PYTS , LBYTS ) , ( 2 1 , WORl 

3 ( 22 , T 1 ,LTI ),( 23, T2,LT2>,( 24, T3,LT3) ,(25, T4,L 

4 ( 26, T5 ,LT5 ),(27,T6,LT6) ,(28, OURYT,LOUBYT) , ( < 

5 ( \!,PRTIAL,LPWTIA) , ( 3 1 ,ADDRS, LADDRS), (32, WO! 

6 ( ,3,IN3/T,LINBYT) ,<34,BYTCTT<  I ) , LBYTCN ( I )), 

0 (40,ADBA'RD,LADB1"R) 

7 ,(43,ORG( 

8 ( 44, PAGE, LPAGE), <45, DEBUG, LDEB'JG)  ,(46, BCT, LI 

9 ( 47,  AD3, LADB)  ,(  48, ENDADR, I.ENDAD) 

♦ DISK  AREA  0UPUT1  (20)  ,0UPUT2(20)  ,OUPUT3(20)  ,()UPUT4(20)  ,( 

1 0UPUT6(20)  , I*JPUT(  1 00) 

M()DE=ON 

ENDAl)b=-l 
ORGCOR--HI 93 
JUMP  AROUNDA RE A I 
DI  SPA«  AREA  "DISPA''| 

AROUNDAHEA » » OPNDI SP  DISPAl 
THAT  IS  USED  FOR  DISPLAY  OUTPUT. 

DISPLS  " ",  16,PEGINHEADER,ENDHFADER-1 * 

JUMP  ENDHEADEkt 
3LKI 

3EGI NHEADER* • » 

DATA  (("********") 8, 3338)2, 

"DATA  EDITING  MODULE  1 VERSION  1.1", ODOA * I 5, 

( ( »**-*****★»  )8, 3338)21 
ENDHEADEkt *» 

★CALL  RDPRM 

SUBROUTINE  RDPRM  INITIAI  IZES  VALUES  TO  BE  USED  BY  THIS 
FOLLOWING  MODULES. 

*CALL  GETEYT 


G< 
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HAS 

,7 


JUST  BEEN  HEAD  IN.  NOW  TO  LOOK  IT  UP. 


ID  WAS  NOT  ON  OUR  LIST. 


A 

A 

A 

A 


THE 


FROM 


THE  HEADER  ID 
*DO  20  ARRAY= I 
★IF  C ARRAY . FO . 7 ) GO  TO  850 
IF  WF  GET  THAT  FAR,  THE  HEADER 
I'6*CNTRL(  I I .ARRAY) 

*IF(T6.EQ.I NHYT) GO  To  25 
THAT  MEANT  WF  FOUND  IT. 

20  ★CONTINUE 
25  *CONTIHUF 

★IF  (DEBUG.LT.  I ) CO  TO  35 
DISPI.S  « ",  I6.BI.EI-H 
SKIP, El  I 

35  ★CONTINUE 
40  r<S»CNTRL(l  .ARRAY) 

★IF  (T6.EQ.O)GO  To  45 

; rHAT  CHECKED  TO  SEE  IF  THERE  IS  TIMING  INFO  TO  GET  FROM 
- HFGINNING  OF  THIS  RECORD. 

★CALL  CNVTIM 

tiuPbegihning^of  Si  SLAND  CHNVERTS  ,T  T0  “eciseconds 

45  ★CONTINUE 

‘ )LDT  I .'/(★)  = OT I MEA  ( A RRAY) 

"5™  1*1  P^eissio!  ,N  nLDTI" HHE“  ™E  LAST  REC[,RD 

r6=CNFRL(4, ARRAY) 

I NPTB= I NPTB+T6 

S”STss^7{„TS T0  THE  BE0INNIRC » A «*  **»• 

SCANS-0 

50  *IF  (DEBUG.LT. I ) GO  TO  55 
DISPLH  "T-SCAN" ,0| 

55  T6=CNTRL(5, ARRAY) 

I NPTB= I NPTB+T6 
Tft-CNTRL< I .ARRAY) 

★IF (T6.EG. I )G0  TO  65 
IF  RE  GET  HERE  THERE 
★CALL  CNVTIM 
T6=CNTRL( 6, ARRAY ) 

I NPTB= I NPTB+T6 

65  *IF  ( DEBUG.LT. I ) GO  To  75 
DISPLH  "TIME”, 16, TIME, TIME! 

DISPLH  "OLDTIM" , I 6.0LDTIM.0LDTIMI 
75  *C()NTINIJE 

MOW  CHECK  TO  SEE  IF  A TIME  SCAN  IS  MISSING. 

★ IF ( . AMY.  ( ( OLDTIM (*)  .ECl.O ) ) )G()  TO  200 
★IF( .ANY. ((TIME(*).LE.0LDTIM(*)+I5) ))G0  TO  200 
★IF( . ANY. ( (TIME( ★) ,GT.OLDTIM(*)+25 ) ) )G()  To  150 
IF  WF  GET  HERE  THERE  WERE  I OR  2 TIME  SCANS  MISSING.  WF  HAVP  To 
rILL  IN  THE  MISSING  TIMF  GAP(S)  WITH  THE  LAST  TIME  GAP,  WHICH  HAS 


IS  A TIME  WORD  WITH  THIS  SCA. 
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' BEEN  CAREFULLY  SAVED  IN  OUTBUF.  WE  DO  So  A BYTF  AT  A time  iimtti 

; all  channels  are  done,  this  is  not  the  fastest  wJy  to do the um 

: COMP  LI  ''ATFD 1 A HR I T°R  v A?UP  ApD  5ASILY  DEBUGGABLE.  THE  SITUATION  IS 

: ADBOUT( ARRAY)  BIT  RY  THF  FA°T  ™AT  S"ME  °F  THE  BYTES  ARF  IN 

SAVADB=ADBOUT( ARRAY) 

SAVBCT=bYTCNT( ARRAY) 

T6=OIJPTW-| 

SAVPTW=T6 
T6=CNTRL( 3, ARRAY) 

T5=SAVBCT 

BYTS=T6-T5 

■ CHANNELS^M?NuISTHE  N^bIr^BUEs”^?  “UAL  T"  ™F  ^ "F 

1 LWORDS=LBYTSSSHRE2ENT  To  FoRCE  DEALL0CATI()N  of  REGISTERS. 

LT I * LW  ORDS . S ML . 2 * *' 

PRTlAL«BYTS-TfK  STATEMENT  T°  FoRCF  DEALLOCATION  OF  REGISTERS. 

T6=SAVPTW 

ADDRS=T6-WORDS 

ADDHESS  IN  OUTBUF  OF  PARTIAL  WORD. 

J*  JuW,,ULD  LIKE  TO  HAVE  OUTBUF  DIMENSIONED  ••  OUTBUF! 81 92.8)"  AND 
SIMPLY  ACCESS  "OUTRUF( ADDRS, ARRAY) " BUT  CFD  FORCES  ME  TO 
DIMENSION  OUTBUF  "OUTBUF<*,64,8) « AND  HAVE TO DO  SOME  ARITMMPTir 
HERE  TO  CALCULATE  THE  PRoPPER  INDICES.  ARITHMETIC 

Tl =ADDRS 

L‘l2  - OFF  Af,N . t"INT  T°  F0RCE  DEALU,CATI1)N  <*  ROISTERS. 

LT2=LT2.AND.LTI 

LTI=LTI .SHL.6 

^Tsoutrufc^^t^  array?)ENT  T0  F0RCE  DEALL0CATI(,N  (,F  registers. 

★IF  (DEBUG.LT.*!  )*GO  TO  105 
DISPLH  "TIME-GAP", Of 
DISPLH  "SAVADB" , I 6,SAVADB,SAVADB» 

DISPLH  "SAVBCT", I6,SAVBCT,SAVBCT* 

DISPLH  "SAVPTW", | 6,SAVPTW,SAVPTW I 
DISPLH  "BYTS" , I 6,BYTS,BYT5| 

DISPLH  "WORDS", I 6, WORDS, WORDS! 

DISPLH  "PRTIAL", I 6,PRTI AL, PRTI ALl 
DISPLH  "ADDRS" , I A, ADDRS , ADDRSi 
DISPLH  "IT", 16, IT, IT! 

05  *CONTINUE 

NOW  WE  COME  TO  ANOTHER  CFD  ABOMINATION.  I WANT  TO  DIVIDF  ONC 
IN  MEMORY  by  4.  I CAN  EITHER  MOVE  IT  To  A VECTOR  OR 
f^^LE  WITH  THE  MODE  OR  SOMETHING  AND  Do  A STRAIGHT  DIVIDE 
OR  I CAN  MOVE  IT  TO  THE  CU  AND  JUST  SHIFT  IT.  I OPT  FOR  THE  LATTER 
AND  CONTINUE  TO  OPT  FOR  THAT  THRUOUT  THE  PROGRAM  THIS  IS  RETAIISt 
SOMEDAY  CFD  MAY  ALLOW  ME  TO  DO  WHAT  I WANT  OR  I CAN^PUT  IT  ?N c!®| 
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C STATEMENTS  MYSELF. 

\ . DUMMY  ASK  STATEMFNT  TO  FORCE  DEALLOCATI ON  OF  REGISTERS. 

LPRTI A-LPRTI A .SHL.2 
LIT* I I T . RTR.  PkTl  AL 

\ -.DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 

PI NTI  ( *)*TI ME(  ★) -OLDTI V( ★ ) 
r6=PINTl ( I ) 

★DO  140  Tl* I »T6t I 0 

C ONCE  FOR  EACH  MISSING  TIME  GAP 

C FIRST  DO  THE  PARTIAL  WORD. 

ri  ME  ( * ) XTI ME  ( ★ ) ♦ 1 0 
C L()UBYT=T!ME 

A SLIT(O)  =TIMFl 

A LOAD(O)  SCO! 

A CSHR(O)  161 

A STL(O)  LotJRYTi 

★CALL  PUTBYT 
A 5LIT(0)  TIMFl 

A LOAD(O)  $C0» 

A LIT( I J-OFFFF*  I 61 

A CAND(J)  SCIi 

A STL(O)  LOIJRYTI 

★CALL  PUTBYT 
★ DO  I 10  T2=l .PHTIAL.I  6 
C ONCE  FOR  EACH  BYTE. 

LIT*I.I  T .RTL.I6 
LOURYT*OFF .TURN  ON.. LAST.  I 6 
L0UBYT-L0URYT.AND.LIT 
I I 0 *CALL  PUTBYT 
C NOW  DO  THF  FULL  WORDS. 

T6-SAVPTW 

★DO  I 20  T2=ADl)RS+l  fT6 

C NOW  WE  HAVF  TO  DO  THF  MESSED  UP  ADDRESSING  BECUASE  OF  THF  CFD 
RESTRICTION  ON  DIMENSIONING. 

• DUMMY  ASK  STATEMENT  To  FORCE  DEALLOCATION  OF  REGISTERS. 
LT3-LT2 .SHR.A 
LT4=0FF.TtJRN  ON.. LAST. 6 
L 14=1  T4.AN’D.LT2 

A - DUMMY  ASK  STATEMFNT  To  FORCE  DEALLOCATION  OF  REGISTFRS. 
W')RD=()UTBIJF  ( T4  fT3t  ARRAY ) 

*00  120  TS=  I , 4 
LNORU-LWokD.RTL.16 
L0UBYT*()FF.TURN  ON . . LAST.  I 6 
LotJBYT-LOlJPYT.  AND.LWORD 
120  *C A LL  PUTBYT 

C NOW  PUT  OUT  THE  BYTES  IN  SAVADP. 

WORD-SAVADR 
T6-SAVBCT 
★DO  140  T2-I ,T6 
LOUBYT-OFF.TURN  ON.. LAST. I 6 
LOUBYT-LOUBYT. AND.LWORD 
★CALL  PUTBYT 


9* 


I <KF RH>CF  ).DFM I | | 


MOM  8- JlJL-74  Hi 4 6AM 
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A B1JMCH  OF  TIM"  GAPS  MISSING.  WF  NFFl) 


THAN  2 TIME  GAPS.",0D0AH6, 


C 

c 


140  L;iord*lw)RD.rTL.I6 

*30  To  280 
150  *CONTINUF 

IF  W E GET  HERE,  THEHE  WY.'ERE 
C ONLY  PUT  OUT  A DIAGNOSTIC 

i 1’  "•  '^bcapdiagnostic.egapdiagnostic-i  l 

A DISPL:!  "OLDTI M"  , I 6.0LDTI  V , OLDTIMi 

A DISPI.H  "TIME",  I 6,  TIME,  TIME  I 

A JUMP  EGAPDIAG.NOSTICl 

A Hl.Kl 

A B3APDI AGNOSTIC i» | 

A DATA  f ("********»' )G,  3338)  10,  "MORE 

A (("********") 8, 3338) I Ol 

A EGAPDI AGNOSTIC i i | 

I 75  *CONTINUE 
200  *CONTINUE 

C :J0W  T()  CHECK  FOR  A TIME  REVERSAL. 

*IF(  .ANY . ( (TIME ( * ) ,GT.oLDTIM( *) +5  ) ) )rj()  To  250 
IF  fit-  CET  HERE  THERF  IS  A TIME  REVERSAL.  IF  THERE  IS  A TIMF  wriDn 

liHni3CARDEo!S  SCAM  IS  THWn,1,N  °UT*  f,THERrtISE»  THE  WHOLE  RECORD 
T6*CNTRL( I .ARRAY) 

♦IF  (T6 ,E0. I ) Go  To  210 
•*E  HAVE  A TIME  if.ORD  PER  SCAN. 

?NPTB=PINTlT^?*C^T^L*  ^ RPAY)+CNTRL(  7, ARft AY) ♦CBTPL( 3,a  RRAY) 

NOW  SKIP  AROUND  PROCESSING  OF  THIS  SCAN. 

♦GO  TO  410 

recorjNDLE  THRf)WING  AWAY  THE  m)LF-  REC(,RD  FnR  time  WORD  PER 

210  PINTI(*)=CNTRL(4,aRRAY)*CNTRL(9,ARRAY)*(CNTRL(S.ARRAY)+ 

?NPTP* P I NT?< I ) +CNTRL { 7 * A RRAY  } ] ♦CNTRL ( 3 , A RR  AY ) ♦ I N PTB 

I6*SCANS 
T6=T6+ I 00 
SCANS= T6 
♦GO  TO  410 

THE^DATA71  VF  IS  ALL  TAKEfJ  CARE  °F*  plWST  WE  OUTPUT  TIME.  THEN 
250  *C()NTINUE 
: Lo'JHY  r=TIMF(  I ) 

A 3LIT(0 ) TIME  I 

A LOAD(D)  $C0| 

A CSHR(O)  161 

A STL ( 0)  LOUBYTl 

♦CALL  PUTBYT 
A SLIT(O)  TP'El 

A LOAD(O)  SCOl 

A LIT(I)  =OFFFFi I 61 

A CANC(O)  SCII 


C 

C 


C 

C 


ALL 


S' 


< 


o o > > o o o o o 
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A STL(O)  LOUBYTt 

★CALL  PUTBYT 
()LUTIM(*)>TIME(*) 
r I ME ( * ) *T1 ME ( ★ ) + 1 0 
r6=TC  rEPS( ARRAY) 
r 6=  r i 

rSTEPS (ARRAY )=T6 
T6=CNTRL<  3, ARRAY) 

★DO  400  Tl  = | ,T6 

C THATS  ONCE  FOR  EACH  CHANNEL. 

★CALL  GETBYT 
OU BYT* IN BYT 
★CALL  PUTBYT 
400  ★CONTINUE 
410  T6*SCANS 
T6=T6+I 
SCANS=T6 

T6*CNTRL<7f ARRAY) 

I *IPTB*  INPTB+T6 
T6*CNTRL(9f ARRAY) 

T5*SCANS 

★ IF (T5.GF  .T6)G()  To  420 
★GO  To  50 

420  r6=CNTRL(8, ARRAY) 

I NPTH= INPTB+T6 
OTIMEA(ARRAY)  = ()LDTIM(  I ) 

★GO  To  10 
850  ★CONTINUE 

★IF(INHYT.EQ.O)GO  TO  855 
A DISPLil  "BADHEADi"t2t 
855  ★CONTINUE 

ADEN  RE  GET  HERE  ALL  THE  DATA  HAS  BEEN  PROCESSED.  TIME  To  EMPTY 
THE  AJB  BUFFERS,  KRITE  OUT  THE  CORE  BUFFERS  , PUT  THE  HEADERS  IN 
THE  OUTPUT  FILES  AND  LAY  BACK  AND  QUIT. 

FIRST  EMPTY  THE  ADB  BUFFERS. 

★DO  86J  ARRAY= I ,6 
BCT=4-3YTCNT (ARRAY) 

•o  DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
LBCT*LBCT.SHL.4 

« DU*MY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 

LA DB=LADBOU( ARRAY) .SHL.BCT 

HERE  RE  HAVE  TO  DO  SOME  MORE  FANCY  ADDRESSING  BECAUSE  OF  THE  CFD 
RESTRICTION  ON  DIMENSIONS. 

Tl  =OIJPTR 

A N DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTFRS. 

LT2-LTI .SHR.6 

LT 6= OFF. TURN  ON.. LAST. 6 

LT1= LT6.AND.LT I 

A -v  DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 

OUTBUF (T I ♦ I ,T2+I ,ARRAY)=ADB 
860  ★CONTINUE 


11 


- 


MM 


3»  >■ 


■ 
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C NOW  WRITE  OUT  ALL  OF  THF  BUFFERS ( 6 AT  THIS  TIME). 
r6=OUPAGE(  I ) 

*WRITE(  64  ,()UTF!IJF(  1,1,1  ) ,OUPUTI  (T6) , 4) 

♦WAIT  64 
T6*01JPAGE(2) 

♦NR I TE(  64 , OUTBUF (1,1,2), OUPUT2 ( T6 ) , 4 ) 

♦WAIT  64 
T6=OlJPAGE ( 3) 

*WRITE(  64 , OUTBUF ( 1 ,1 ,3)  ,OUPUT3(T6)  ,4) 

♦WAIT  64 
T6=OUPAGE(4) 

*WRITE(  64  , OIJTB'JF  (1,1,4),  OUPUT4(T6) , 4) 

♦WAIT  64 
T6=OUPAGE ( 5 ) 

♦WRI TE ( 64 , OUTBUF  (1,1,5),  OUPUT5 ( T6 ) , 4 ) 

♦WAIT  64 
T6=OUPAGE ( 6) 

♦WRITE(64,OUTBUF( 1,1,1 >,0UPUT6(T6),4) 

♦WAIT  64 

C NOW  FOR  THE  HEADER  PAGE.  FIRST  CLEAR  A BUFFER. 

♦DO  930  Tl=l ,16 
930  OUTBUF < *,TI , I )=0 
C NOW  FILL  THEM  IN  ONE  AT  A TIME. 

OUTBUF  (I  ,1,1 )=CNTRL(  11,1) 

C THAT  WAS  THE  HEADER  ID. 

OUTBUF (2, 1 ,1 )=TSTEPS ( I > 

♦WRITE ( 64,  C)UTBUF(  1,1,1)  ,0UPUT1  ()),)) 

♦WAIT  64 

OUTBUF ( I ,1,1  )=CNTRL( 1 1 ,2) 

OUTBUF  (2, 1 , I )=TSTFPS  ( 2 ) 

♦ WRITE(  64, C)UTBUF(  1,1,1)  ,()UPUT2(  I ) , I ) 

♦WAIT  64 

OUTBUF ( 1,1,1  )=CNTHL( 11,3) 

OUTBUF (2, 1 ,1  )=TSTEP5( 3) 

♦WRITE!  64 , OUTBUF ( 1,1,1) ,0UPUT3(  I ) ,1  ) 

♦WAIT  64 

OUTBUF ( I ,1,1  >=CNTRL(  11,4) 

OUTBUF (2, 1 ,1  )=TSTEP5(  4) 

♦ WRITE!  64, OUTBUF ( 1,1,1) ,0UPUT4(  | ) , ) ) 

♦WAIT  64 

OUTBUF (I ,1,1 )”CNTRL(  I 1,5) 

OUTBUF (2, I , I )=TSTEP5  ( 5 ) 

♦WRITE  ( 64 , OUTBUF  (1,1,1)  ,0UPUT5(  l),l  ) 

♦WAIT  64 

OUTBUF ( 1,1,1 )=CNTRL(  11,6) 

OUTBUF (2, 1 , 1 )=TSTEPS  ( 6) 

*WR ITE(  64, OUTBUF ( I ,1,1 ) ,0UPUT6( I ) , I ) 

♦WAIT  64 

UISPLS  11  I 6,BFINAL,EFI  NAL-I  » 

JUMP  EFINALl 
BLKl 


12< 
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A 

A 

A 

A 

A 

A 


BFlNALu 

DATA  ( ("*★*★★*★*») 8, 3338)  10, 

"MOD  I GOING  TO  END  OF  JOB" , GDOAt  I 6, (("★★★★*★**« )8, 3338) 1 Ol 


EFINALi i t 


CLSDISP  DISPAl 
★STOP 
★END 


m o r-  ao  o o 
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C 

C 

C 

C 

c 

c 


A 

A 

A 

A 

A 

A 

A 

A 

A 


DEM2  - DATA  EDITTING  Mnmn  p 3 


★PE  INTEGER 


★PE 

I 

★PE 

I 

★PE 

★CU 

1 

2 
3 

★CU 

1 

2 
3 


CHGOODC  80 ) ^IlTEGD (80 ) * S I TES  { ( *♦ 640 } ♦ ABUFF2 < 704 00 ) , ■ 

? Hit  J:t!MS?i80 * ♦ R0UFF 1 ( 4096 ) , ROWSUM , RBUFF2  ( 70400) 

INTEGER  ADBBUF(S)  .COREPT,  BYTE • ADBWRD  ARRAY  DFRUG  Tw^ 

(1VLAP,NCHAN,NSITH,NR()WS,DIFFH,DIFFW.NEN  ’lLDB0*pTTSr»M''  ■ 
INDEX!  .INDEX2,  INDEX3, INDEX4,T!  It2,T3!t4  T5  T6  CH*TpArp'>*" 
OFFSET,  I NBYT,NGDCH,TWSZR,NGDST|NGDRFRF3PFNGT^oPArF^T7* 
l«0,cal  ^^|®U*®’*j-^<,REP,LASTls|LBYTE:LADBW:LARRAySiBUG°E,T7- 

‘EXTERNAL  GTDATA^O I ATA‘^r * ENGDR, LNGT *LT7 * 

♦caIEkNAL  GTDATA tC I 6T64.C64T32. ROWSUM  PIINPPT  cxotaa 

‘C0MM0N/MAIN2/NBUFF I , BUFF2 ,BUFF3, FINSCN,COMP*TOTSCN  PINT!  PINT? 

I TIME,OTIME,TNTIME,PEN,CNfRL,GLCHFT,VARFT,PRFAU^PREArJ^'  ' 
eouival^n^^ 

1 ! BUFF*?! ! >:tWFF3( ! .’!  •)‘CHC<K,Dt  1 > I*  * I 

★EQUIVALENCE  ( I .ADBBUFC I ) .LADBBUC  I.) ) , (9,C0REPT.LC0REP)  CIO  RYTP 

2 f I1!’ w?SD.,}-5?!Sr5  * ( ? 2 * A RrIy . L A RRA Y ) I J ?3?t5Iz . I 

3 FIS??’* ( f ;°MnA^L  )VLAP } ’ ( 1 5 • WHAN.LNCHAN > , C I 6, NS  ITE.  - 

j LNSITE) t( I7.NR0WS.LNR0WS) t( I8.DIFFR.LDIFFR)  CI90IPPW 

4 LDIFFW).(20.NEW,LNEW)f(2l  OLD  L OLD)  C22.GAP  LOAP)  fS*  " 

*unnc  n2EA  INDM2C20) ,()UTDM2(4I )fCONPRM(  I ) 


JUMP  PASTAREAI 
DISP2*  * AREA  "DISP2,,| 

PASTAREA  * 1 1 

OPNDISP  DISP2I 

DISPLS  , I6.BHEAD2.EHEAD2-I I 

JUMP  EHEAD2I 

BLICl 

BHEAD2  1 1 DATA 

(("“******«  )8 ,0DOA.  I 6)  I 0, ‘START  EXECUTION  DATA  EDITING  MODULE  2“, 
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A ODOAi I 6, "VERSION  2 .0” ,ODOAi  I 6, ( («★★★★★★★★" )8 ,ODOA 1 1 6)2 1 

A EHEAD2 it| 

★READ(64, 1 BUFF  I (I  ) , INDM2( I ) , I ) 

★WAIT  64 

*WRITE( 64 , 1 BUFF  I ( I ) ,0UTDM2<  I ) ,1  > 

♦WAIT  64 
Tl  = 1 BUFF  I ( I ) 

*D0  20  ARRAY= I , 7 
★IF (ARRAY .EQ.7)G()  TO  1105 
T2=CNTRL( I , ARRAY) 

★ IF(TI  .EQ.T2)C()  TO  25 

C THAT  MEANT  WE  FOUND  IT. 

20  ★CONTINUE 
25  ★CONTINUE 

FI NSCN (*)  = I BUFF  I (2) 

*READ( 64, 1 BUFF  I ( I ) ,CONPRM(  I ) , I ) 

★WAIT  64 
i)EBUG=  I BUFF  I (I  ) 

TWSZ=I BUFF  I ( 2) 

OVLAP= IBUFF 1(3) 

GLCHFT ( ★ ) = RBUFF I (4) 

VARFT(  ★)  =RRUFF  1(5) 

LOFREQ-IBUFFI (6) 

HIFRE()*I  BUFFI  ( 7) 

COMP( ★)= I BUFF  I (8) 

NCHAN=CNTRL( 3, ARRAY ) 

NS  I TE=CNTRL(  4 ARRAY) 

A /.  DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
LTWSZR*LTWSZ.SHR.6 

A a DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
PI  NT I (★)=NCHAN*TWSZR 
NHOWS-PINTI (I  ) 

DI  FFW*TWSZ-()VLAP 

A H DUMMY  ASK  STEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
LDI FFR*LDI FFW.SHR .6 
LTI=OFF .TURN  ON. .LAST. 6 
LDIFFW=LDIFFW.AND.LTI 

A X DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
LASTI6=0FF .TURN  ON. .LAST. 16 
BYTE-4 
ADBWRD=8 
C0REPT*4096 
0PAGE=2 
I PAGE*2 
NEW- 2 
OLD*  I 
GAP*  I 

*IF(DEBUG.LT.  I )G()  TO  30 
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IF  WE  DETECTED  A TIME  GAP  LAST  TIME. 


A DISPLH  " I N I T « " , 2 1 
3D  ★CONTINUE 
♦CALL  GTDATA 
riME(*)=INBYT*A5536 
♦CALL  GTDATA 
FI  ME (* )=TI ME (*)+INBYT 
I 00  ♦CONTINUE 

TSCANS»0NS  ™E  L(X,P  GONF  ™RU  °NCE  F°R  EACH  C(,MPLETE  TIME  WINDOW 

r'WTIMEU)*TIMEU> 

* Ln»0FF.TURNS0l5TELAST  I”  F°HCE  DEALL(,CATIoM  0F  REGISTERS. 
LNEW=LNEW.AND.I.fi 
Lf)LD«L()LD.AND.^n 

A NEW-NEW-fASK  STATEMENT  T°  F()WCE  DEALLOCATION  OF  REGISTERS. 
OLD=OLD+l 

♦I F( DEBUG.LT. I )G0  TO  105 
A DISPLH  "TWD" ,0l 

* I F ( DEBUG . LT . 2 ) GO  To  105 
A DISPLH  ,21 

105  ♦CONTINUE 
C NOW  CHECK  TO  SEE 

♦IF (GAP .EO. I )G0  TO  200 

C WE  CANT  OVERLAP  IF  THERE  WAS  A GAP. 

♦IF (DEBUG.LT. I ) GO  TO  115 
A DISPLH  "OVLAP" ,0l 

I 15  ♦CONTINUE 

TWTI  ME  ( ♦ ) »TWTI  ME  ( ♦ -OVLAP^  1 0 
PI NTI ( ♦ » = 0 

♦IF ( ( PEN( ♦) ,GT. 64-DI FFW))PI NT >(♦)=! 

: THETPRECEDINGHW)WTHAT  VALUES  W(,!JTED  ACROSS  A ROW  BOUNDARY  Go  To 

Tl =NRoWS-DI FFR 
♦IF(TI .EQ.OIGO  TO  140 
♦DO  130  INDEX) » I ,T1 
I NDE X2=I NDFX I *DI FFR 

ri^T21(.*w^DE?,1  *NEi^  ^*EUFF2  ( **DI  FF-W,  I NDEX2+PI NS  I ( * ) ,OLD) 

riDDifr LLrucEpDATA  WANTED  FW0M  THE  (,LD  TIME  WINDOW,  P!  ilF 
BY  G(X)d  SaTA ‘LATER  AGE  IS  °KAY  RECA,,SE  11  W,IE  FF  OVERWRITE' 
130  ♦CONTINUE 
140  ♦CONTINUE 

TSCANS-OVLAP 

* I F ( DER  UG . LT . I ) GO  TO  145 
DISPLH  "E  OVLAP", 0» 

♦IF(DEBUG.LT.2)G0  TO  145 
DISPLH  "BIJFF2I",  IB,BUFF2,BUFF24255l 
DISPLH  "BUFF22", I 6,BUFF2*550^A4,BUFF2+550^54*255 | 

145  ♦CONTINUE 

200  ♦continue 


* »»< 

. i 


L 


oo  o»  o onon 


« 


<KERR>CFU.DEM2l70 


MON  I 5-JUL-74  I 0142  AM 


PAGE  I *3 


C NOW  ITS  TIME  TO  READ  IN  AND  MOVE  A TIME  STEP. 

★IF (DEBUO.LT. I ) GO  TO  205 
A DISPLH  "TSTEP" , Ol 
205  *CONTIN!JE 
CH=  1 

★DO  2 BO  TI=TSCANS+I ,TWSZ 

C THATS  ONCE  FOR  EACH  TIMESTEP  TO  GET  FOR  THIS  TIME  WINDOW. 

ToTSCN (*)=TOTSCN ( ★)  ♦ I 
OFFSET=0 

★IF (NEW.E0.2) OFFSET* 35200 

THAT  MAKES  ADDRESSING  A LITTLE  BIT  EASIER.  WE  CAN  TREAT  BUFF2 
( ABUFF2)  AS  I DIMENSIONAL.  OFFSET  ACTS  LIKE  (NEW, OLD) • ITS  REALLY 
BECAUSE  CFD  WILL  NOT  ALLOW  A 2 DIMENSIONAL  ARRAY  WITH  THE  FIRST 
DIMENSION  OTHER  THAN  64. 

★DO  240  T2= I , NCHAN 
★CALL  GTDATA 
INDEX  I *CH+OFFSET 
ABUFF2( INDEXI  ) = INBYT 
OFFSET*OFFSET+TWSZ 
240  ★CONTINUE 

WE  JUST  DID  ONE  TIME  STEP. 

TSCANS=TSCANS+ I 
CH*CH+ I 

★IF (DEBUG. LT.2)G0  To  245 
DISPLH  " I -CHAN" ,0l 
★CONTINUE 

NOW  WE  CHECK  FOR  A GAP. 

T6*T0TSCN(  I ) 

★IF ( .ANY . ( (T6.EQ.FINSCNI* ) ) ) )G0  TO  300 
EOF  IS  REALLY  JUST  AN  INFINITE  GAP. 

GET  THE  NEXT  TIME  WORD  AND  CHECK  FOR  A GAP. 

★CALL  GTDATA 
T6*TIME( I ) 

OTIME( *)=TIME( ★) 

TIME(*)=INBYT 
★CALL  GTDATA 

TIME(*)*TIME(*)*65536+INBYT 
★IF ( DEBUG.LT. 2 ) GO  TO  255 
A DISPLH  "TIME", I 6, TIME, TIME  I 

255  ★CONTINUE 

★ I F ( . ANY.  ( (TIME(*)-T6.GT.  15)  ))G()  TO  300 
260  ★CONTINUE 

GAP=0 

★GO  TO  400 
300  ★CONTINUE 

C IF  WE  GET  HERE  WE  HAVE  A GAP. 

★IF(DEBUG.LT. I )G0  TO  315 
A DISPLH  "GAP" ,0l 

★IF (DEBUG.LT. 2) GO  TO  315 
A DISPLH  ,21 

315  ★CONTINUE 

★ IF  (TSCANS.NE  .TWSZ)G()  TO  325 
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C 

c 

c 


GET  HERE,  THERE 
THE  TIME  WINDOW 
GAP. 


IS 

IS 


A GAP,  BUT  WE 
ALREADY  FULL. 


DONT 

JUST 


NEED 

MARK 


anymore  data 

THF  FACT  THAT 


m 


c 

c 

c 

c 


325 


400 


c 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


330 


335 


C 

C 

n 


C 

C 


c 

c 


350 


c 

c 


IF  WE 
SINCE 
HAD  A 
GAP=  I 
★GO  TO 
★CONTINUE 

hinS*  n owPLirraW 

T,ME  H'mm-  1N  CASE  be" 

★IFIGAP.NE. I )G0  TO  335 
WH(X)PS,  AN  IRRECOVERABLE  ERROR. 

JUMP  EGAPMESSAGEf 
BLKj 

BGAPMESSAGF*  *DATA 
C )8,0D0A*  16)10, 

"IRRECOVERABLE  TIME  GAP ." ,ODOA*  I 6, 

(("********" )8,0D0Ai 16) I 01 
EGAPMESSAGEm 

DISPLS  , I 6,BGAPMESSAGE,EGAPMESSAGE-I  I 
DISPLH  "TIME*", I6,TIME^TIME| 

DISPLH  "OTIME*",  I 6,0TIME,  ()TI  MEl 
DISPLH  "TWT I ME  * " , I 6, TWTI ME ,TWTI MEl 
★IF (DFBUG.LT. I )G0  TO  330 
DISPLH  ,2t 
★CONTINUE 
★GO  To  1000 
★CONTINUE 

IF  WE  GET  HERE  ITS  TIME  T()  ACTUALLY  FILL  IN  A TIMF  GAP  f i dft  tut 
CURRENT  WINDOW  HAS  TO  BE  "SHIFTED  FORWARD" . THE  “(IUNt’to  JoIe 


NUMBER  OF  MISSING  SCANS-TWSZ-TSCANS 


ROUTE. 

ROUTE. 


IS  THE 
GAP*  I 

T3=TWSZ-TSCANS 
TWTIME(*)=TWTIME(*)-T3*I0 

%LT4=LT3!sHRT6TEMENT  T0  F(>MCE  DEALL0CATI,W  f,F  REGISTERS. 

LT5= OFF ! TURN  ON.. LAST. 6 
LT3=LT3.AND. LT5 
T3  IS  THE  NUMBER  OF  WORDS  TO 
T4  IS  THE  NUMBER  OF  ROWS  To 
PINTI ( *)=0 

*IF((PEN(*).LE.T3))PINTI  (★)*! 

THATS  ONE  FOR  ALL  THE  PE'S  THAT 
BOUNDARY. 

★DO  350  INDEX  I =T4+I , NROWS 
INDEX2  = INDEXI  -T4 

★CONTINUE  ,NE^  )=B^Fp2  < ★-T3, 1 NDEX  I -PINTI ( ★) ,NEW) 

NOW  COMES  THE  DIFFICULT  PART.  WE  HAVE  To  MOVE  DATA  FROM  fHF  l«sr 
TIME  WINDOW  WITHOUT  MOVING  ANY  GARBAGE,  SINCE  THIS  TIMF  WF  WOUI 
BE  OVERWRITING  G(X)D  DATA.  THE  ROUTE  AMOUNT  IS  TSCANS.  THE'nUMRFR 


ARE  CONNA  SEND  DATA  ACCROSS  A 


Ron 


j 
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C OF  ITEMS  TO  BE  MOVED! BEFORE  WE  ALWAYS  DID  THEM  ALL)  IS 
C TWSZ-TSCANS. 

C n=ROUTF  AMOUNT!  ROWS)  . 

C T4=ROUTE  AMOUNT! WORDS) 

C Tb=N UMBER  TO  DO!R()WS). 

C T6=NUMBER  TO  DOIWORDS). 

A H DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
LT3=LTSCAN .SHR .6 
LT 4= OFF. TURN  ON. .LAST. 6 
LT4=LT4 . AND.LTSCAN 
T6=TWSZ-TSCANS 

A X DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
LT5=LT6.SHR.6 
LT7=OFF.TURN  ON. .LAST. 6 
LT6=LT6. AND.LT7 

A X DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 

★ IF  ! DEBUG  .LT.  I )G()  To  355 
A DISPLH  " RT3" ,2 1 

355  ^CONTINUE 

C NOW  DO  IT  ONCE  FOR  EACH  CHANNEL,  ONCE  FOR  EACH  FULL  ROW,  AND 

C THEN  ONCE  FOR  EACH  PARTIAL  ROW. 

PI  NT  I ! *)=0 

★ I F ! ! PEN ! *)  .GT.64-T4 ) ) P INTI (★)  = I 

C NOW  WG  WOULD  LIKE  TO  DO  "DO  370  Tl « 1 ,NROWS,TWSZR"  BUT  CFD  INSISTS 

C THAT  THE  INCREMENT  BE  A CONSTANT,  So  WE  WILL  CONSTRUCT  THE 

C EQUIVALENT  LOOP. 

T I = I 

★DO  370  T7=l ,NCHAN 

★ IF  !T5.  EQ.O)G()  TO  365 

C IMPROPPER  DO  L(X)PS  MUST  BE  AVOIDED. 

★DO  360  T2=l ,T5 
I NDEXI =T2+TI -I 
I NDEX2= I NDEX I +T3 

BUFF2!  ★,  I NDEXI  ,NEW  ) = BUFF2  ! *+T4,  INDEX2+PINT!  !★)  ,()LD) 

360  ^CONTINUE 
365  *CONTI NUE 

C NOW  TO  DO  THE  LAST  ROW  OF  THE  CHANNEL.  IT  IS  ROW  T5-H  . GET  TO 

C PLAY  WITH  THE  MODE  THIS  TIME. 

I NDEX I =TI +T5 
I NDEX2= INDEX  I +T3 
C ROUTE  AMOUNT  IS  THE  SAME. 

MODE*! PEN!*)  .LE.T6) 

3UFF2!*, INDEX! ,NEW)=BUrF2!**T4, INDEX2+PINT1 ! ★) ,OLD) 

M()DE=ON 
Tl =T1 +TWSZR 
370  ★CONTI.NUF 

C THATS  EVERYTHING.  WE  NOW  HAVE  A COMPLETE  TIME  WINDOW. 

★IFIDEBUG.LT. I )G0  TO  375 
A DISPLH  "ERACKUP" ,2 1 

★IFIDEBUG.LT. 2)G0  TO  375 
A DISPLH  "BUFF2I " , I 6,BUFF2 ,BUFF2*255l 


CJOOOOO  ' ' o o C'OO 
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A375  y°(I)^j{Ju”BUFF22M»ld*BUFF2+64*550«BUF  '2+64*550+255  | 

C THATS  IT.  WE  ARE  BACKED  UP. 

400  *C()NTINUE 

TIME  hIND()W  IN  16  BIT  FORMAT  IN  BUFF2(-,-,NEW) . WF 
?!bir-L!£VE  IT  IN  NEW  AND  CAN  OVERWRITE  THE  DATA  IN  BUFF2(-.-,0LD) 
SINCE  IT  WILL  NEVER  BE  USED  AGAIN.  WE  WILL  CONVERT  TO  64  BIT 
FLCWTTNG  POINT  FORMAT  AND  'MOVE  FROM  RUFF2(  NEW,-,- ) TO  "OLD". ONCE 
THE  DATA  IN  4,OLD*'  HAS  FOUND  ITS  WAY  TO  BUFF3,  IT  WILL  BE  OVER 
WRITTEN  NEXT  TIME  AROUND. 

*D0  405  I NDEX1 = I , 550 
BUFF 2 ( ★, INDEXI .OLD) =0.0 
405  ^CONTINUE 

★DO  410  INDEX  1 = 1 ,NR()WS 

410  *CONTINUETM(BUFF2<**  INDEXI  ,NEW>  «BUFF2(*. INDEXI  .OLD)) 

?™AniS  N0W  IN  BUFF2(*.-.OLD).  IT  WILL  NOW  BE  DEGLITCH  ED,  MEAN 
SQUARE  CALCULATED  AND  CHECKED  AND  THEN  FFT'ED. 

PINT! ( ★ )=0 
PINT2(*)=0 

★IF ( ( PEN( *) .EQ. I ) )P INTI (*)  = l 
★IF ( ( PEN( ★) ,E0. 64 ) ) PI NT2( ★) *1 
ALLMSQ(*)=0.0 
T2=l 

WE  WOULD  LIKE  To  DO  "DO  500  CH=1  .NROWS.TWSZR4'  BUT  CFD  INSISTS 

Km  J!<Ll^RE!!ENT  BE  A C0NSTANT»  SO  WE  WILL  CONSTRUCT  AN 
hUUI  VALENT  L(X)P. 

CH*I 

★DO  500  T7=I,NCHAN 
★IF (DEBUG. LT.2JG0  TO  420 
PINT! ( 1 )=CH 

\ OISPLH  "CH", 16, PINTI .PINTI I 

420  *CONTINUE 

★DO  430  T3=0,TWSZR-I 
INDEXI =CH+T3 

PRFAL1 (*)=ABS(BUFF2(*, INDEXI .OLD) -BUFF 2 (★- 1 , INDEXI -PINTI (★) .OLD)) 
PREAL2 (*)=ABS(BUFF2( ★-! , INDEXI -PINTI ( *) .OLD) -BUFF2( ★+! .INDEXI* 

I PINT2(*) .OLD)) 

★IF (DEBUG.LT. 2) GO  TO  422 
^ DISPLF  "PREAL1  ", I6.PREALI .PREALI +631 

^ OISPLF  "PRFAL2" , I 6.PREAL2 , PREAL2+631 
422  ★CONTINUE 

PHEAL2 ( ★ )=PREAL2 ( ★ ) *GLCHFT( ★ ) 

★ IF  (T3.  EQ  .0)M()DE=M()DE  .AND.  ( PEN(  ★)  .NE.I  ) 

★IF ( T3. EQ .TWSZR-1  ) MODE-MODE.  AND .(  PENT*)  .NE.64) 

*IF((  PREALI  (★)  . GT.PREAL2(*) ) ) BUFF2 (★,  INDEX i ,()LD)  = 

1 ( BUFF2 (★-! , INDEXI  -PINTI ( ★) ,0LD)+BUFF2(*+I , INDEXI +PINT2(*) , OLD) ) 

2 /2 . 0 

★IF (DEBUG. LT.2)G0  TO  430 
SETC(O)  El 
DISPLH  "MODE",  1 1 
PREALI  (*)  = BUFF2(  ★,  INDEXI  ,()LD) 


0;>  > onon 
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430 


440 


450 


C 

C 

C 

C 

c 


470 


500 


iC? 


DI SPLF  "B JFF2-" , 1 6.PREAI. I , PREALI ♦63l 

♦CONTINUE 
MODE=ON 

NOW  REMOVE  THE  BIAS. 

PREALI (♦)=0.0 
♦DO  440  T3=0,TWSZR-I 
INDEX I =T3+CH 

♦CONTniJE>SPREALI  U,*R()WSUM(BUFF2(a.INDEXI  .OLD)) 

PREALI (♦ )=PREALI (♦)/FLOAT(TWSZ) 

♦DO  450  T3=0, TWSZR-I 
INDEXI =T3*CH 

,ndexi  ... 

DI SPLF  "BIAS" . I 6 .PREALI , PREALI +631 
PREALI  (♦  ) = B1JFF2  (*, INDEXI  .OLD) 

♦CONTI N UE A"BIAS",,6*PREAL,»PREAL,+63» 

Now  LETS  COMPUTE  THE  MEAN  SQUARE  FOR  EACH  CHANNFL  ANn  cno  tuc 
ENTIRE  TIME  WINDOW.  PREALI  WILL CONTAIN IS SOWR^  Fn^ 

THE  PARTICULAR  CHANNEL.  ALLMSO  WILL  CONTAIN  THF  MEAN  SOUARF  Fna 

PREALI (*)=o.O 
♦DO  470  T3»0, TWSZR-I 
INDEXI »T3+CH 

♦CONTI NUE)“PHEALI (*,+R0WSUM(BUff2(*.INDEXI ,0LD>~2) 

PREALI  (*)  = PREALI  (♦)/FL()AT(TWSZ) 

CilMSQ!  T2  )»PREALI ( I ) 

I"2*T2«T 

ALLMSOI ♦)=ALLMSO( * ) ♦PREALI (♦) 

CH*CH*TWSZR 
♦ CONTINIIF 

sg 

CRITERIA  EASIER  AND  EASIER  UNTIL  AT  LEAST  HALF  THE  CHANNELS  PASS 
w TH  A LIMIT  OF  TEN  TIMFS  ON  OUR  PATIENCE.  CHANNELS  PASS. 

ALLMSQ (*) =ALLMSO( *) /FLOAT! NCHAN ) 

♦IFf DEBUG.LT.2)00  To  502 
D I SPLF  "ALLMSO", I 6, ALLMSO, ALLMS0*63l 
3 I SPLF  "CHMSQ" , I 6.CHMS0.CHMSQ+50I 
♦CONTINUE 
LT3«LNCHAN.SHR.I 
4 D'JM A*  AS<  STATEMENT 

thats  half  the  number 
TVARFT(a)»VARFT(a) 

PUT  VARFT  IN  THE  PE'S 
♦DO  550  TI-I.IO 
*IF(DEbUG.LT.2)C0  To  504 
DISPLH  "VARLOOP*,0l 


TO 

OF 


FORCE  DEALLOCATION  OF  REGISTERS. 
CHANNELS,  THE  NUMBER  THAT  NEED  To  PASS. 


SO  WE  CAN  GET  To  IT  EASILY. 


o o n o n n o 
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C 

C 

c 


550 

560 


504  ♦CONTINUE 
NCDCH-0 

*D()  540  CH=  I , NCHAN 
CHGOOD(CH)»0 
aREALi (*)»CHMSO(CH) 

PREALI (*)=PREALI ( ♦)/ALLMSQ( ♦ > 

^ F *i^‘/(!nSrA!-,U>,GT*TVARf:TU?>>>GO  TO  510 
IF  ( .ANY  . ( ( PREALI  (♦)  ,LT.  ( I .O/TVARFTf ★ miirn  t « tin 
CHGOOO(CH)  = I •u/ivakm(*)))))00  TO  510 

NGDCH= NGDCH+ 1 
510  ♦CONTINUE 
540  ♦CONTINUE 

*IF(NGDCH.GT.T3)G(>  To  560 
T»/ARFT(*)»TVARFT(*)*1  .25 
♦CONTINUE 
♦CONTINUE 

SOTlWcSl  pTea'cTa^Y0'^™’  "E  HW * SEPARATE 

.if??SraAyB!!nAI.i',,g,AItI';^av  21  n‘«sa»-array  3 

LASA  DATA.  CHANNELS  ARE  ARRANGED  VVVV  nmmm  trccc 

*IF(NANY  (TcoMpf*?  SEnM!J!TERESTlNG  CHANNELS.*  *# 

: 0MEANS*Je  G^-NE-0»>°»T»580 

♦DO  570  CH=NSITE+I,NCHAN 
CHGOOD(CH)*0 
570  ♦CONTINUE 
♦GO  TO  800 
oHO  +C0NTINUE 

d?sRJlh’^r«80»,'™  pb,icessed  here.  n„t  implemented  yet. 

♦GO  To  1 1 00 
600  ♦CONTINUE 

*IF(ARRAY.NE.2)GO  To  650 

J^^ATjONHWSTL0E  DWEN"EaLmCe5T|23?P3i™EE  AN°LES  AN°  S°“E 
■ T" 

: .,-S!:J;^-fSIEUNC  11 

I NDEX3= I , ( INC  TWSZR*3) 

I NDEX4= | , (INC  3) 

INDE X2= I 
INDEX3=I 
I NDEX4= I 

T I =TW5ZR+TWSZR+TWSZR 
♦DO  620  I NDEX 1 « I ,NSITE 
T2=0 

T6=CHG(X)D(  INDEX4) 


r> 

V./ 

A 


3 COMPONENTS  PER  SITE. 


* . i< 


o n o r>  o o o 


I 


<<E’ RR>f F J .DEM2I  70 
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PACE  HP 


★IF (T6 . EG .0)00  TO  610 
T6*CH'j(K)D(  INDEX4-H  ) 

★ IF  CT6.  EO.O  )G()  TO  610 
I7>*CHG(X)D(  INDEX4+2) 

★IF(T6.E0.0)GO  TO  610 
f2* ! 

C ALL  THREE  COMP  )NENTS  MUST  BE  G(X)D  FOR  A SITE  TO  BE  C(K)D. 

★00  60-j  T 3=0,TWSZR-I 
T4=I NUEX2+T3 
T5=  I NDEX3+T3 
T6= I NDFX  3+T3+TWSZR 
T7= I N0FX3+TWSZR+TWSZR 

BIJFF2  ( * ,T4  , OLD  )=  .5  773b*( BUFF2( *,T5 , NEW ) + BUFF2(*,T6, NEW)* 

I BUFF2 ( *,T7, NEW  ) ) 

60S  ★CONTI  OIF 

I NDEX2=1NDFX2*TWSZR 
610  CHGOO.D ( I NDFX  I )=T2 
INDEX3* INDEX3+TI 
I NDF  X4= I N0FX4  + 3 

620  ★continue 
★go  To  aoo 
630  ★CONTINUE 

C HORIZONTAL  MOTION  PROCESSES  HERE.  NOT  IMPLEMENTED  YET. 

A DISPLrl  " ARG630"  , 2l 
★GO  To  I I 00 
650  ★CONTINUE 

★IF (ARRAY .NE. 3) GO  70  700 

NORSAR  DATA.  IT  IS  ARRAMGED  VNEVNEVNE...  WE  DONT  HAVE  TO  DO 
ANY  COMPUTING,  JUST  REARRANGING. 

★ I F ( . ANY . ( (COMP ( ★) .NE .0 ) ) ) GO  TO  690 
0 MEANS  VERTICAL. 

I NDEX I = I ,2 ,3  ..  .NS1TE  (INC  1 ) 

I NDFX2*  I , UNC  TWS7R) 

I NDEX 3= I , (INC  TWSZR*3) 

INDEX4=| , (INC  3) 

I NDEX2= I 
INDEX  3=1 
I NDEX4= I 

T I =TWSZR+TWSZR+TW5ZR 
★DO  680  I NDEX 1=1 ,NSITF 
T2=0 

T7*CHG(X)D(  INDFX4) 

★ I F (T7 . EG  .0 )G()  TO  670 
F2=  I 

★DO  660  T3=0,TWSZR-I 
T4=INI)FX2+T3 
T5= I NDEX3+T  3 

3UFF2  ( *,T4  ,()LD  ) = BUFF2(  *,T5  ,NEW) 

660  *C()NTINUE 

I NDEX2= I NDEX2+TWSZR 
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670  CHG()OD(  INDEX  I )=T2 
INDEX3=INDEX3+I 
INDEX4-INDEX4+3 
680  ^CONTINUE 
★GO  TO  300 
690  ^CONTINUE 


C HORIZONTAL  MOTION.  NOT  IMPLEMENTED  YET. 

A DISPLH  "ARG690" ,24 

★GO  TO  1100 
700  *CONTINUE 

ALL  KNOWN  ARR/.YS  HAVE  BEEN  CHECKED  FOR. 

DISPLH  "ARG700" ,2| 

★GO  TO  1100 
800  ★CONTINUE 

WE  NOW  HAVE  ALL  THE  DATA  CONVERTED  TO  64  BIT  FLOATING  POINT 

RPSn/vPn^Mn6^?  Sulfur? * VARIANCE  CHECKED,  MOTION  COMPONENTS 
RE$i)LVED  ANn  BAD  CHANNELS  MARKED. CHGOOD(-)  IS  NOW  REALLY 
FTFFF-**  ^HE  FIRST  NSITE  ENTRIES  ARE  THE  ONLY  ONES  WE  ARE 

SETLNGO^E?nSSp  JNDICATE  WHICH  SITES  ARE  GOOD.  WEENOW 

»F^NGSSI  T°  THE  NUMBER  0E  GOOD  SITES,  NGDCH  TO  THE  NUMBER 

VFCTnp  «tNcDSI  Ff)R  VERTICAL  MOTION)  AND  THE 

PhS  ?*?  THAT  SITES(I)  INDICATES  WHICH 

PHYSICAL  SITE  LOGICAL  SITE  I REALLY  IS. 

NGDST=0 


A 


C 

C 

C 

C 

C 

C 

C 

c 

c 


★DO  820  INDEX  1 * 1 , NSITE 
T7=SITEGD( INDEX! ) 

★IF(T7.EQ.O)GO  TO  810 
NGDST-NGDST+I 
SITES(NGDST) -INDEX! 

810  ^CONTINUE 
820  ^CONTINUE 

NGDCH-NGDST 

★IF( .ANY . ( (COMP(*)  .EO. 1 ) > JLNGDCH-LNGDCH.SHL. I 
★IF (DEBUG.LT. 1 )G0  TO  830 
A DISPLH  "SITEGD",!8,SITEGD,SITEGD+70» 

A DISPLH  "SITES" , I 6,SITES,SITES+30l 

330  ★CONTINUE 


C 

c 

c 

c 


NOW  ITS  TIME  FOR  FFT.  FIRST  WE  HAVE  TOCONVERT  TO  32-BIT 
FLOATING  POINT,  SINCE  THATS  HOW  FFT  EXPECTS  THE  INPUT 
PI  NT  I ( ★) =TWSZR*NGDCH  * 

NGDRsPINT! ( I ) 

★DO  850  INDEX! = I ,NGDR 

★CALL  C64T32( BUFF2(*, INDEX  I ,OLD ) ) 

850  ★CONTINUE 

NOW  FFT.  NGDCH  GIVES  THE  NUMBER  OF  FFT'S  TO  DO. 

TWSZ  GIVES  t::e  size  of  each  FFT. 

STARTING  ADDRESS  IS  BUFF2( I , 1 ,OLD) . 

ALL  THIS  IS  PASSED  IN  COMMON  TO  RUNFFT. 

★CALL  RUNFFT 

NOW  TO  CONVERT  BACK  TO  64  BIT  FLOATING  POINT. 

★DO  870  INDEX  1 = 1 ,NGDR 


> O o O O O o 
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C 

C870 

C 

C 

C 

C 

c 

c 

c 

c 

c 


★CALL  C32T64( BUFF2(*, INDEX1  ,OLD)) 
★CONTINUE 

TO  BUFF3. 

BUFF 3 IS  (EACH  PE) * 


NOW  WE  GO 
FORMAT  OF 
WORD  I i 
W0RD2 » 
WORD  3 i 
WORD  28 i 


TWTIME 
NGDST 
SITES(-) 
DATA 


(1  WORD) 

(I  WORD) 
(25  WORDS) 


F ( LOFREQ ) (CHI . . . CH (NGDST ) ) 


F(HIFREQ)(CHI . .. CHI NGDST) 
BF3PE=BF3PE+ I 
I BUFF3( BF3PE, I )=TWTIME(I ) 

I BUFF3 ( BF3PE ,2 )=NGDST 
★DO  920  INDEX  1 = I , NGDST 
I BUFF3(BF3PE, I NDEXI +2 )=SI TES ( INDEX  1 ) 
920  ^CONTINUE 


★IF(TWSZ.EG.64)TI =6 
*IF(TWSZ .EQ . 1 28 )TI  = 7 


★IF(TWSZ.EQ.256)TI=8 
★ IF (TWSZ. EO .5 1 2 )TI  = 9 
INDEXI  *28 

LNGT=LNGDCH.SHL.TI 

T2=0 


★IF ( 0LD.EQ.2 )T2*35200 
T5*L0FREQ 
T6=HIFREQ 
★DO  960  F=T5 ,T6 

WE  WOULD  LIKE  TO  DO  "DO  950  CH=0,NGT-1 ,TWSZ"  BUT  CFD  INSISTS  THAT 
THE  INCREMENT  BE  CONSTANT  SO  WE  WILL  CONSTRUCT  AN  EQUIVALENT 
Loop  • 


CH«0 


★DO  950  T7= I , NGDCU 
I NDEX2=CH+F+T2 

blJFF3(  BF3PE,  INDEXI  ) = RBUFF2(  I NDEX2) 
I NDEXI “INDEXI + 1 
*IF( INDEXI .LE.640JGO  TO  930 
BUFF 3 IS  OVER  FLOWING. 

DISPLH  " ARG930" ,2  t 
★GO  TO  1 1 00 
930  ★CONTINUE 


CHCH+TWSZ 
950  ^CONTINUE 
960  ★CONTINUE 

★I F ( BF3PE .LT. 64) GO  TO  1000 
C HAVE  TO  WRITE  OUT  BUFF3, 

BF3PE=0 

★WRITE ( 6 A , BUFF3 ( 1,1), 0UTDM2 ( OPAGE ) , 40 ) 
★WAIT  64 
0PAGE*0PAGE+40 
★ IF(DE31JG.LT.  I ) GO  TO  965 
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r * 


A DISPLH"BUFF3"  ,2 1 

★IF (DEBUG.LT. 2) GO  TO  965 
A DISPLH  , I 6,BUFF3,BUFF3+1 02 3 I 

965  *C()NTINUE 

C ZERO  OU  BUFF3.  NOT  REALLY  REQUIRED,  BUT  USEFUL  ANYWAY. 

*00  970  INDEX1-1 ,640 
3UFF3( *, INDEX! >*0.0 

970  ★CONTINUE 
I 000*C()NTINUE 

*IF(.ANY.( (TOTSCN (*) .NE .FINSCNf *) ) ) )00  TO  100 
I I 00*C()NTINUE 

C GOING  TO  END  OF  JOB. FIRST  WRITE  OUT  REMNANTS  OF  BUFF 3. 
*WRITE( 64, BUFF3( 1,1 ) , 0UTDM2 ( OPAGE ) ,40 ) 

★WAIT  64 

★IFCDEBUG.LT. 2 ) GO  TO  1105 
A DISPLH  "EOJ" , I 8 , BUFF 3, RUFF 3+ 1 02 3 1 
I 1 05 ★CONTI NUE 
A JUMP  EFINALPRIMTl 

A ULKl 

A BFINALPRINT * * ! 

A DATA  ( ("★★★★★★★★") 8, ODOA * 1 6) 10,"DEM2  GOING  TO  END  OF  JOB", 
A ( ( " ★★★★★★★★" ) a , ODOA « I 6 ) I 0 1 

A EFINALPRINTin 

A DISPLS  , I 6, BFINALPRINT, EFINAI.PRINT-U 

A CLSDISP  DISP2I 

★CONTINUE 
★STOP 
★END 


oooooooo 
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FKCOMB 

WRITTEN  FOR  ILLIAC  BY 
PROGRAM  READS  IN  DATA 
AND  ARRANGED  WITH  ONE 
JETECTS  SEISMICEVENTS 


ANN  KERR  MAY  1974. 

THAT  HAS  BEEN  FFT'D 
TIME  WINDOW  PER  PE  AND 
BY  SEARCHING  A THREE  DIMENSIONAL 


REAL 


SPACE, ONE  DIMENSION  OF  FREQUENCY  AND  TWO  DIMENSIONS 
OF  WA/E  NUMBER. 

DECLARATIONS! 

★PE  INTEGER  I NBUF (*,640) , CNTRL! ★, 6) , NCHAN( *) , PINT!  (★) , OFFSET! ★) . 

LOCATE! ★) ,NPTS! *) ,C0UNT2(*) .C0UNT3! ★) ,LOC2D(*,25 ) , 
I.0C3D! ★,25)  ,TWTI  ME!  ★) , ADJF!  ★) 

POWER! *,25) , FMAX  (*,25),FKX(*,25),FKY(*,25),RINBUF(*,640), 
X ! *,25 ) , Y! ★, 25 ) , FFT! *,612) , KERNEL! *,25 ) , 
XCOORD!*),YCOORD!*),PREALI !★),  COSK!*) ,SINK(*) , 
COSDK! ★) ,SINDK(*) .REAMER!*) ,FPMAX(*) ,KXMAX(*), 

KYMAX!*) ,DELX! *) ,DELY(*) .KXSEP! *) .KYSEP!*) ,KSEP!*) , 

VEL!*) ,AZ!*) .SIGNAL!*) .FSTAT!*) .SUMSC!*)  . 

TEST!*) ,K!*)  .CHANAV! *) .TPOWER!*) .FREQ!*) 

ADKX! 4 ) , ADKY (4 ) , YPOl NT! 50) , YMAX(50) , DX(500> , DY(500) 

BEAM !*) ,TPOW!*> .DELTAK! *) ,P 

PREAL2!*)  .RPOWER! *,25)  , I POWER! *,25 ) .RTPOW!*)  .ITPOW!*) 
n M A X 

LOFREQ.HIFREQ, DEBUG, SM.TI  , T2, ARRAY, PAGE, I .N.MNCHAN, 
MNPTS, NPOINT, SWITCH, NFREQ.IGO, LINE, LINES, INDEX, IP, 
TWIN, SAM, IFRFQ.J.NFREQI , REFINE, IND.YTOP 
YPMI ,SIGN,NTIMES,LINEPI 
DELTX! 3000) ,DELTY ! 3000 ) ,DIST 
DELTAX, DELTA Y.KX.KY 

DKX,  LOWER,  UPPER,  L I NEP,  HDKX , BORDER  ,TW!)H 
DELTAF, RADIUS, ANGLE 


2 

*PE 

1 

2 

4 

5 

6 
7 

*PE 

★ PE 

★ PF 
★PE 
*CU 


REAL 
REAL 
REAL 
INTEGER 
INTEGER 


REAL 

REAL 

REAL 

REAL 


2 
3 

★PE 
★PE 
*CU 

★cu 

*CU  LOGICAL  MODE3.NMODE 

★ EXTERNAL  MA X.FNGR ID, REALE, I MG, GRID, CHECKR, OUTPUT 

★COMMON /MAINFK/INBUF,CNTRL,NCHAN, PI NTI , OFFSET, LOCATE, NPTS.C0UNT2, 

1 C0UNT3, LOC2D.LOC3D, POWER, FMAX, FKX.FKY.X.Y, KERNEL, 

2 XC(X)RD,  YCOORD,  PREALI  ,COSK, SINK, BEAM, 

3 TPOW, DELTAK, RPOWER, I POWER,  COSDK, SINDK, 

4 BEAMFR.FPMAX.KXMAX, KYMAX, DELX.DELY, KXSEP, KYSEP, KSEP 

b TWTIME, TPOWER, VEL, AZ, SIGNAL, FSTAT, SUMSQ, TEST, K, 

6 CHANAV, FRFQ, ADJF, DX, DY  ,P .YPOINT, YMAX, ADKX, ADKY, 

7 KX  KY  DELTX  DFLTY  ’ 

★EQUIVALENCE!  INBIJF!  I ,1  ),RINRUF(  I ,1  ) ) . ! INBUF!  I ,28),  FFT!  1,1)) 
★EQUIVALENCE  ! I .LOFREQ) ,!  2, HI FREQ) , ! 3 .DEBUG) , ( 4 ,SM) , ( 5 ,TI ) , !6,T2) , 

* (7, ARRAY), !8, PAGE), (9,1), (10, N),  (I2.MNCHAN), 

2 (I  3, MNPTS), (I  4, NPOINT),! 15, SWITCH), (I  6, ICO), 

j ( 17,  INDEX), (18, IP), (19, DKX), 

4 (20, LOWER), (2  I, UPPER), (22, LINE), (23, LINES), 

5 (24, HDKX) ,(25,  BORDER) ,(26, TWOH) , 

6 (27, DELTAF), (28, RADIUS),  (29, SIGN), 

7 (30.M0DE3) , ( 31  ,NMODE ), (32 .TWIN) ,( 33, ANGLE) ,( 34, SAM) , 

★* 

8 ( 35.NFREQ) , ( 36, 1 FREQ ) , ( 37, J ) , ( 38 ,NFREQ1 ) , ( 39, REFINE) 

9 ,(40, NTI MES) , (41 , IND) , (42.YT0P) , (43, YPMI ) , ( 44, LI NEP) , 

0 (45.LINEPI) 

★DISK  AREA  CONPRM ( 1 ) , STCORD (1 ) , FK I N ( 8 1 ) 

MODE-ON 


*;:r 
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OPNDISP  FKDISPt 

JUMP  EHEADl 

FKD1  iP « AREA  "FKDISP"! 

BHEADi BLK l 

DATA  ( ( " **♦★****•' ) 8 , ODOA  • I A)  2, 
"START  EXECUTION  FKCOMB" , ( ( "♦*****, 
EilEADiDISPLS,  I 6,BHEAD,EHEAD-I I 
♦READ!  64,  I NBUF ( 1,1)  ,C()NPRM!  I ) , I ) 
♦WAIT  64 


>">8,0D0Ai  I 6)  2 I 


TWIN=INPUF!2, 1 ) 

DKX=RI NBUF  <12, 1 ) 

LOFREQ* I NBUF (6,1 ) 

HIFREO* I NBUF ( 7, I ) 

L()WER=RINBUF!  I 3,  I ) 

UPPEP-RI NBUF ( I 4,1) 

PREALI  ( *)=R I NJ'JF  ( IS,  I ) *0.01  74533 
ANGLE=PREALI  < I ) 

REFI NE=I NBUF ( I 6,  I ) 

SAM= I NBUF ( 17,1) 

NFREOHI FREQ-LOFREQ+1 
NFREGI *NFREO-l 
JEBUG=INBUF( 1,1) 

SKIP  , ETl 
T«WDS  I i 
ET»SLIT(0)  T l 
STORE! 0)  TWINl 
DISPLH  "TWIN", I 6,T,Tl 
STORE! 0)  DKX « 

DISPl.F  "DKX" ,1 6,T,T| 

STORE! 0)  LoFREOl 
DISPLH  "I.0FREQ"  , I 6,T,Tl 
STORF(O)  HIFREOl 
DISPLH  "HIFREO" , I 6,T,T» 

STORE! 0)  LOWER  I 
DISPi.F  "LOWER",  16, T,T| 

STORE! 0)  IJPPERl 
DISPLF  " UPPER" , I 6,T,Tl 
STORE! 0)  ANGLE  I 
DISPLH  "ANGI  F", I 6,T,Tl 
STORE! 0)  RFFINFl 
DISPLH  "REFINE", 16, T,Tl 
STORE! 0)  SAMi 
DISPLH  "SAM", I 6,T,Tt 
STORE! 0)  NFREOl 
DISPLH  "NFREQ" , I 6,T,T| 

STORF! 0)  DEBUG » 

DISPLH  "DEBUG" , I 6,T,T* 

♦READ! 64 , I NBUF! 1,1 ) .STCORD! I ) , I ) 
♦WAI T 64 

XC00RU!*)*RINBUF!*, I ) 

/COORD ! * )“RI NBUF! * ,2 ) 


I 
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JISPLF  "XC(X)RD",  I6,XC(X)RD,XC(X)RD+I6| 
DISPLF  "YC(X)RD", 1 6,YC(X)RD,  YCOORD+I  6» 
*READ(64, INBUF( 1,1) ,FKIN(  I ) , I ) 

♦WAIT  64  * 

r I = I NJUF ( 1,1) 

*00  10  ARRAY=I ,7 
*IF(ARRAY.EQ.  7)G()  To  9000 
THAT  MEANS  UNKNOWN  HEADFR 
r2=CNTRL( I .ARRAY) 

*IF(TI  ,E0.T2)G()  TO  15 

that  meant  we  found  it 
I 0*C()NTINUE 
I b*C()NTlNUE 
SLIT(O)  Tl 
STORECO)  ARRAY t 
DISPLF  "ARRAY  N()",I6,T.T| 

PAGE=2 

A BUNCH  OF  DEBUG  PRINT  OUT  0 F INITIAL  VALUES 
?0*CONTINIJE  L 

50*C()NTIMUE 
MODE-ON 

♦READ (64,  INRIJF(  I , I ) ,FKIN(  PAGE)  ,40) 

♦WAIT  64 
PAGE=PAGE+40 
T'WTI  ME  ( *)  = INBUF ( *,  I ) 

NCHAN( *)  =1 NBUF ( +,2) 

PI  NT  I (*)=MAX(NCHAN(*)) 

MNCHAN=P INTI ( I ) 

*IF<.ALL.(<NCHANU).EQ.O)>)  GO  TO  9100 
P I NT  I ( ♦ ) = 1 

♦DO  60  Tl = I ,25 

MODE* ( I NBUF (♦, PI  NT I ( ♦) +2 ) . EQ.TI ) 

X(  *, PI NTI  (♦)  ) = XC(X)RD(TI ) 

Y ( * , PI  NTI  (♦)  ) = YCOORD(TI  ) 

PINTI  ( ♦)  = PINTI ( ♦ ) + 1 

MODE=ON 

6OC0NTINUE 

♦IF (DEQUG.LT. I )G0  TO  70 
DISPLH  "C(X)RDIN»tO? 

♦IF (DEQUG.LT. 2) GO  TO  70 
QISPLF  "X-0" , I 6,X,X+3| 

OISPLF  " X— I ",  I 6, X+64 , X+64+ 3l 
DISPLF  "X-2" , I 6, X v2+64, X+2*64+3| 

PffPLF  "X-3" , I 6, X+3*64,X+3*64+3l 
DISPLF  "Y-0", l6,X,X+3? 

DISPLF  "Y-l " , I 6,Y+64 , Y+64+31 
DISPLF  "Y-2" , I 6, Y+2*64, Y+2*64+3? 

DISPLF  "Y-3" , I 6, Y+3*64,Y+3*64+3» 
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\ ISPL.li  "NCHAN",I  A,NCHAN,NCHAN+3l 

70*C()NTIi4UE 
OFFSET !*>-0 

* ).)  I0:)0  iFREQ=I.OFKEQ,HIFREQ 
SWITCH  =1 
♦CALL  GRID 
rPo:vFrt(*)=o.o 
RTPOW!  * > = o . o 
I TPOW!  *>=0.0 
LOCATE (*)=  I 
*00  I 00  N=  I ,MNCHAN 
"ODF=<  N.LE.NCHAN ( * ) ) 

*:ALL  ,<F  ALE!FFT!*,OFFSET!*)*N),PkEALl  <*)) 

♦CALL  I VO  ! FFT !*, OFFSET! *)+N) ,PREAL2!*) ) 

, j ^ WN  EL  C * , N ) = ♦ ( A . 2 8 3 1 85.30*  ( DELTX  ( I )*X(*f  NJ^DELTY!  I )*Yf  *,f. 

C0SK!*)=C05! KERNEL! *,N)> 

51  UK ! * )=SIN! KERNEL! *,N  ) ) 

i<Po,iF,;  ( * ,N  ) = PREAL  I ! ★)*C0SK'*)-PREAL2!*)*SINK!*> 

I F : ...T ,<  ( * , N ) = PREAL I !*)*SI NK! * ) + PREAL2 ! * ) *C()SK(*) 

RTPOW!  ★) =RTP0W ! * ) ♦RPoWER  ! *,  N ) 

I TPOW!  ★)= ITPOW ! * ) ♦ I PowFR ! *.N ) 

1 00*C0NTI  JfJE 
m;)DE=on 

. P()f'>  < ★ ) = I TPOW  ( ★)  **2+RTP()W!  *)  *+2 
* I F ! DEb  UG . LT. 2 ) GO  TO  110 
OISPLF  "FIRST", 0| 

)ISPLF  " RTPOW" , I A,kTP0W,kTP0W>3l 
D ISP  LI-  " I TPOW" , I A,  ITPOW,  ITP()W»3| 

'PLF  " TPOWFR" , I A,TPONFR,TPOWER^3» 

ISPLF  "DFLTX! I ) " , I A , DELTX , DFLTX t 
'I  5PLF  “ UELTY ( I ) " , I 6,  LFLTY,!)FLTY| 

I 10  ★CONTINUE  * 

RESTORE  MODF 
ri=‘JPTS!  I ) 

*)0  300  NP0INT-2.TI 
RTPOW'! *)=0.0 
ITPOW! *)«0.0 
M3 1 T I = MODE 
ISO  *0(>NTIMUF 

★1)0  2 00  N=  I ,MNCHAM 
)ODF»MODE. AND. ! N.LE.NCHAN!*) ) 

)ELTA\  (*)  =+A.283l  HS30*! OELTX ! NP()I NT) *X! *.N ) ♦ 

I DELTY ! NPOl NT)*Y!  *,N ) ) 

CoSDK! *)»COS !DELTAK! * ) ) 

SINDK! *)=SIN!DELTAK( *) ) 

PkEALI ( * >=RPONEk! *,N )*COSDK! *>-I POWER!  *,N>*SINr)K!*) 

1 ER  ( * ' jj  } =pPRA^/ 1 * * fi  > *S  1 Nn*  ( * > * I POWER ! *,  N ) *COSDK ! * ) 

RTPOW!  *)=RTP!)W!  *)  ♦RPOWER  !*,N) 

I TPOW!  *)  = ITPOW  ♦ I POWER  !*,N) 
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200*C0NTINUE 

MODE-MBITI 

TPOW!*)=RTPOW!*)**2+ITPOW!*)**2 


RESTORE  MODE 

MODE«!TP!)W!*).GT.TPOWER!* 
TPOWER  ( ★ )-TPOW  ( * ) 

LOCATE ( ★ )=NP0I NT 
MODE-MBITI 
3O0*CONTINUE 

FPMAX!*)-TPOWER!*) 


)) 


A 

A 

A 

A 

A 


A 

A 

A 


MODE-ON 
SWITCH-2 
★CALL  GRID 
★DO  350  N=1 ,MNCHAN 
MODE-IN.LE.NCHAN!*)) 

★CALL  H EALE !FFT!*, OFFSET! *)+N) ,PREALI  (★)) 

★CALL  IMG  !FFT!*, OFFSET! *)+N),PREAL2!*)) 

KERNEL! *,N)«  ♦ ! 6. 283 1 8530* ! KXMAX! *)*X! *, 

1 Y!*,N))) 

COSK!  *)=C()S!  KERNEL!*  ,N  ) ) 

S INK !* )-S IN! KERNEL! *,N ) ) 

RPOWER !*,N)-PREAL1 !*)*C0SK!*)-PREAL2!*)*SINK! ★> 

I POWER ! *,N )-PREAL1 !*)*SINK!*)+PREAL2!*)*C0SK!*) 
350*C0NTINUE 
MODE-ON 

★IF! DEBUG.LT. 2) GO  TO  360 
DISPLF  "COARSE", 01 
DISPLH  "NPTS" , I 6.NPTS.NPTSI 
DISPLH  "LOCATE",  1 6,  LOCATE,  LOCATED  I 
DISPLF  "KXMAX" , 1 6,KXMAX,KXMAX+3» 

DISPLF  "KYMAX" , I 6.KYMAX.KYMAX+3* 

PREAL1 !*)-RPOWER!*,N) 

PREAL2!*)=IP0WER!*,N) 

DISPLF  "RPOWER", I6.PREALI ,PREAL1  +3l 
DISPLF  "IPOWER", I6,PREAL2,PREAL2+3| 

DISPLF  "FPMAX" , 1 6,FPMAX,FPMAX+3l 


NJ+KYMAX! 


*)★ 


360  ^CONTINUE 
★CALL  FNGRID 
FMAX!*, I FREQ) -FPMAX! ★) 

FKX!*, I FREQ)® KXMAX! ★) 

FKY!*, I FREQ) -KYMAX! ★) 

OFFSET! ★)«OFFSET!*)+NCHAN!*) 


I000*C0NTINUE 
★CALL  CHECKR 

C DISPLAY  MAXIMUM  AND  ASSORTED  PARAMETERS 
★CALL  OUTPUT 


★GO  TO  50 
9000*C0NTINUE 
A DISPLH  "BAD  HE AD", 0 I 

91  00*C()NTINUE 
A JUMP  PBANNERl 


23< 
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A BANNEiiiBLKl 

• )A f A (("♦♦♦*♦*♦♦'• J3.0C0AI I 6)2,"FKCOMB  EOJ", 
A ((»*★******••  >8t0D0Al16»2* 

A I [JISPLS  , 1 6, BANNER, PBANNER-I  I 
A .LSD ISP  FKQISPl 

♦STOP 

♦END 


MM* 
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★BLOCK  DATA 

★PE  INTEGER  I NBUFU, 640)  , CNTRL( ★, 6)  , NCHAN(  ★) , PI NTI  ( ★)  .OFFSET!  ★ ) , 

1 LOCATE! ★)  ,NPTS( *)  ,C0UNT2(*)  ,C0UNT3!*)  ,L()C2D(*,25) , 

2 LOC3D(*t25) ,TWTIME!*),ADJF<*) 

★ PE  REAL  POWER! *,25 ) , FMAX!  *,25 ) ,FKX! * ,25 ) ,FKY(  *,25  ) ,RINRUF(  ★,  MO)  , - 

1 X!*,25),Y!*,25),FFT!*,612), KERNEL ! *,25 ) , 

2 XC(X)RD(*)  ,YC(X)RD!  *) , PREALI  ! *) , C!)SK ! *) ,SINK( *) , 

4 COSDK!  ★)  ,S INDK! *)  ,BEAMER(* ) , FPMAX! *) , KXMAX! *) , 

5 KYMAX! *)  ,DELX!*)  ,DELY!*)  ,KXSEP(*) .KYSEP! *), KSEP! *) , 

6 VEL! *)  ,AZ!*) .SIGNAL!*) ,FSTAT(*) ,SUMSQ! ★)  , 

7 TEST!*)  ,K(*)  .CHANAV!*)  ,TP()WER(*)  ,FREQ!*) 

★PE  REAL  ADKX!4)  ,ADKY!4)  , YP()INT(50)  ,YMAX(50)  ,DX!500)  ,DY!500) 

★PE  REAL  BEAM(*),TPOW!*) ,DELTAK<*)  ,P 

★PE  REAL  PREAL2! ★ ) .RPOWER! *,25 ) , I POWER!  ★, 25  ) , RTPOW ! *) , ITPOW! *) 

★PE  INTEGER  MAX 

★PE  REAL  DELTX!  3000) ,DELTY!  3000)  ,DI  ST 
★PE  REAL  DELTAX, DELTA Y,KX,KY 

★COMM()N/MAINFK/INBUF,CNTRL,NCHAM,PINTI , OFFSET, LOCATE, NPTS.C0UNT2,  - 

1 C0UNT3, LOC2D,L(£3D, POWER, FMA X, FK X, FKY ,X,Y, KERNEL, 

2 XC!X)RD,  YC!X)RD,  PREALI  ,COSK, SINK, BEAM, 

3 TPOW,DELTAK, RPOWER,  I POWER , COSDK,SINDK,  - 

4 BEAMER,  FPMAX,  KXMAX,  KYMAX, DELX, DELY, KXSEP,  KYSEP,  KSFI\- 

5 TWTIME.TPOWER, VEL, AZ, SIGNAL, FSTAT.SUMSQ, TEST, K, 

6 CHANAV,FREQ,ADJF,DX,DY,P,YPOlNT,YMAX,ADKX,ADKY, 

7 KX  KY  DELTX  DELTY 

★EQUIVALENCE! INBUF!  I , I ) , RINBIJF!  I , I ) ), ! INBUF!  1 ,2R)  , FFT!  1,1)) 

★DATA  CNTRL/54227, 63*0,4961 9,63*0,60366,63*0, I 92*0/ 

★END 
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♦BLOCK  DATA 

♦PE  INTEGER  CNTRL(*,6> ,OUTBUF<*, 64,6) tPINTl <*>, INBUF'*, 128) . 

I TIME!*) , ()LDTIM(*) , 

1 SAVBCT,SAVPTN,OUPAGE<6>,TSTEPS(6> .SCANS, 

2 0UPTNA<6>,  OTIMEAC6)  ,ORGADR,  INRIJF I (8192) 

♦ COMmon/ma  I N/CNTRL#OUTBUF  .INBUF,  PI  NT  I ,TIME,OLDTIM,SAVBCT,SAVPTW, 

I ESTEPS .SCANS , OUPTWA  , OUPAGE , OTI NEA , ORGADB 

♦ EQUIVALENCE  ( INBIIFI I , I ) ( INBUF  I ( I ) ) 

♦DATA  UUPAGE/6*2/ 

♦DATA  CNTRL/I ,1 ,51,16,1,0,1 03,446, 1 0, I ,54227,53*0. 

1 G, 2, 57, 0,0, 6,29, 95, 15,2,4961 9,53*0, 

2 0, 1 ,66,0, 1 ,72,12,0,10,1  ,59366,53*0, 

3 10*0,65536,53*0,10*0,65536,53*0,10*0,65536,53*0/ 


f 


* UI  rvj 
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★BLOCK  DATA 

r ,nteoer  ii: : it i „ - 

*PE  INTEGER  .ABUFF2<704ro, .- 

CHMS0(80).RBUFFI  (4096)  ,R0WSUM,RBUFF2(70400) 
★COMMON/MAIN2/NBUFFI  ,BUFF2,BUFF3,FI NSCN,C()MP,TOTSCN,PINTl  PIMT2 
p » PEN , CNTRL ,GLCHFT , VARFT, PRFAL I .PREAI  2 * 

*E0U!VALE^E^?NBUFF^n*^T^RBUFFE(T?'?BUFFU?)T^?iiuFF2?l  ^1  _ 

★DATA  PEN/I«2*3, 4.5,6, 7,8, 9, 10, II, 12,13, 14, 15, 16, 17, 18, 19, 20, 

2 po’in*??*^'2^26'27'28*29*30,31  *32»  33,34,35,36,37,38, 

3 - 

★DATA  CNTRL/54227 ,0,51 ,1 7,60*0, 

1 49619,0,57,19, 60*0, 

60366,0,66,22,60*0, 

END  65536, 63*0, 65536, 63*0, 65536, 63*0/ 


Best  Available 

Copy 

for  page  29 
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★SUBROUTINE  CHECK R 

W'JRP()>E*  VERIFY  Tif.':)  AML  THREE  DIVENSJ  ONA'  HX I .T'  - v"  • 

"A  XI HI’.'-  r;'!R  i.'AC! ; FREQUENCY. 

CL ARATI ON. 5 i 

INTEGER  IN'JIJF(*,640)  ,CNTRL<*,6)  ,NCHAN(  *)  ,p|TTl  (*)  , . ’ 

LOCATE!*)  , N PTS ( * ) .COUNT2!*)  ,CnUNT3(*)  ,1  :iC2i  , 

LOC  3D ( * , 25 ) , TWTI M E ( * ) , ADJF ( * ) , I MDF X ( * ) 

* AI  POWER!  * , 25 ) , FMAX  ! * *25  > , FKX ! * , 25 ) , FKY  < * , 25  ) , R I NBUF  ( * <5 

I X(*,25),Y(*,25>,FFT(*,*I2) .KERNEL!*, 25) , 

XC(X)RI ) ! * ) , YC (X)RD( * ) , PREA L I ( * ) , K!*)  (*) 

COSIJ.L!*)  ,GP'DK(*>  .REAMER!*)  ,F!>MAX(  * ) , -v  < /. y ' - 
K YMAX(  * ) :!.X  ^ *)  , OELY!  *)  ,KXSED ! * ) , CYC:  '(*),  ; * ; , 

/E!.(  *)  , AZ(  *)  , IGNAL!  *)  .FSTAT ( *)  ,SB '(5Q{ *) 

/ TEST!*) ,K(*) .CHANAV! *) .TPOWER! *) ,'  RFC!*) 

* REAL  ADKXT4 ) ,ADKY! 4) ,YP0INT!50) ,YMAX(50) .DXL500)  ,DY(500) 

* \L  V.  ( * ) ,TC ::•(*) , UELTAK  ( * ) ,C!  :!•  .■;;*,  ; . , , 

*■  EA1  PREAI  2!*)  ,RPOWER(*,25),lPOWER(*,25),RTPOW(*)  , i TPO*  { * > 

» ' I.  :r  "AX 

* • • Lg;r;  u.i  IF  n ri'-C , DE3UG.SM , TI  ...  . ■ . 

‘ NPT5 . NPO I NT,  SN I TCH , NFREO , I ' ; , I’.-. 
TWIN.SAM.IFREQ.J.NFRFQl .REFIN  - , 

Y2"  I .SIGN, MI  MES ,L  I NI  P I 
’ :/(  < >J),  ■ LTY! 3100),  r . 

. LTAX, DELTA Y,KY,  Y 

.I-INEP.H  •, 

• - I A , i i , A‘-:JLF: 

* ' ' L!  il  -Aj  M();)E3,0"()Db 

* ,i  , /• : . i ,1  ■(;,  i.. 

*  I/MAINFK/INPUF.CNTRL.NCHAN, PI NT)  .OFFS  )CA  o 

. i,LDC2Q,LO(:3:;,  .■  ■ ■. 

XC.  • O.YC:  : 

1 rPO  i.DELTAK.RPOKER,  IPO-.V-  /, 

; 

i 1 , ' * A.  , . ■ < , s i 

CM  A )AV  ,Ff:'"'i,ADJF  ,D> , )Y,  ,Y 

7 KX  KY 

' ■ n.i >,  i , i . 

< I *L  .hi f . . , 

( 7,  A Ri : A Y ) , ! f , -A  If  >,(■/,  i ),(  I , . i , 

( I 3,  '•"G*T5)  , f ! 4 , 'IPoI  NT)  , ! I 5 , - ..‘i  : . 

3 ( I R I P ) ( i 9 )KX ) 

(20, LOWER), (21, UPPER), (22, 1.IN  >,  . / 

(24,t|DKX),!25,  iOR  , 

A (27, DELTAF), (28, RADIUS),  <20, SIGN), 

! 3'  , ' > 3)  , r 31  , I )')!- ) , (32,  ' 

+ ★ 

( 35  .NFRFQ ) , ( 36,  IFRFO) , ( 37,  I) , ( 3E , 

,(40,NTP.'FS)  ,(4|  ,IND),(42,YT  :i  ,(\ 

( 45 .LINFP I ) 

* ISK  ARf; A CONPRM!  I ),STCORD(  I ),FKIN(HI  ) 

* I E ( i/T  UG.LT.  I )C()  TO  10 


2b* 
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A DISPLH  "CHECKR” ,0| 

10  *C0NTINUE 
C0UNT2(*)=0 
C0IJNT3(*)=0 
PREALI  ( *)=0. 5*PKX 
HDKX=PREALI ( I ) 

★DO  100  I=L()FREQ+I  ,HIFREQ-I 
M0DE3* ( FMAX ( ★,  I)  ,GT.FMAX(*,I-I )) 

MBITI* (FMAX( *,  I ) .GT.FMAX (★, I +1  )) 

MODE 3-MODE3. AND. MBITI 
MODE* MODE .AND.  .NOT. M0DE3 

★ IF ( ( FMAX(  ★ , I -I ) .GT.FMAXf  ★ , I ) ) ) INDEX ( *)  = I -1 

★ I F ( ( FMAX(*, I + l ) ,GT.FMAX(*,I  )))  INDEX(*)=I+I 
KXSEP(*)=FKX(*,I )-FKX( ★, I NDEX( ★) ) 

KYSEPC *)=FKY (★,  I )-FKY( ★, INDEX!  ★ ) ) 
KSEP(*)=SORT(KXSEP<*)**2+KYSEP(*>**2) 

NM()DE=  ( KSEP( *)  .LT.HDKX) 

M()DE=M()DE  .AND.  .NOT.NMODE 
AD JF ( ★ ) = I — I 
★DO  560  J=l  ,2 

DELY  ( * ) =FKY  ( ★ , AD  JF  ( * ) ) -FKY ( ★ , I ) 

DELX (*)=FKX(  *, ADJF(*  ) ) -FKXC* , I ) 

MBIT2=M0DE 

OFFSET  (*)  = ( ADJF  (*)-L()F  REQ)  *NCH  AN  ( ★ ) 

HTPOWUI-0.0 
ITPOWt *)=0.0 
★DO  550  N=  I .MNCHAN 
MODE-M ODE . AND . ( N . LE . NCHA N ( ★ ) ) 

★CALL  REALE ( FFT ( * , OFFSET ( ★ ) +N ) , PREAL I (★)) 

★CALL  IMG  (FFTU,0FFSET(*)+N)fPREAL2(*)) 

DELTAK(*)=  +< 6.2831 8530*(FKX(*f ADJF ( ★ ) )*X (★, N) + 

I FKY ( *,ADJF(*) )*Y(**N ) ) ) 

COSK(*)=COS(DELTAK(*)) 

5 1 NK ( ★ )=S I N ( DELTAK  ( ★ ) ) 

3TP0W( *)  = PREAL I ( ★)*C0SK(*)-PREAL2  ( ★)  *S  I NK(  *) +RTPOWC ★ ) 

I TPOW(  *)  = PREAL  I ( *)*S  INK(*)  +PREAL2C  ★)  ★C()SK(  ★)  + ITPOWC  ★ ) 
550*COHTINUE 
MODE' MB I T2 

CHECKC ★, J ) =RTP()W  (★)  ★*2  + ITP()W(  ★)  **2 
ADJF (*)=I + I 
560*CC)NTINUE 
MODE* ON 

MBIT2* (FMAX( ★, I ) .GT.CHECK(*f « ) ) 

MBITI* (FMAX (★,  I)  .GT.CHECK(*,2)) 

MBITI -MBITI . AND.MBIT2 
MBITI =MBITI .AND.. NOT.NMODE 
MODE3-MODE3. OR. MBITI 
NMODE*  .NOT  .M0DE3 
MODE-HMODE 

C0UNT2 (*)-C  UNT2( *)+l 
L0C2l)(  C0UNT2  (★))*! 

MODE* MODE 3 

COUNT3(*)=COUNT3(*)  + l 
L0C3D(  ★,C()UNT3(  ★ ) ) = I 
MODE-ON 
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A SETC(O)  G| 

A 5 ETC ( 1 ) Hi 

A IJISPLH  ''MODES'',  I | 

A iJlSPLF  "CHFCKI  ",  I 6,CHECK ,CHECK+9| 

A DISPLF  "CHECK2", I 6,CHECK+64,CHECK+73l 

I OO*  CONTINUE 
* PE TURN 
★END 
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A 

A 

A 

A 

C 

C 

A 

A 

A 

A 


★PE 


*PE 

I 

★PE 

I 

★PE 

*CU 


2 

3 

*CU 

1 

2 
3 


★SUBROUTINE  CI6T64( IN, OUT) 

★PE  REAL  IN (★) .OUT ( ★) 

INTEGER  NBUFFI  ( ★, 64 ) .FINSCN (★ ) ,COMP(  ★)  .TOTSCNU)  ,PINTI  (★) 

PINT2 (*),TIME(*),OTIME(*) ,TWTI ME( ★) , PEN (★) ,CNTRL( ★ 6) 

TOE!S^Il:X^Sci>®T^RFT?5?0•2,•^UFF3C*•MO>•PREAI-,^*,•  - 

REAL  CHMSQ(80)  , RBUFFI  (4096)  ,RC)WStJM,  RBUFF2  ( 70400) 

INTEGER  A[}BRUF( |fi) i.COREPT,  BYTE,  ADBWRD,  ARRAY,  DEBUG,  TWSZ  , 

TNnpyf,TwnryoNfiJcvoR?WS,DIFFR,DrFFW’NEW’oLD»GAP*TSCAN5«- 
INDEXI  , INDEX2, INDEX3, INDEX4 ,TI  ,T2 ,T3 ,T4 ,T5 .T6.CH 

r , °fFSET.INBYT,NGDCH,TWSZR,NGDST, NODR, F.BF3PE.NGT. PAGE  T7 

LOGICAL  LADBBU(8) ,IC0REP,LASTI6,LBYTE,LADBWR,LARRAY,LDERUG, ’ 

-™SZ . E(A  nLA  £ • E1?CH • LNS 1 TE * LNROWS , LDI  FFR ,LDI  FFW , LNF  W , - 

LOLD, LGAP ,LTSCAN  ,LTI  ,LT2,LT3,LT4,LT5,LT6,LCH,L0FFSF 

LINBYT,LF,LNGDCH,LTWSZR,LNGDST,LNGDR.LNOT  LT7 

★EXTERNAL  GTDATA ,C64T32 ,ROWSUM, RUNFFT.C32T64  * 

★C0MM0N/MAIN2 /NBUFFI  ,BUFF2,BUFF3,FINSCN,C0MP,T()TSCN,PINTI  ,PINT2 

1 TI MF,()TIME,TWTIME,PEN,CNTRL,GLCHFT,  VARFT.PREALI  PRF.AL2  * 

2 ALLMSQ,TVARFT,LOFRFQ,HIFREQ,S  ITEGD, SITES. CHMSQ  ’ ’ 

★EQUIVALENCE  (NBUFFI (1,1)  ,RBUFF 1(1) , I BUFF  I ( I ) ) , (BUFF2(  1,1,1  ) 

1 ABUFF2(  I i ,RBUFF2(  I ))  ,(CHGOOD(  I ) ,SITEGD(  I )) , * 

2 (BUFF3(  I , I ),I BIJFF3(  1,1))  * 

★EQUIVALENCE  ( I ,ADBBUF( I ) ,LADFBU( I ) ) , (9,C0REPT,LC0REP) , ( I 0, BYTE 

I~Tw^7 ! * J •’ i * ^"wP^^R)  * ^ 1 ^ * ARRAY , LA  RRAY ) , ( I 3, TWSZ,  - 

( 1 ? «()yf-Ap«LOVLAP) , ( I 5, NCHAN.LNCHAN) , ( I 6.NSITF,  - 
1 7.NR0WS.LNR0WS) ,( I8,DIFFR,LDIFFR) ,(  I9.DIFFW,  - 
♦()Ld.LOLD),(22.GAP.LGAP>,(23,  - 

IngfK  w } ♦ (25»  INDEX2) , (26,  INDEX3)  ,(27.- 

INDEX4) ,(28,T! ,LT1 > , ( 29 ,T2 ,LT2 ) , ( 30, T3, LT3) , ( 31  ,T4, 

lacier?2?15*1  T5)f  (3’3*T6,LT6)  ♦ ( 34.CH.LCH)  ,(35, OFFSET, 

V^?^^3r^5^T-LINBYT)-(37-F-LF)^38-NGDCH.LNGBcH),  - 
^M?on2:L7>,SZR)  ’(40*NGDST’LNG0ST)  •(4t  ♦ngdr.lngdr)  , - 

?^V-!Grr4-^GT),(44,LAST,6)  ’(45*DERUG»LDEBnG)  •- 

^ rAuhJ  1 1 4 7 f T7 f I.T7 ) 

INDM2 ( 20 ) , 0IJTDM2 (40)  ,C()NPRM(  I ) 

T < i 


2 

3 

4 

5 

6 

7 

8 
9 
0 
I 

★DISK  AREA 

★IF (DEBUC.LT. I )G0  TO 
IJISPLH  "Cl  6T64"  , 0| 
★IF(DEBUG.LT.3)G0  TO 
LDL(O)  SD49t 
LDA  IN (0) t 


10 


DISPLH  "IN" , 32 1 
10  ★CONTINUE 

NOW  TO  CHECK  WHICH  ARRAY  WE  HAVF. 

★ I F ( A RRAY . GT. 2 ) GO  TO  100 

LASA  AND  ALPA  HERE.  SIMPLE  14  BIT  TWOS  COMPLEMENT. 
LDL(O)  $D49| 

LDA  IN (0 ) I 

SHAR  =2»  % RIGHT  JUSTIFY  IT. 

AND  =3FFFH6l  % GOT  RID  OF  ANY  GARBAGE  BITS. 


* .■« 


< 
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LOU  SA l % SAVE  IT 

X NOW  TO  TAKE  CARE  OF  THE  SIGN. 
LIT(O)  = OCOOE  OOOOOOOOOOOO  i 1 6 1 

% WITH  A I 


THE  HIGH  ORDER  BIT  ()T  THE 
^TS  -(2**13)  WHICH  IS  THE 


SHAR  = I 3 1 
SHAL  =47 l 

LEX  SCOl 
NORM  | 

SANl 
LOS  $ A I 
LDA  SRI 
SHAL=35 l 


% MANTISSA  THATS  -(2**13)  WHICH  IS  Tl 
% VALUE  OF  THE  SIGN  BIT. 
x ISOLATE  THE  SIGN  BIT. 

XP'JT  IT  IN  THF  HIGH  ORDER  BIT  OF  THE 
% MANTISSA. 


EITHER  - ( 2 *★ I 3 ) OR  0.0 


X SAVE  SIGN 


LUA  SRI  % GET  THE  ORIGINAL. 

SHAL=35 1 X SIGN  BIT  WENT  IN  EXPONENT  SOON  TO  BE 

OVPRWRT  TTPN 

LIT(O)  = 04000000000000000*  I 6| 

LEX  SCOl 

NORM  I X ABSOLUTE  VALUE  IS  IN  FLOATING  FORMAT. 

ADRN  SSl  X ADD  IN  THE  SIGN. 

STA  OUTl  X GOT  IT. 

*G0  TO  500 

*IF(ARRAY.NF.3)G0  To  200 

NORSAR  DATA.  4 BITS  OF  GAIN  CODE  AND  I2-RIT  TWOS  COMPLEMENT 
MANTISSA. 

LDL(O)  SD49I 
LDA  I M ( 0 ) I 

LDS  SA  I ,6  SAVE  IT  IN  $S. 

« FIRST  LETS  DO  THE  SIGN. 

AND=0800il6l 

SHAL=36l  % PUT  THE  SIGN  BIT  IN  THE  HIGH  ORDER  BIT 

% OF  THE  MANTISSA. 

LIT(O)  =OCOOCOOOOOOOOOOOOt 161 

LEX  SCOl  % NOW  ITS  EITHER  -(2**11)  OR  0.0  . 

NORM  I 
SANl 

LDR  SA I x SAVE  IT  IN  SRl 

LDA  SSl  « RESTORE  THE  ORIGINAL  TO  DO  THF  REST 


AND=7FFi I 6 1 % LOW 

SHAL  =37|  % LEF 

LIT(O)  =0400B 00000000000 1 1 61 
LEX  SCOl 

NORM!  X THA 


X SAVE  IT  IN  SRl 

X RESTORE  THE  ORIGINAL  To  DO  THF  REST 
X OF  THE  MANTISSA. 

% LOW  ORDER  II  BITS. 

X LEFT  JUSTIFY  THEM  IN  THE  MANTISSA. 


ADRN  SRl 
LDR  SAl 
LDA  $S| 
SHAR  =12 1 
SANl 

ADM  =400C» 


% THATS  THE  MANTISSA  WITHOUT  THE  SIGN. 

X THATS  THE  MANTISSA  WITH  THE  SIGN. 

% SAVE  IT  IN  SR. 

% RESTORE  THE  ORIGINAL  TO  DO  GAIN  CODE. 

X ISOLATE  THE  4 PITS  OF  GAIN  CODE. 

% MAKE  IT  NEGATIVE. 

X WE  WANT  EXPONENT  TO  BE  THAT  OF 
X 2**(  I O-GA IN ) . THATS  EXPONENT  ()T  2**!0. 
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A 

A 

A 

A 

A 


200 

A 

C * 
500 


A 

A 

A 

510 


SHAL  =48| 

SAB  =161 
MLRN  SR | 
LDL(O)  SD50I 
STA  OUT( 0) I 
★GO  TO  500 


X THATS  GOT  THE  GAIN  CODE. 

% MULTIPLY  IN  THE  MANTISSA  AND  SIGN. 


★CONTINUE 

DISPLH  "ARG200»f2| 

ARRAY  WAS  OUT  OF  RANGE. 
★CONTINUE 

*IF(DEBUG.LT. I ) GO  TO  510 
DISPLH  "ECI6T64"f0» 

★IF (DEBUG.LT. 3) GO  TO  510 
LDL(O)  $D50« 

STA  OUT(O) I 
★CONTINUE 
★RETURN 
★END 
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★SUBROUTINE  C32T64 ( IN ,OUT) 

★PF  REAL  I N ( * ) , OUT(*) 

★PE  INTEGER  NBUFF I ( ★ , 64 ) , FI NSCN(* ) ,COMP( ★ ) , TOTSCN( ★) , PI NTI (★) 

I PINT2(*)  ,TIMF(*) ,OTIME(*)  , TWTI  ‘A F ( ★ ) . PEN  (★ ) , CNTBL( *, 6 ) 

★PE  REAL  GLCHFTI*)  ,VARFT(  ★)  ,BUFF2(*,550,2)  ,8UFF3(  *,640)  .PREALI  (★)  . - 
I PREAL2(*) ,ALLMSO(*) ,TVARFT(*) 

★ PE  INTEGER  LOFREQ,  HI  FREO , I BUFFI  ( 4096 ),  IBUFF3  (★,  640)  , ARIFF2  ( /0400)  - 

I CHGOOD( 80) ,SITEGD(80) ,S ITES (80) 

★PF  REAL  CHMSQ(80)  ,RBUFFI  (4096)  ,ROWSUM  ,RRt)FF2  ( 704  00) 

★CU  INTEGER  ADRRUF( 8) ,COREPT,  BYTE,  AD3WRD,  ARRAY , DEBUG, TWS/ , 

1 OVLAP ,NCHAN,NS ITE ,NROWS ,DI FFR ,01  FFW,NEN,0LU,0AP,T5CA0L.- 

2 INDEXI ,INDEX2,INDEX3,INDFX4,TI ,T2,T3,T4 ,T5,T6,CH, 

3 , OFFSET, INBYT,NCDCH,TWSZR,NOD5T,NODP,F,BF3PF,NGT,PAGF,T/ 

★CU  LOGICAL  LADRBU(8)  ,LC0REP,LASTI6,LBYTE,LADPWR,LARHAY,L!)FRUG,  - 

1 LTWSZ, LOVLAP,LNCHAN,LNSITE,LNR:)iVS,LDIFFR,LDlFFW,LNF.-.  - 

2 L()LD,  LGAP  ,LTSCAN  ,LTi  ,LT2  , LT3, LT4  , LT5  , LT6 , LCH ,LOFFSF , - 

3_  LINBYT,LF,LNGDCH,LTWSZR,LNGDST,LN0DR,LNGT,LT7 

★EXTERNAL  GTDATA.CI  6T64,C64T32,R0WSUM,R!JNFFT 
★C0MV0N/MAIN2/NBUFF I ,BUFF2,BUFF3,FINSC.N,C0MP,T0TSC0,PINTI  .PIMT2. 

1 TIME, OTIME,TWTIMF,  PEN, CNTRL,GLCHFT,VARFT,  PRFALI  ,PREAL2,’ 

2 ALLMSO.TVARFT, LOFREQ, HIFREQ ,S ITEOD, SITFS ,CH«SO 
★EQUIVALENCE  (NBUFF I ( I , I ) , RBUFFI ( I ) , I PUFF  I ( I ) ) , (BUFF2 (1,1,1), 

1 ABUFF2( I ) ,RBUFF2( I ) ), (CHGOOD( I ) ,SITEGD( I ) ) , 

2 ( RUFF 3 ( I ,1  ) , IBUFF3(  I , I ) ) 

★ EQUIVALENCE  ( I , ADRRUF( I ) ,LADRBU( I ) ) , ( 9, COREPT.LCOREP ) ( I O.BYTF 

1 LBYTE) , ( I 1 ,ADBWRD,LADBWR)  ,(  | 2 .ARRAY.LARRAY) ,(  I 3,TWS7,  - 

2 LTWSZ ) , ( I 4 , OVLAP,  LOVLAP ) , ( I 5 , NCHAN , I.NCUAN  ) , ( I 6,  NS  ITF,  - 

3 LNS  ITE ) , ( I 7,  N ROWS , LNROWS  ),(I8,DIFFR,LDIFFR),(  I 9,DI FFW  - 

4 LDIFFW)  (20,NEW,LNEH),(2I ,OLD,LoLD) ,( 22, GAP, LGAP) , (23,  - 

5 TSCANS.LTSCAN)  , (24,  INDEX  I ) , (25, 1.NDFX2  ) , (26,  INDFX3)  , (27,- 

6 INDEX4) , ( 28.TI , LTI ) , ( 29,T2,LT2) , ( 30,T3,LT3 ) , ( 31 ,T4 , 

I J-t1^*<.32,T5,LT5),(33,T6,LT6),(34,CH,LCH),(35,0FFSL:T, 

8 LOFFSE) , ( 36, 1 NRYT ,LINBYT) , ( 37 , F, I F) , ( 38, NCDCH ,LNGDCU ) , - 

9 ( 39 , TWSZR  , LTWSZ R ) , ( 40,  NGDST,l.NGDST)  , ( 4 | ,NGDP  ,LNGDR ) 

^ j^,RF3PE)  ,(43^NGT^LNGT),  (44,LASTI6)  ,(45,DEPUG,Lr>FBUG)  ,- 

★DISK  AREA  I NDM2 (20) , 0UTDM2 (40)  ,C()NPRM(  I ) 

*IF(DEBUG.LT, I ) GO  TO  10 
DISPLH  "C32T64" , 01 
★IF (DEBUG.LT. 3) GO  TO  10 
LDL(O)  $D49» 

LDA  I N ( 0 ) 9 
DISPLH  "IN”, 32 l 
★CONTINUE 
LDL(O)  $D49» 

LDA  IN(0)» 

LDR  SA I SAVE  IT. 

RAB  =01  ELIMINATE  THE  SIGN  SO  WE  CAN  DO  THE 


SHAR  =561 
SBM  =40i I 61 


/o  SAVE  IT. 

« ELIMINATE  THE  SIGN  SO  WF  CA! 
% EXPONENT. 

% ISOLATE  THE  EXPONENT. 

% SUBTRACT  THE  3?  BIT  OFFSET. 
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A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


ADM  =4000*161 
SHAL  =48» 

LDS  $ A I 
LDA  SRI 
SHAR  =63* 

SHAL  =63* 

OR  $S| 

LDS  SAt 
LDA  SRI 
SHAL  =40| 

SHAR  =161 
OR  SSI 

LDLCO)  SD50I 
STA  OUT( 0) I 

♦IFCDEBUG.LT. I )G0  TO  110 
DISPLH  "EC32T64" , Ol 
♦IFCDEBUG.LT. 3)G0  TO  110 
LDLCO)  SD50I 
LDA  OUTCO) l 
DISPLH  "OUT" ,321 
♦CONTINUE 
♦RETURN 
♦END 


X ADD  THE  64  BIT  OFFSET. 

X PUT  IT  IN  THE  64  BIT  EXP.  FIELD. 

X «UVF  TT 

X NOW  FOR *THE  SIGN. 

X SIGN  BIT  IS  ISOLATED. 

X NOW  WE  HAVE  EXPONENT  AND  SIGN. 

X SAVE  IT. 

X NOW  FOR  THE  MANTISSA. 

X ISOLATE  THE  MANTISSA. 

X PUT  IT  IN  64  BIT  EXP  FIELD. 

X DONE. 
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★SUBROUTINE  C64T32(  IN) 

★PE  REAL  I N ( * ) 

*PE  INTEGER  NBUFF I ( ★ , 64 ) ,FI NSCN( ★ ) ,COMP( ★ ) ,T()TSCN ( ★ ) , PINTI  ( ★) . 

' nr.  , , P INT2 (* ) , TI  MF( ★)  , OTI ME( ★)  ,TWTI ME( ★ ) , PEN ( ★ ) , CNTRL( *, 6) 

*PE  REAL  GLCHFT(*) , VARFT(  ★) , BUFF2 (★, 550,2 ) , BUFF 3( *,640) ,PREALI (*), 

I PREAL2(*) ,ALLMSQ(*) ,TVARFT(*) 

*PE  INTEGER  LOFREO,H I FREQ, I BUFF  I ( 4096) , I RUFF3 ( *,640) , ABUFF2( 7G400)  . 

I CHG(X)D(  80)  ,SITEGD  ( 80) , SITES(  BO) 

*PE  REAL  CHMSQ(  80) , RBUFFI  ( 4096)  , R0WSUM,RBUFF2  ( 704  00) 

*CU  INTEGER  ADBBUF(B)  ,C()REPT,  BYTE, ADBWRD, ARRAY, DEBUG, TWSZ, 

1 OVLAP,  NCHAN ,NS ITE  , NROWS ,DI FFR,  DI  FFW, NEW, OLD , GAP,TSCANfj  ,■ 

2 INDEX! ,INDEX2,INDEX3,INDEX4,TI  , T2 ,T3,T4 ,T5 ,T6, CH, 

2 , M_T/Nir  OFFSET, INBYT,NCDCH,TWSZR,NGDST,. NODR, F,BF3PE,NGT, PAGE, T/ 
*CU  L(X3ICAL  LADBBU(B)  , LCOREP,  I.ASTI  6,LBYTE,LADBWR,LARRAY,LDERUG, 

1 LTWSZ . LOVLAP,  LNCHAN , LNS I TE , LNROWS  ,1.01  FFR,  LD I FFK , LNFW , ■ 

\ F(iM2itGrAM-P£AN  .LTI.LT2,LT3,LT4,LT5,LTA,LCH,LOFFSE  ■ 

3 LINBYT,LF , LNGDCH , LTWSZR , LNGDST , LNGDP , LNGT , LT7 
★EXTERNAL  GTDATA.CI  6T64,ROWSUM,HUNFFT,C32T64 

★COMM ON /MA I N2 /NBUFF I , BUFF2  , BUFF 3,  FINSC.N,COMP  ,TOTSCN,PINTl  P I NT2 
! TIME, OTIME.TWTIME, PEN, CNTRL.GLCHFT, VARFT, PRFALI .PREAI  2. 

2 ALLMSQ, TVARFT ,L()FREQ,  HIFREQ,SITEGD,SITES, CHMSQ 

★ EQUIVALENCE  (NBUFF 1(1,1) .RBUFFI ( I ) , I BUFF  I (!)),( BUFF2( 1,1,1), 

1 ABUFF2(  1 ) ,RBUFF2  ( I ) ) , (CHGOOD(  ! ) ,SITEGD(  I ) ) , 

2 (BUFF3C 1,1), IBUFF3C I , I ) > 

★EQUIVALENCE  (I ,ADBBUF( I ) ,LADBBU( I )), (9, COREPT, LCOREP ) , (10, BYTE 
' LBYTE),(| I , ADBWRD, LADBWR) , ( 1 2 , A RRAY , LA  RRAY) , ( I 3.TWSZ, 

i f-™5!3  • ( 1 4i°vi-AP,LOVLAP),(  I 5,  NCHAN,  LNCHAN) , ( I 6,  NS  ITE, 

3 LNS ITE) ,(  I7,NR0WS,LNR0WS),(I8,DIFFR,LDIFFR) ,( I9.DIFFW, 

i J-^fFJf3.(20,NEW,LNEW),(2l,OLD,LOLD),(22,GAP,LGAP),(23,  - 

5 TSCaNS.LTSCAN) , (24, INDEX! ) , ( 25 , INDEX2 ) , (26, INDEX3) , (27,- 

6 INDLX4) , ( 28 ,TI , LT 1 ) , ( 29,T2,LT2 ) , ( 30,T3, LT3) , ( 31 ,T4 , 

7 LT4 ) , ( 32, T5,LT5 ) , ( 33,T6,LT6) , ( 34 , CH.LCH) , ( 35, OFFSET, 

o 3,64-ic?YT • L I NBYT  > , ( 37 , F , LF ) , ( 38 , NGDCH , LNGDCH ) , - 

9 (39, TNSZR, LTWSZR), (40, NGDST, LNGDST), (41 , NGDR, LNGDR ) , 

0 ( 42 , BF3PE ) , (4  3, NGT , LNGT) , ( 44, LAST! 6) , (45, DEBUG, LDEPUC) 

1 ( 46, PAGE) , ( 47, T7, LT7)  ’ 

★DISK  AREA  INDM2  (20) , 0UTDM2(  40)  ,C()NPRM(  I ) 

*IF(DEBUG.LT. I )G0  TO  10 
DISPLH  "C64T32" ,0| 

★ IF ( DEBUG.LT. 3) GO  TO  10 
LDL(O)  $D49| 


A SAVE  IT. 

» GET  RID  OF  THE  SIGN  BIT  FOR  NOW. 
- ISOLATE  THE  EXPONENT. 

S SUBTRACT  OUT  THE  64-BIT  OFFSET. 
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A ADM  =40H6l  x ADD  IN  THE  32-BIT  OFFSET. 

A SHAL  =561  % PUT  IT  IN  32-BIT  OUTER  EXPONENT  FIELD. 

A LDR  $A t X SAVE  IT. 

A LDA  $S|  X NOW  FOR  THE  SIGN  BIT. 

A SHAR  =63l 

A SHAL  =63l 

A UR  $R|  % $a  HAS  EXPONENT  AND  SIGN. 

A LDR  $ A I % SAVE  IT. 

A LDA  SSI  x NOW  FOR  THE  MANTISSA. 

A SHAL  =161 

A SHAR  =401  XTHE  32-BIT  MANTISSA  FIELD. 

A OR  SRI  X THATS  ALL  OF  IT. 

A LDL(O)  $D49| 

A STA  IN (0) I 

*IF (DE8UG.LT.  I ) GO  TO  110 
A DISPLH  "EC64T32" #0l 

★IFCDEBUG.LT. 3)G0  TO  110 
A LDL(O)  SD49I 

A LDA  IN ( 0) I 

A DISPLH  "()UT(  IN)".32l 

I I 0 *CONTI NUE 
★RETURN 
★ END 
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★SUBROUTINE  CNVTIM 

;PE  1NTEGE  m^iiiroLDTi’L1™^'*'64’6’  •p,nti  <#)  ■ ,nbuf<*-  '«>  • 

a 0UPTWAI6)  SAVBCT.|JVPT«inUpAOEtrt> .TSTEPS(6) .SCANS. 

icu  INTEGEft 

3 DEBUG , BCT . ADByENDADB ' » I ^DYT , CUB  YT.  ORCCnfl^PAGE  * 

*CU  LOGICAL  UDBBU18) ™.LINPTN,LSAVAD,LADBOU(A)  .LOUPTW,. 

i innoc  DWORDS,  LT I » LT2 , LT  3 , LT4 1 LT5  , LT  , LOUB  YT  LIT  LPRTTa  . 

3 ^^LT?:ti^^li^CN‘A,’UDBWfl-L,'RGc'):Lp^ 

★EXTERNAL  RDPRM , GETBYT, PUTBYT 

♦COMMON /MAIN/CNTRL.nu™  •T,“E.OLDTIU,SAVBCT,SAVPTW,  . 

;EGU,  VALENCE,  J.  ADBBUF «,  I.LADBBU,,  ! > . I P.ARR^APRAY, . 

3 <22V|’"lti  I •HTS,,-^?'?US^BVTS>  • '21 . WORDS, LWOHDS) , - 

4 2A.’JJ:tfJ  ’ f^^^KL,^(?^I3•LI3).(25,T4.LT4). 


3 22  T I T iV  ( jT  t .’ri! ' ; ;LDn 5 ' ’ 1 ^ 1 ' K(,f,DS • LWORDS ) , - 

4 tii ri  ’ ?M2»LT2)»(24*T3.LT3)f(25,T4tLT4). 

I j|o:^Xe:^S!l?:Sf:A^?S2SSsk?V?S.**£IS*ifefcgT*  • - 

t ; ( -BYTCNT‘  ■ > • ; r:BD  *LW[,HD » • : 

9 \ Tl  • ^»E,’  LP*CE  >•  «5 , DEBUG,  LDEBUO.C A^BCTaBCTK  "R0C" * ’I 

i AREA  nupu^jfgj : Vnpu^o?' ’C,UPUT3(20) .OUPUT5 C20) , - 
*IF(DEBUG.LT. I )GO  TO  10 
UISPLH  "CNVTIMl", Ol 
> T6*CNTRL(2, ARRAY) 

★IF(Trt.NE . I ) GO  TO  I 00 

CNTRL( 2, ARRAY)  EQUALS  I ONLY  IF  TIME  IS  ALRFADY  TM  the  cnnu  „c 
DECI SECONDS  FROH  THE  BEGINNING  OF  THF  YEAH  AS^READ  ?Rnu  TOP  TAPP 
NO  CONVERSION  IS  NECCESSARY.  R ™R  TAPR* 

WE  NEED  ONLY  TO  MOVE  IT. 

★CALL  GETBYT 
LDL(O)  LINBYTl 
CSHL(O)  I6» 

LDA  $C0| 

STA  TIME | 

★CALL  GETBYT 
LDL(O)  LINBYTl 
LDA  $C0» 

OR  TIME? 

STA  TIMEl 

★IF (DEBUG.LT. I )G0  TO  20 
DISPLH  "FMT#!",0» 

★IF (DEBUG.LT. 2) GO  TO  20 


aT*S^ 

J 


I 
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A DISPLH  "TIMEJ',2| 

20  *G0  TO  500 
100  T6«*CNTRL  ( 2 , A RRAY ) 

★IF(T6.NE.2)G0  TO  200 

C TIME  IS  IN  EBCDIC  IN  THE  FORM  DDDHHMMSS.  GET  IT  A CHARACTER  AT  A 

C TIME. 

★CALL  GETBYT 
LT6»0FF .TURN  ON.. LAST. 4 
LTI-LINBYT.SHR.I2 
LT2-LINBYT.SHR.8 
LT2-LT2.AND.LT6 
LT3=LI NBYT.SHR.4 
LT3=LT3.AND.LT6 
LT4=LI NBYT.AND.LT6 

A *DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 

TIME(*)«(T1*1 00+T2* I 0+T3 ) *2  4+T4* 1 0 
★CALL  GETBYT 
LT6« OFF. TURN  ON.. LAST. 4 
LTl-LINBYT.SHR.12 
LT2-LI NBYT.SHR.8 
LT2-LT2.AND.LT6 
LT3-LI NBYT.SHR.4 
LT3-LT3.AND.LT6 
LT4-LI NBYT.AND.LT6 

A ADUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
TIME(*)«((TIME(*)+TI )*60+T24l0+T3)*60+T4*l  0 
★CALL  GETBYT 
LT1-LINBYT.SHR.12 

A %DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
TIME(*)«(TIME(*>+T1  >*10 
★IF  (DEBUG.LT. 1 ) GO  TO  110 

A DISPLH  "FMT#2",0l 

★IF(DEBUG.LT.2)G0  TO  110 

A DISPLH  "TIME*", I6,TIME,TIME» 

110  *G0  TO  500 
200  ^CONTINUE 

C IF  WE  GET  HERE  SOMETHING  WIERD  IS  GOING  ON,  SINCE  IT  AINT  A VALID 

C FORMAT.  PRINT  OUT  A LITTLE  MESSAGE  AND  SEE  WHAT  HAPPFNS. 

A DISPLH  "ARGTIME",2» 

500  ^CONTINUE 

★ I F ( DEBUG . LT . 1 ) GO  TO  510 

A DISPLH  " ECNVTIMt*' ,01 

510  ★ RETURN 
★END 


CTO< 
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★SUBROUTINE  FNGRIL 

PURPOSE* ONCE  A POWER  MAXIMUM  HAS  RFEN  FOUND  FOR  EACH  FREQUENCY 
ON  THE  COARSE  GRID  REFINEMENTS  OF  THE  MAX  ARE  PERFORMED  ON  A 
SQUARE  GRID.  THE  SPACING  ON  THE  SQUARE  GRID  IS  DKX/6. 

THE  SEARCH  STARTS  IN  THFEAST  DIRECTION  AND  PROCEEDS 
CLOCKWISE.  THE  DIRECTION  OF  THE  SEARCH  IS  CHANGED  WHEN  A 
MAX  IS  FOUND  IN  THA  DIRECTION. 

DECLARATIONS* 

★ PE  INTEGER  IN3UF( *,640 ) ,CNTRL( ★, 6) , NCHAN( ★) ,PINTI ( ★) , OFFSET( ★)  , 

1 LOCATE! ★) ,NPTS ( ★) ,C0UNT2 ( ★) ,C0UNT3( ★) , LOC2D(*,25 ) 

2 L()C3D(*,25)  ,TWTIME(  ★)  , ADJF(*) 

★PE  REAL  POWER! *,25) , FMAX (*,25),FKX(*,25),FKY(*,25),RINBUF(*,640),  - 

1 X( *,25 ) , Y( *,25 ) , FFT! ★, 61 2) , KERNEL (★,25 ) , 

2 XCfX)RD( ★) , YCOORD!  ★) , PREALI  ( ★) , COSK ( ★) ,SINK( ★)  , 

4 COSDKt ★) ,SINDK (★) , BEAMER ( ★ ) ,FPMAX( ★) , KXMAX (★) , 

5 KYMAX!*) ,DELX(*) , DELY(*) , KXSEP(*) ,KYSEP(*) ,KSEP(*) , 

6 VEL!*) , AZ(*) , SIGNAL! ★) ,FSTAT(*) ,SUMSQ(*)  , 

7 TEST(*) ,K(  ★) ,CHANAV(  ★) ,TPOWER(*) .FREQ!*) ,TDKX(*) 

★PE  REAL  ADKX ( 4 ) . ADKY ( 4 ) , YPOI NT! 50) , YMAX ( 50) ,DX(500) , DY (500) 

★PE  REAL  BEAM  (*)  ,TPOW(*)  .DELTAK!  ★)  ,NRPOWR( *,25  ) ,NI  POWR!  *,25  ) ,P 
★PE  REAL  PREAL2!*)  .RPOWER(i<,25  ) , I POWER!  *,25 ) .RTPOW!  ★)  , ITPOW!  ★) 

★PE  INTEGER  MAX 

★CU  INTEGER  LOFREO, HIFREO, DEBUG, SM, T1  ,T2, ARRAY, PAGE, I .N.MNCHAN, 

1 MNPTS.NPOINT, SWITCH, NFREQ.IGC), LINE, LINES,  INDEX,  IP, 

2 TWIN, SAM, IFREQ, J.NFREQ1 , REFINE, I ND.YTOP 

3 YPMI,SIGN,NTIMES,LINEPI 
★PE  REAL  DELTX! 3000), DELTY!  3000), DIST 
★PE  REAL  DELTAX, DELTA Y,KX,KY 

★CU  REAL  DKX, LOWER, UPPER, LINEP.HDKX, BORDER, TWOH 
★CU  REAL  DELTAF, RADIUS, ANGLE 
★CU  LOGICAL  M0DE3.NM0DE 
★EXTERNAL  MAX, GRID, REALE, IMG 

★COMMON/MAINFK/INBUF,CNTRL,NCHAN, PI  NT I .OFFSET, LOCATE, NPTS, COUNT?,  - 

1 COUNT3,LOC2D,L()C3D,  POWER, FMAX, FKX.FKY,  X,  Y,  KERNEL, 

2 XCtXJRD,  YCOORD,  PREALI  , COSK, SINK, BEAM, 

3 TPOW, DELTA K, RPOWER, I POWER,  COSDK.SI NDK,  - 

4 BEAMER, FPM AX,  KX.MAX, KYMAX, DELX.DEI.Y,  KXSEP , KY5EP , KSEF' , - 

5 TWTIME.TPOWER, VEL, AZ, SIGNAL, FSTAT.SUMSQ, TEST, K, 

6 CHANAV, FREQ,  ADJF.DX  ,DY,P,YP()INT,  YMAX, ADKX,  ADKY, 

7 KX  KY 

★EQUIVALENCE! INBUF! I , I ) , R INBUF (1.1 ) ) » ( I NRUF ( I ,2 B), FFT! I, I )) 

★ EQUIVALENCE  ( I ,L()FREQ) , ( 2, HI  FREQ) , (3, DEBUG) , (4,SM  ) , (5  ,TI  ) ,(6,T2)  , - 

1 (7, ARRAY) , (8, PAGE) , (9, I ) ,( 10, N) , (12, MNCHAN ) , - 

2 ( I 3.MNPTS) , ( I 4.NP0INT) , ( 1 5, SWITCH)  ,(16,100), 

3 ( 17, INDEX), (IB, IP), (19, DKX), 

4 (20, LOWER) , (21  .UPPER) ,(22, LINE), (23, LINES) , 

5 (24, HDKX) ,(25, BORDER) ,(26, TWOH) , 

6 (27, DELTAF), (2B, RADIUS),  (29, SIGN), 

7 (30, M0DE3), (31 ,NMODE), (32, TWIN), (33, ANGLE), (34, SAM),  - 

★* 

(35.NFREQ) , ( 36, IFREQ ) , ( 37, J ) , ( 38 .NFREQI ) , ( 39, REF INE)  - 
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9 , (40,NTIMES)  , ( 41  , IND) , ( 42,  YTOP) , (43,YPMI  ) , ( 44, UNEP)  , 

0 (45.LINEPI) 

★DISK  AREA  CONPRMC  I ) ,STC»RD(  I ) ,FKIN(8I  ) 

TDKX(*)=DKX 

★DO  500  NTIMES=I , REFINE 
TDKX ( ★)=TDKX ( *)/6.0 
PREALI  ( ★)=TDKX(*) 

ADKX ( I ) = PREALI ( I ) 

ADKX(2)=0.0 
ADKX (4 )=0.0 
ADKY( I )=0.0 
ADKY ( 3)=0.0 
ADKY ( 4 )=PREAL I ( I ) 

PREAU  ( ★ )=-PREALI  ( ★) 

ADKX(3)=PREALI (I ) 

ADKY (2 )=PREALI ( I ) 

★DO  475  1=1 ,4 
MODE=ON 

DELX (*)=ADKX ( I ) 

DELY(*)=ADKY ( I ) 

RTP()W(  *)=0.0 
I TP()W ( *)  = 0.0 
MB  I T I = MODE 
★DO  250  N= I , MNCHAN 
MODE=(N.LE.NCHAN(*>> 

DELTAK( ★)  =+6.283l 85 307* ( DELX (*)*X( ★, N) +DELY( *)*X(  *,N) ) 

COSDK( ★ ) = C05 ( DELTAK(  ★) ) 

SINDK(*)=SIN(DELTAK(*)) 

NRPOWRI ★ ,N)  = RPOWER( *,N)*COSDK( ★) -I  POWER ( ★»'I)*SINDK( ★ ) 

NI POWH  ( ★, N ) = RPOWER( *,N )*SINDK(  ★ ) + 1 POWER!  ★,'lI)  *COSDK(  ★ ) 

RTPOW(  *)=RTPOW( *)+NRPOWR(  *,N) 

ITPOW!  ★)  = ITPOW(*)+NIP()WR(  *,M) 

250*C0NTINUE 

MODE=MRITI 

TPOW!  *)=RTPOW(  ★)  ★★2+ITPOW!  ★)  irk? 

400*C()NTI  NUE 

★IF (DEBUG.LT.2)G0  TO  405 
A DISPLH  "FINE", 01 
PINTI ( I )=I 

A DISPLH  "I", I 6, PINTI , PINTI I 

A DISPLF  "DELX", 16, DELX, DELX* 

A DISPLF  "DELY"  , I 6,DELY,DFLY| 

A DISPLF  "RTPOW" , I 6, RTPOW , RTPOW+3* 

A DISPLF  " ITPOW" ,16, ITPOW, ITPOW+3* 

A DISPLF  "TPOW", I6,TPOW,TPOW+3| 

A DISPLF  "FPMAX"  , I 6,FPMAX,FPMAX+3l 

405*C0NTINUE 

M()DE=MODE . AND.  ( TPOW! ★)  .GT.  FPMAX  ( ★ ) ) 

★ I F(  . NOTANY . ( MODF ) )G()  TO  475 
MB IT I* MODE 
★DO  410  J=l , MNCHAN 
RPOWEH  ( ★,  J )=NRPOWR( ★, J ) 

I power ( ★, j )=ni powr( ★, j) 


CT'  ; < 
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4 1 0*CONTINUE 
M()DE*MBITI 

<XMAX(*)=KXWAX(*)+DELX(*) 
r£YMAX(  *)  = KYMAX  ( * ) +DELY(*) 

FPMAX(  ♦ )*TPOW( ★) 

TPOW  (*)=0. 

RTPOWt  *)=0.0 
ITPOW<*)=0.0 
*D()  450  J=l  .MNCHAN 
M()DE=M()DE .AMD . ( J .LE.NCHAN(* ) ) 

DELTAK ( *)=+6.283 I 853070*( DELX ( *) ♦X ( *, J ) +DELY ( *)*Y( *, J ) ) 
COSDK ( *)=COS ( DELTA K(  * ) ) 

SI NDK ( *)=SI N( DELTAK ( * ) ) 

NRPOWR  ( * , J )=RPOWER(* , J )*COSDK<  *)  - I POWER ( *,  J)*SINDK(*) 

N I POWR ( * , J ) = RPOWER ( * , J ) *S I NDK ( ★) + 1 POWER ( * , J ) *C05DK ( * ) 
RTPOWC  *)  = RTP()W ( ♦ ) +NRPOWR( J ) 

I TPOW ( * ) = I TPOW ( ♦ ) +N I POWR  ( * , J ) 

450*C()NTINUE 

M()DE=.MBITI 

TPOW  ( * ) = RTPOW ( * ) **2+ ITPOW( ♦)  Jr* 2 
♦Go  TO  400 
475  *C0NTINIJE 
MODE* ON 
500*C()MTIMUE 

*IF( DEBUG.LT. 2) GO  TO  510 
A DISPLF  " EFNGRID" ,0l 
510  ★CONTINUE 
♦RETURN 
♦F' 
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★SUBROUTINE  GETBYT 

GETBYT  IS  SUPPLIED  A POINTER  TO  THE  BYTE  WANTED  VIA  INPTB.  IT 
FIGURES  OUT  IF  THE  BYTE  IS  IN  ADB  OR  CORE  OR  ON  DISK.  IT  GETS  THE 
DEMANDED  BYTE  AND  LEAVES  IT  RIGHT  JUSTIFIED  IN  » I NBYT".  ALL 
crr>NTnR  fTART  AT  ZER0  RATHER  THAN  ONE  BECAUSE  OF  THE  SHIFTING, 

fc  1 v/  • U ON  fc  • 

★PE  INTEGER  CNTRL(*,6) , OIJTBUF  (★  ,64, 6)  ,PINTI  (★) , INBUF(  ★,  128), 

I T I ME < ★ ) ,OLDTIM(*) , ’ * 

' SAVBCT,SAVPTW,OUPAGE( 6) ,TSTEPS( 6) , SCANS, 

2,  OUPTWA( 6)  , OTIMEAt  6)  .ORGADB,  INBUFI  (8192) 

★CU  INTEGER  ADBBIJFf  8)  , ARRAY , I NPTB , INPTW,SAVADB, ADB0UT(6) , OUPTW. 

1 BYTS, WORDS, Tl ,T2,T3,T4,T5 ,T6,  IT,PRTIAL, ADDRS 

2 WORD,  BYTCNT ( 6 ) , ADBWRD , I NBYT , OIIBYT , ORGCOR , PAGF , 

3 DEBUG, BCT, ADB, ENDADB 

★CU  LOGICAL  LADBBU(8) ,LARRA Y , LI NPTB ,L I NPTW, LSAVAD, L ADBOU( 6 ) , LOUPTW , > 

1 LBYTS , LWORDS , LT I ,LT2,LT3,LT4,LT5,LT6,LOUBYT,LIT,LPRTIa!« 

2 L ADDRS, LWORD , LI NBYT ,LBYTCN (6) ,LADBWR,  LORGCO,LPAGE 

3 ,LDEBUG,LBCT,LADB,LENDAD 
★EXTERNAL  RDPRM , PUTBYT , CNVT I M 

★COMM  ON /M A I N/CNTRL , OUTBUF , I NBUF , P I NTI , T I ME , OLDTI M , S A VBCT, SA VPTW . - 

I TSTEPS,  SCANS,  OUPTWA,OIJPAGE,OTIMEA,  ORGADB 
★EQUI  VALENCE(  I ,ADBBIJF(  I ),LADBBU(I  ) > , ( 9,  ARRAY, LARRAY) , 

' C 10, INPTB, LINPTB), 

( II ,INPTW,LINPTW),( 12, SAVADB, LSAVAD), 

( I 3,ADBOUT(6) ,LADBOU( 6> ) 

:('94.(.)UPIW  «L()UPTW) , (20, BYTS, LBYTS), (21  .WORDS, LWORDS) , - 
, } 23'T2»LT2)  ♦ ( 24,T3,LT3) , <25  ,T4,LT4) , 

( 26,T5,LT5 ) , ( 27,T6,LT6) , ( 28,0IJBYT,L0UBYT) , (29,  IT, LIT) , - 
\ o? * J-PRTI A ) , ( 3 1 , ADDRS , L ADDRS ) , ( 32 , WORD , LWORD) , - 
( 33, I NBYT ,LINBYT) , ( 34, BYTCNT ( I ),LBYTCN(I )), 

(40, ADBWRD, LADBWR) 

,,,  # ,(  43, ORGCOR, LOROCO),- 

( 44 , PAGE, LPAGE ) , ( 45 , DEBUG ,LDEBUG) , ( 46 , BCT.LBCT) . 

( 47 , ADB , L ADB ),( 48 , ENDADB , LENDAD ) 

★EQUIVALENCE  (INBUF( I ,1 ), INBUFI (I )) 

★DISK  AREA  OUPUTI (20) ,0UPUT2( 20 ) ,0UPUT3( 20) ,0UPUT4( 20) ,0UPUT5 (20) • - 
I OUPUT6(20),INPUT(50) 

★IF( DEBUG.LT. I )G0  TO  5 
DISPLH  J,GETBYT"  ,0* 

★I F( DEBUG.LT. 3) GO  TO  5 
DISPLH  "GETBYT I", 2 I 
★CONTINUE 

FIRST  MAKE  WORD  P()INTERCBYTE  POINTER/4. 

LINPTW-LINPTB.SHR.2 

* DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 

★IF (I NPTW .LT. ENDADB )G0  TO  30 
IF  WE  GET  HERE  ADB  MUST  BE  REFILLED. 

★ I F ( I NPTW.LT.ORGCOR+8 I 92) GO  TO  20 
IF  WE  GET  HERE  CORE  MUST  BE  REFILLED. 

PAGE  WANTED  IS  BYTE  POINTER/4096  ♦ I . 

LPAGE-LINPTB.SHR.I2 


I 

1 

2 

3 

4 

5 

6 
6 

7 

8 
9 
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* n°9^Mn.^K,STATEMENT  T0  F()RCF  DEALLOCATION  OF  REGISTERS. 

P AGE=PAGE+I 

*READ(64,INBUF( I , I ) , I NPUT( PAGE) ,8) 

LORGCO-LINPTB.SHR.I2 
L()RGC()*LORGCO.SHL.  I 0 

/.DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS 

f^fCOR  NOW  HAS  THE  WORD  ADDRESS  OF  THE  FIRST  BYTE  IN  INBUF. 
*rtAIT  64 

A BIT  OF  DEBUG  OUTPUT  NOW 
♦IF (DEBUG. LT. I )GO  TO  20 
DISPLH  "REFILLCi",Ol 
♦IF (DEBUG.LT. 3) GO  TO  20 
DISPLH  " I NBUF * M , I 8 , 1 NBUF , I NBUF+64M- 1 1 
20  ♦CONTINUE 

AT  THIS  POINT  WE  REFILL  ADRBUF. 

LT6-LI NPTB .SHR .5 
LT6= LT6.SHL.3 

/.DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 

0RGADB-T6 

ENDADB-T6+8 

ORGADB  NOW  HAS  THE  WORD  ADDRESS  OF  THE  FIRST  BYTE  IN  ADB  ENDADR 
HAStTHErWORD+ADDReSS  of  TOE  RYTE  AFTER  THE  LAST  BYTE  IN  AM? 

WITHIW  INBUF  WE  WANT  TO  REFILL  ADRBUF  FROM. 
*TRANSFER(8)  ADBBUF ( I )* INBUF I (T5) 

♦IF (DEBUG.LT. I ) GO  TO  30 
DISPLH  "REFI LLA*" ,2| 

30  ♦CONTINUE 

"HEN  WEGET  HERE  THE  BYTE  WANTED  IS  IN  ADBBUF.  JUST  GOTTA  FIGURE 

( JU  i nn  tnfc  • 

T6= ORGADB 

ADB W RD« I NPTW-T 6+ 1 

LT6-OFFMTURNK()NTALASTN2  T°  F°RCE  DFALL()CATION  OF  REGISTERS. 

LT6=  L I NPTB .AND. LT6 
LT6-LT6.SHL.4 

T6  IS  USED  TO  DETERMINE  THE  AMOUNT  TO  ROTATE  THE  WORD  IN  ADR 
-J”  JUSTIFY  BYTE  WE  WANT.  NUMBERING  THE  BYTES  0,1. 2, 3 

THE  BYTE  WANTED  IS  ( I NPTB  MOD  4).  TO  RIGHT  JUSTIFY  THA~  BYTE  WE 
ROTATE  LEFT  BY  ( C I NPTB  MOD  4 >-#- 1 )♦!  6 OR  ( I NPTB  MOD  4)*I6+I6.T6 

HAS  JUSTED  BEEN  ASSIGNED  ( I NPTB  MOD  4 ) *1 6 SO  WE  HAVE  ONLY  To  ADD 

16. 

LT5=  OFF  < i yji,  ( ON..  LAST.  I 6 

-A  DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 
LINBYT-LADBPU(ADBWRD) .RTL.T6+1 6 
LI NBYT=LINBYT. AND.LT5 
*IF(DFBUG.LT. I )G0  TO  40 
DISPLH  "EGETBYTt" ,0» 

♦IF (DEBUG.LT. 3)G0  TO  40 
DISPLH  "GETBYTi " ,2» 

40  ♦CONTINUE 

I NPTB= INPTR+ 1 
♦RETURN 
♦END 
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★SUBROUTINE  GRID 
DECLARATIONS! 

★PE  INTEGER  I NBUF! *,640) ,CNTRL( ★, 6) , NCHAN! ★) , PINT! (★), OFFSET! ★) , 
LOCATE!*) ,NPTS(*) ,COUNT2<*) .COUNT3! ★) ,LOC2D(*,25 ) , 
LOC3D! *,25 ) .TWTIME! ★) ,ADJF ! ★) 

POWER! *,25) ,FMAX(*,25) ,FKX( ★ ,25 ) ,FKY ( *,25 ) , RINBUF !★, 640) 
X(*,25),Y(*,25),FFT(*,6I2) , KERNEL! ★, 25) , 

XC  CXJRD  < ★ ) , YC(X)RD(*) , PREALI  (★),  C()SK<  ★)  .SINK!  ★)  , 

COSDK ( ★) ,SINDK(*) ,BEAMER(*),FPMAX<*) ,KXMAX(  ★) , 

KYMAX!*) ,DELX( ★) ,DELY< ★) ,KXSEP(*) .KYSEP!*) ,KSEP(*) , 
VEL!*) ,AZ<*) .SIGNAL! ★) ,FSTAT( ★) , SUMSQ! ★) 

TEST (★) ,K( ★) ,CHANAV< ★) .TPOWER!*) .FREQ!*) 

ADKX! 4) . ADKY ! 4 ) . YPOI NT (50) ,YMAX(50) ,DX(500) ,DY(500) 
BEAM(*).TP0W(*).DELTAK(*),P 

PREAL2!*) .RPOWER! *,25) . I POWER! *,25) .RTPOW!*) . ITPOW!*) 

MAX 

LOFREQ. HI  FREQ. DEBUG, SM.TI , T2, ARRAY, PAGE, I .N.MMCHAN, 
MNPTS, NPOINT, SWITCH, NFREQ.IGO, LINE, LINES, INDEX, IP, 
TWIN, SAM, IFREQ.J.NFREQI .REFINE, IND.YTOP 
YPMI.NTIMES.LINEPI 
DELTX ( 3000 ) , DELTY ( 3 000 ) , D I ST 
DELTA X, DELTA Y.KX.KY 

DKX, LOWER, UPPER, LINEP.HDKX, BORDER, TW OH, SIGN 
DELTAF .RADIUS, ANGLE 


1 

2 

★PE  REAL 

1 

2 

4 

5 

6 
7 

★PE 
★ PE 
★PE 
★PE 
★CU 

1 

2 
3 

★PE 
★PE 
★CU 


REAL 
REAL 
REAL 
INTEGER 
INTEGER 


REAL 

REAL 

REAL 

REAL 


★CU 

★CU  LOGICAL  M0DE3.NM0DE 
★EXTERNAL  MAX, FNGRID, REALE, IMG 

★COMM ON/M A I NFK/ 1 NBUF , CNTRL , NCHAN , P I NT I , OFFS ET, LOCATE , NPTS , Cf  )UNT2 , 

1 C0UNT3, LOC2D.LOC3D, POWER, FMAX.FKX.FKY.X.Y, KERNEL, 

2 XCOORD.YCOORD, PREALI  , COSK, SINK, BEAM, 

3 TPOW.DELTAK, RPOWER, IPOWER,  COSDK, SINDK, 

4 BEAMER.FPMAX.KXMAX, KYMAX, DELX.DELY.KXSEP.KYSEP.KSFP, 

5 TWTI ME .TPOWER , VEL , AZ , S I GN AL , FSTAT, SUMSQ .TEST, K , 

6 CHANAV, FREQ, ADJF,DX,DY,P,YP0INT,YMAX, ADKX, ADKY, 

7 KX.KY, DELTX, DELTY 

★EQUIVALENCE! I NBUF! I , I ), RINBUF! t , I )), ( I NBUF! I ,28) , FFT! 1,1)) 
★EQUIVALENCE  ! I , LOFREQ ) , ! 2, HI FREQ) , (3, DEBUG) ,(4,SM),!5,TI ),!6,T2) , 

I (7, ARRAY), !8, PAGE), <9, I), <10, N),  <I2,MNCHAN), 

! I 3, MNPTS) ,!  1 4, NPOINT),! 1 5, SWITCH) ,! 16, IGO) , 

! 1 7 . INDEX) . ! 18, IP), (19, DKX), 

(20, LOWER), (2  I .UPPER), <22, LINE), (23, LINES), 

( 24 , HDKX ) , ( 25 , BORDER ) , ( 26 , TWOH) , 

(27, DELTAF), (28, RADIUS),  (29, SIGN), 

( 30.M0DE3) , ! 3 1 ,NMODE) , (32, TWIN) , < 33, ANGLE) , ( 34, S AM ) , 


8 (35.NFREQ) , ( 36, 1 FREQ ) , ( 37, J) , ( 38, NFREQI  ) ,(39, REFINE) 

9 , (40.NTIMES) .(41 , IND) , (42.YTOP) , (43.YPMI ) t ( 44.LINEP) , 

0 (45.LINEPI) 

★DISK  AREA  CONPRM! I ) .STCORD! I ) ,FKIN(8I ) 

★IF(SWITCH.EQ.2)G0  TO  500 
PREALI (*)=DKX*(SQRT( 2.0) /2 .0 ) 

BORDER=PREALI( I) 
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PREALI  (*)=FL()AT(  SAM  )/FL()AT(  TWIN) 

QELTAF=PREALI ( I ) 

PREALI (*)=ANGLE 

*IF( .ANY. ((ANGLE. GT. 0.0) ) )G0  TO  20 
PREALI  (*)  = (|  .O/LOWER )*FL()AT(  IFREQ)*DELTAF+BORDFR 
RADI U5=PREAL 1(1) 

KY=0 .0 
KX=0 .0 
★GO  TO  30 

20  RAIM  US-PREAL  f?^°WER-  * *0/UPPER  * *°* 5*FL()AT(  I FREQ)  ★DELTAF+RORDER 

PREALI ( ★ )=RADIUS*SIN( ANGLE ) 

.<X=PREAL  I ( I ) 

PREALI (★)-RADIUS*COS( ANGLE) 

<Y=PREAL I ( I ) 

30  PREALI (★)*2.0*SQRT(2.0)*DKX*.75 
DELTAY=PREALI ( I ) 

PREALI ( ★ )= (DKX*SQRT( 6.0)/2.0>/2.0 
DELTAX=PREALI (I ) 

C CHANGE  IN  X IS  HALF  THE  BASE  TRI ANGLE, CHANGE  IN  Y IS  TWICE 
C THE  HEIGHT  OF  BASE  TRIANGLE. 

PINTI ( ★)*IFIX( RADIUS/DFLTAX)-I 
LINE=PINTI (I ) 

PREALI  ( *)=FL()AT(  LINE)*DELTAX 
P-PREALI ( I ) 

PINTI (*)=2*LINE+I 
LINES«PINTI (I ) 

I G()=  I 

PINTI (*)=(L!NE/2)*2 
T I =P INTI ( I ) 

★ IF(LINE.EQ.TI  ) IG(*0 
★DO  100  INDEX=ltLINE 

PREALI  (*)=SQRT(  RADIUS**2-(  P-FLOAK  I NOEX- 1 ) *DFLTAX ) **2  ) 
DIST=PREALI ( I ) 

★ IF ( IGO.EQ  .0)G()  TO  50 

PINTI ( ★)=IFIX( (DIST-DELTAY/2 ,0)/DEITAY) 

YTOP=PINTI (I ) 

PREALI (★)=2.0*FL0AT(YT0P)+2.0 
YPOI NT ( I NDEX )=PREAL I ( I ) 

PREALI  ( ★ ) =DELTAY /2 . 0+DELTA Y*FL()AT ( YTOP ) 

YMAX( I NDEX)=PREALI  ( I ) 

IG()=0 
★GO  TO  70 
50  *C()NTI  NUE 

PINTI ( ★)=IFI X(DIST/DELTAY) 

YT()P=P  I NTI  ( I ) 

PREALI  ( ★)  =2 ,0*FL()AT( YTOP ) + 1 .0 
YPOINT(INDEX)=PREALl (I ) 

PREALI  (★)=FLOAT( YTOP)*DFLTAY 
YMAX( I NDEX)=PREALI ( I ) 

IG()=  I 
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70  IND=LI NES-INDEX+I 

YPOINT(IND)=YPOINT< INDEX) 

YMAX( I ND)=YMAX( INDEX) 

100  ^CONTINUE 

PREALI (*)=KX-P 
DELTX( I )=PREALI ( I ) 

PREALI (*)=KY-YMAX( I ) 

DELTYC I )=PREALI  ( I ) 

DXCI )=DELTX< I) 

DY( I )=DELTY( I ) 

LINEPI=LINE+I 

PI  NT  I (★)  = IF  I X( RADI  US/DELTA Y) 

YT()P=PINTI  C I ) 

PREALI (*)=FL0AT(2*YT0P+I ) 

YP0INT(LINEP1 )=PREALI  ( I ) 

PREALI  C *)=FLOAT( YTOP )*DELTAY 
YMAX(LINEP) )*PREAL1 ( I ) 

1 = 1 

SIGN=I . 

INDEX=0 

150  INDt X=INDEX+I 

PINTl ( *)=IFI X( YPOINTI  INDEX ))-l 
YPMI=PINTl (I ) 

★DO  160  IND=I,YPMI 
I-I  + l 

DELTX( I )=0.0 
PREALI (*)=SIGN*DELTAY 
DELTY( I ) -PREALI ( I ) 

PREALI  (*)=DX<  I-l  )+DELTX<  I ) 

DX i I )=PREALI ( I ) 

PREALI (*)-DY(I-l )+DELTY< I) 

DY  ( I )*PREALI ( I ) 

160  ^CONTINUE 

★IF( INDEX. EO.LINES)GO  TO  200 
1 = 1 + 1 

DELTX( I )=DELTAX 

PREALI (*)=(YMAX( INDEX+I )-YMAX( INDEX))+SIGN 
DELTY( I )=PREALI ( I ) 

PREALI (★)*DX( I-l >+DELTX(  I) 

DX ( I )=PREALI ( I ) 

PREALI (*)=DY(I-I )+DELTY( I) 

DY( I )=PREALI ( I ) 

PREALI (★)=-SIGN 
SIGN=PREALI ( I ) 

★GO  TO  150 
200  ★CONTINUE 
NPTS(*)=I 
★GO  TO  515 
5 00  ★CONTINUE 

★DO  510  Tl=l ,64 
T2=L0CATE(TI ) 

KXMAX(TI )=DX(T2) 

KYMAXITI )=DY(T2) 
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5 10  ^CONTINUE 
515  ★CONTINUE 
520  ^CONTINUE 
★RETURN 
★ END 


I  <KE RR>CFD .GTDATA 1 9 MON  IO-JUN-74  6»4IPM 


PAGE  I 


★SUBROUTINE  GTDATA 

C GTDATA  GETS  ONE  16-BIT  BYT  FROM  THE  INPUT  FILE.  THE  BYTES  ARE 

C DELIVERED  SEQUENTIALLY  STARTING  AT  THE  BEGINNING  OF  THE  DISK  AREA 

C INDM2.  THERE  IS  A CORE  BUFFER  (BUUFI (-))  AND  AN  8 WORD  ADB  BUFFER 

C SO  ALL  MEMORY  ACCESSES  ARE  DONE  VIA  THE  BIN  INSTRUCTION.  THE  BYTI- 

C ENDS  UP  RIGHT  JUSTIFIED  IN  M IN BYT" . 

★PE  INTEGER  NBUFFI  ( ★ , 64 ) , FI NSCN(*> ,COMP( ★) ,T0TSCN(*) , PINTI (★) , 

I P INT2  ( ★) , TI ME  < ★ ) ,OTIME( ★) ,TWTIME(*> ,PEN(*> ,CNTRL( *,6) 

★PE  REAL  GLCHFT(*) ,VARFT(*) ,BUFF2(*, 550,2) ,BUFF3(*, 640) .PREALl (★) , - 

I ppCil  3 f + ) il  I TViDCTf+l 

★PE  INTEGER  L()FREq|hIFREQ  , I BUFF  I ( 4096) , I BUFFS  (*,640) , ABUFF2(  70400) , - 
I CHGOOD ( 80 ) , S I TEGD ( 80 ) ,SITES(80) 

★PE  REAL  CHMSQ(80) ,RBUFFI ( 4096) , ROWSUM, RBUFF2 ( 70400) 

★CU  INTEGEP  ADBBUFf  8) ,COREPT,  BYTE, ADBWRD, ARRAY, DEBUG, TWSZ, 

1 OVLAP,NCHAN,NSITE,NROWS,DIFFR,DIFFW,NEW,OLD,GAP,TSCANS,- 

2 INDEX  I , INDEX2, INDEX3, INDEX4.TI ,T2,T3,T4,T5,T6,CH, I PAGF,  - 

3 OFFSET, INBYT,NGDCH,TWSZR,NGDST,NGDR,F,BF3PE,NGT,OPAGE,T7 
★CU  LOGICAL  LADBBU(8) ,LC0REP,LASTI6,LBYTE,LADBWR,LARRAY,LDEBUG, 

1 LTWSZ,LOVLAP,LNCHAN,LNSITE,LNROWS,LDIFFR,LDIFFW,LNEW,  - 

2 LOLD.LGAP ,LTSCAN  ,LTI ,LT2,LT3, LT4 , LT5 ,LT6,LCH,L0FFSE,  - 

3 L I NBYT.LF , LNGDCH , LTWSZR , LNGDST, LNGDR , LNGT, LT7 
★EXTERNAL  Cl 6T64,C64T32 , ROWSUM, RUN FFT.C32T64 

★C0MM0N/MAIN2/NBUFFI , BUFF2.BUFF3, FI NSCN, COMP, ToTSCN, PINTI ,PINT2, 

1 TIME, OTIME.TWTIME, PEN, CNTRL, GLCHFT, VARFT.PREALI .PREAL2, 

2 A LLMSO, TVARFT , LOFREQ , H I FREQ , S I TEGD .SITES, CHMSQ 
★EQUIVALENCE  (NBUFFI (1,1) .RBUFF I ( I ) . I BUFF  I ( I ) ) . < BUFF2  (1,1,1), 

1 ABUFF2( I >,RBUFF2(I ) ) , ( CHGOODf I ) ,SITEGD( I )), 

2 (BUFF 3(1, l),I BUFF3I 1,1)) 

★EQUIVALENCE  ( I ,ADBBUF( I ) ,LADPBU( I ) ) , ( 9 , COREPT, LCOREP ) , ( I 0, BYTE, 

1 LBYTE) ,( I I , ADBWRD ,LADBWR) ,( I 2, ARRAY.LARRAY) ,( I 3, TWSZ,  - 

2 LTWSZ),(I4,0VLAP,L0VLAP),(I5,NCHAN,LNCHAN),(I6,NSITE,  - 

3 LNSITE) ,( I7,NR0WS,LNR0WS) ,( I B,DI FFR.LDIFFR) ,( I9.DIFFW,  - 

4 LDIFFW),(20,NEW,LNEW),(2I,OLD,LOLD),(22,GAP,LGAP),(23,  - 

5 TSCANS.LTSCAN) , (24, INDEX  I ) , (25, INDEX2) , (26, INDEX3) , (27,- 

6 INDEX4) ,(28,TI ,LTI ) , (29.T2.LT2) , (30.T3.LT3) ,( 31 ,T4, 

7 LT4 ) , ( 32 , T5 , LT5 ) , ( 33, T6, LT6) , ( 34 ,CH,LCH ) , ( 35, OFFSET, 

8 LOFFSE) , ( 36, I NBYT.LINBYT) ,( 37,F,LF) , ( 38,NGDCH,LNCQCH> , - 

9 (39, TWSZR, LTWSZR), (40, NGDST, LNGDST) , (41, NGDR, LNGDR), 

0 ( 42 , FF3PE ) , ( 4 3, NGT, LNGT ) , ( 44 , LASTI 6 ) , ( 45 , DEBUG , LDEBUG ) , - 

1 ( 46.0PAGE ) , ( 47,T7 ,LT7 ) , ( 48, 1 PAGE) 

★DISK  AREA  I N DM2 ( 20) , 0UTDM2 (40) ,CONPRM(  I ) 

★ I F ( DEBUG.LT.  I )G()  TO  10 

A DISPLH  "BGTDATA" ,01 

★ I F ( DEBUG.LT . 4) GO  TO  10 
A DISPLH  ,2l 

10  *CONTINUE 
BYTE*BYTE+I 

★IF(BYTE.LT.5)G0  To  300 
C HAVE  TO  GO  TO  NEXT  ADB  WORD. 

ADBWRD-ADBWRD+I 
★IF( ADBWRD.LT. 9) GO  TO  200 
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C HAVE  TO  REFILL  ADB  BUFFER. 

COREPT*C()REPT+8 

★ IF(COREPT.LT.4097)G()  TO  IOO 

C HAVE  TO  REFILL  BUFFI . 

★ READ! 64, 1 BUFF  I ( I ) , INDM2( I PAGE) ,4) 

★WAIT  64 

IPAGE= IPAGE+4 
CORE PT*I 

★IFtDEBUG.LT. I )G0  TO  20 
A DISPLH  "REF I LLC" ,0| 

★IF (DFBUG.LT. 2) GO  TO  20 
A DISPLH, I 8, IBUFF1 , 1 BUFF  1 ♦255  * 

20  ^CONTINUE 
100  ★CONTINUE 

C CORE  IS  OKAY.  HAVE  TO  REFILL  ADB  BUFFER. 

★ TRANSFERS)  ADBBUFd  )*IBUFF1 (COREPT) 
ADBWRD* I 

★ IFtDEBIJG.LT.  I )G()  To  110 

A DISPLH  " REF  I LLA",Ol 

★ I F ( DEBUG. LT. 2 ) GO  TO  110 
A DISPLH  ,2 l 

110  ^CONTINUE 
200  *CONTINUE 
C NEW  ADB  WORD. 

BYTE* I 

★ IF (DEBUG. LT. I ) CO  TO  2 i 0 

A DISPLH  "NEW  WORD",0l 

★ IFtDEBIJG.LT. 2)G0  TO  210 

A DISPLH  ,2 1 

210  ★CONTINUE 
300  ★CONTINUE 

LADBBU ( ADBWRD) *LADBBU ( ADBWRD) . RTL.  I 6 
LI NBYT=LADBBU( ADBWRD ) .AND.LASTI 6 

★ IFtDEBIJG.LT.  I )G0  To  310 

A DISPLH  "EGETBYT" ,01 

★IF (DEBUG.LT. 4) GO  TO  310 
A DISPLH  ,2» 

310  ★CONTINUE 
★RETURN 
★END 
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★SUBROUTINE  IMG(  IN, OUT) 
♦PE  REAL  IN(*),OUT(*) 
★CU  INTEGER  DEBUG 
★EQUIVALENCE  (3, DEBUG) 
★IF (DFBUG.LT. 5) GO  TO  10 
DISPLH  "IMG"tO$ 

★ I F ( DEBUG.LT. 5) GO  10 
LDL(O)  SD49I 
LDA  I.J(0)« 

DISPLH  "IN" , 32 1 
10  ^CONTINUE 


I 10 


LDL(O)  $D49» 
LDA  IN (0)1 
LDR  $A| 

SHAL  *9* 

SHAR  *57* 

SBM  >40H6I 
ADM  >4000(161 
SHAL  *481 
LDS  $A« 

LDA  $R* 

SHAR  -551 
SHAL  -631 
OR  $Sl 
LDS  SAI 
LDA  $R I 
SHAL  -16* 

SHAH  =401 
SHAL  =24* 

OR  $S* 

LDL(O)  $D50l 
STA  OUT(0)* 


X SAVE  IT. 

X ISOLATE  THE  EXPONENT. 

X SUBTRACT  THE  32  BIT  OFFSET. 

X ADD  THE  64  BIT  OFFSET. 

% PUT  IT  IN  THE  64-BIT  EXPONENT  FIELD. 

X SAVE  IT  FOR  NOW. 

X NOW  FOR  THE  SIGN. 

X START  ISOLATING  THE  SIGN. 

X COMPLETE  ISOLATION  AND  PUT  IT  IN  RIGHT  SPOT. 
X NOW  WE  HAVE  EXPOONENT  AND  SIGN. 

x save  it. 

X RESTORE *()RIGINAL  TO  DO  MANTISSA. 

X START  ISOLATING. 

X REMOVE  OUTER  MANTISSA. 

X COMPLETE  ISOLATION. 

X PUT  THE  THREE  PARTS  TOGETHER. 


★IF (DEBUG.LT. 5) Go  TO  110 
LDL(O)  $D50* 

LDA  OUT( 0) l 
DISPLH  "OUT" ,32 l 
★CONTINUE 
★RETURN 


★END 
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★FUNCTION  M A X ( I ) 

★PE  INTEGER  I ( ★) 

★PE  INTEGER  VAX 
★CU  INTEGER  DEBUG 
★EQUIVALENCE  ( 3, DEBUG) 

★IF (DEBUG.LT. 5) GO  TO  5 
A DISPLH  "B  MAX", Ol 

5 ^CONTINUE 
A LDL(O)  SD49I 

A LDA  1(0)1 

A LIT(0)  = l ,6,1  » 

A LI T( I ) *1 » 

A CLC(3) I 

A COMP( 3) t 

A SETC(2)  El 

A LDEEI  $C3» 

A RTL  $A ,0( I ) l 

A IMG  $R| 

A SETE  -I. AND. El 

A SETE I E.OR.Ei 

A LDA  $R l 

A CSHLd  ) I l 

A TXEFM(O)  ,MAXLO0Pi 

A MAXL(X)P«LDEE1  SC2| 

A STA  MFUNVALl 

★IF (DEBUG.LT. 5) GO  TO  10 
A LDL(O)  $D49| 

A LDA  1(0)1 

A DISPLH  "I",32l 

A DISPLH  "MAX" , I 6,MFUNVAL,MFUNVAL+63l 
10  *CONTINUE 
★RETURN 
★END 


r 

V. 
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♦ PE 

1 

2 

4 

5 

6 
7 

*PE 


REAL 


♦SUBROUTINE  OUTPUT 

PURPOSE*  COMPLETE  CA!.CULATI0NS  ON  THE  POWER  MAXI MIJMS 
3FFORF  PRINTING  OUTPUT. 

AZI MUTH  ARF  CALCULATED  AND  SIGNAL  To  NOISE 
RATIO  AND  FISHER  STATISTIC. 

AND  ™REE  DIMENSIONAL  MAXIMUM  ARE  OUTPUT  IN  SEPARATE 

L i j 1 j i 

DECLARATIONS* 

*PL  INTEGER  I NBUF ( *, 640) , CNTRL( *, 6) , NCHAN( *) , PI NTI (*) ,<)FFSET(*> 
LOCATE (*)  ,NPTS(*)  , COUNT2(*)  ,COUNT3(  *)  ,LoC2D ( * , 25 ) 
LOC3D(*,25),TWTIME(*),ADJF<*) 

POWER!  *,25),FMAX(*,25),FKX(*,25),FKY(*,25),  R I NBUF  ( * , MO) 
X(*,25),Y(*,25),FFT(*,6I2> , KERNEL <*, 25 ) , 

XCOORD(*) , YC()()RD(* ) , PREALI  (★) , PREAL2!*)  ,C;)5K(  *)  .SINK!*) 
C()SDK(*)  ,SINDK(*) , BEAMF.Rf*)  , FPMAX(*)  ,KXMAX(*) 

KYMAX! *) ,DELX( *) , DELY ( *) , KXSEP(  *) , KYSEP(*) ,KSEP(*) , 
VEL(*),AZ(*)  .SIGNAL! ★)  ,FSTAT( ★)  ,SUMSQ(*) 

TEST (★ ) , K( *)  .CHANAV!  ★)  ,TP()WER(  *) , FREQ(  *) 

RFAL  GROUP! (*) .GR0UP2 ! * ) , AFREQ! * , 20 ) , AK ! * ,20) .AFSTAT! *.20) 

1 AAZ!*,20) ,AVEL(*,20) , ACHANA! *,20) ,AFMAX(*,20) , 

2 AGRUPI (*,20>,AGRUP2(*,20),ASIGNA(*,20) 

*PE  REAL  ADKX(4),ADKY(4),YP0INT(50) , YVAX! 50 ) , DX ( 50) ,DY( 50) 

: «TP()W(*)  tDFLTAKU)  ,P,  RPOWER!*,25)  , I POWER!  *,25  ) 

*Kfc  INlhufcK  MAX 

*CU  INTEGER  LoFREQ.HIFREO.DFBUG.SM.TI ,T2 , ARRAY , PAGE, I ,N, MNCHAN, 

1 MNPTS,NP()INT,SWITCH,NFREQ,  I GO,  LINE.  LINES,  INDEX,  I 6 

2 TWIN, ANGLE, SAM, IFREQ,J,NFREQI .RFFINE, IND.YToP 

3 YPMI ,SIGN,NTIMES,LINEPI 
*PE  REAL  DELTX! 3000), DELTY! 3000), DI  ST 
*PE  REAL  DELTAX , DELTA Y , KX , KY 

*CU  REAL  DKX, LOWFR, UPPER, LINER, HDKX, BORDER, TWOH 

*CU  REAL  DEI TAF, RADI  US 

*CU  LOGICAL  M0DE3,NM0DE 

♦EXTERNAL  MAX .FNGRID, GRID , REALE , I MG 

*C()MMON/MAINFK/INBUF,CNTRL,NCHAN,  PI  NTI  , OFFSET,  LOC  ATE,  NPTS, COUNT  2 
C()UNT3,L0C2D,L()C3D,P()WER,FMAX,FKX,FKY,  X,  Y,  KFRNEI  . ’ 

2 XC!)ORD,  Y COORD,  PREAL ! ,COSK,SINK,BFA  ‘ 

3 TP()W,DELTAK,RPOWER,IPOWER,  COSDK.S  I ND* 

t REAMER, FPMAX,KXMAX, KYMAX, DELX, DELY, KXSEP.KYSEP,  r 

-]  TWTIME  ,TPOWER,  VEL,AZ, SIGN AL , FSTAT, SUMSO , TEST,  K, 

7 ^kx^ky  ’ A^F  t DX,DY , P, YPOINT, YMAX, ADKX, ADKY , 

♦ EQUIVALENCE! INBUF! I , I ) , R I NBUF ( 1,1  )) , ( I NBUF! I ,2R) ,FFT( I .1  ) ) 

♦FQUI  VALFNCF  ! I ,L()FREQ) , ( 2,  HIFRFO ) , ! 3, DEBUG) , ( 4 ,SM> , ( 5 ,TI  ) .(6.T2) 

[ ?!.A  RRAXi * ^ » PAGE ),(9,I),(I0,N),  'I  2, MNCHAN ) , 


( I 3, MNPTS) , ( I 4, NPOI NT),!  I 5, SWITCH), (16,100), 

, U7,  INDEX),  (IB.  IP),  (|Q,  DKX), 

(20, LOWER), (21 , UPPER), (22, LINE), (23, LINFS), 

( 24 , HDKX ) , ( 25 , BORDER ) , ( 26 , TWOH ) , 

(27, DELTAF ) , (28, RADIUS) , (2 9, SIGN ) . 

(30,MODF.3),(3I  , NMODE ),  (32,  TWIN),  (33,  ANGLE ),(  34,  SAM), 
(35,NFREQ),(36,IFREQ),(37,J),(38,NFREQI ) , ( 39, REF  I NF ) 
, ( 40 ,NTI  MES  ) , ( 4 1 , 1 ND) , (42,  YT  ()u)  , (43,  YPMI  ) , ( 44  L I NF P ) 
(45.LINEPI)  ’ * 

♦DISK  AREA  CONPRMd  ),STCORD(l  ),FKIN(8I  ) 

♦I F ( DEBUG.LT . I ) GO  TO  10 


3 


8 


0 


o o 
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A DISPLH"OUTPUT'\Ol 

♦ IF (DEBUG.LT.  I )G()  TO  10 
A DISPLH  "C0UNT2" , I 6,C()UNT2tCOUNT2+l Ol 

A DISPLH  "C0UNT3"  , I 6,COUNT3,COUNT3+l  Ol 

A DISPLH  "L0C2DM  ) ",l 6, L0C2D, L0C2D+ I Ol 

A DISPLH  "L0C3D(  I ) ",  I 6,L()C3D,LOC3D+l  Ol 

I 0* CONTINUE 

fMMCuTJIi°ci^TM^CALCULATE  AZIM,JTH  NORTH  AND  VELOCITY 
IN  1C M PER  wECOiiD • 

MODE-ON 

P I NT 1 C *)*MAX(CO0NT2! ★ ) ) 

TI-PINTKI  > 

★DO  200  I-l ,TI 

MODE-MODE .AND. ( I .LE .C0UNT2 (★ ) ) 

MBIT2-M0DE 

PINTI (*)»L0C2D(*,I ) 

FM AX ( ★ t PI NTI ( ★ ) ) -FMAX ( ★ , P I NT  I ( ★ ) ) /FLOAT! NCHAN ( ★ ) *TWI N ) **2 
AZ(*)=ATAN(FKX(*, PINTI (★)) /FKY(* , PINTI (*) )) 

'4BITl-(FKY(*#  PINTI  (*)).LT.O.O) 

MODE-MODE. AND. MBITI 
AZ(*)-AZ(*)+3. 14159265 
M0DE-MBIT2 
AZ(*)-AZ (*)*57 . 3 
C RADIANS  TO  DEGREES. 

★IF <( AZ (★). LT. 0.0 ))AZ(*)=AZ< *>+360.0 
K( *)-SQRT(FKX(*t PINTI (★) ) **2+FKY! ★, PINTI  (*>)**2> 

VEL(*> -FLOAT! PINTI (★)-! ) *DELTAF/K(*> 

CHANAV (*)=0.0 

OFFSET!*)-! PINTI (★) -LOFREQ)*NCHAN( ★) 

MB IT I -MODE 
★DO  100  T2-1.M NCHAN 
MODE-MODE . AND . ( T2 . LE . NCHAN ( ★ > > 

★CALL  REALE ( FFT! ★ , OFFSET! ★) +T2  > , PREAL I ( ★)) 

★CALL  IMG! FFT! ★, OFFSET! ★) +T2) , PREALZ! ★) ) 

CHANAV (★) -CHANAV!* J+PREAL1 !★) ★★2+PREAL2!*) **2 
I 00*C()NTINUE 
MODE-MBITI 

CHANAV ! ★ J-CHANAV ! * ) /FLOAT! NCHAN! ★>  *TWI N*TWIN  > 
SIGNAL(*)»FMAX(*, PINTI !★))/! CHANAV (★) -FMAX <★, PINTI  (★))) 
FSTAT! ★) » ( FLOAT (NCHAN ( ★ ) )-l .0>*SIGNAL(*> 

GROUP I (★>-SQRT((FKX(*, PINTI !★) )-FKX( ★, PINTI (*>-l ) ) **2+ 

' „ „ !FKY(*, PINTI (*))-FKY!*, PINTI (*>-| ))**2) 

GROUPI (★J-DELTAF/GROUPI (★) 

GR0UP2(*)=SQRT( ! FKX! ★, PI  NT I !★) J-FKX! ★, PINTI (★)+!  ) ) **2+ 

' _ !FKY(*jPINTI !*))-FKY(*tPINTl (★)♦! ) )**2) 

GR0UP2 ( ★ ) -DELTAF/GR0UP2 ! ★ > 

AFREQ! ★, I ) -FLOAT! PI NT I (★)) 

AK ( ★ , I ) — K ( ★ ) 

AFSTAT!*, I J-FSTAT!*) 

AAZ! ★, I ) — AZ ! ★ ) 

AVEL(*tI )=VEL(*> 

ACHANA!*, I J-CHANAV! ★> 

AFMAX! ★, I J-FMAX! ★, PI NTI ( ★ ) ) 

AGRUPI (★, I ) -GROUP I (★) 

AGRUP2 !★, I J-GR0UP2!*) 

ASIGNA!*, I J-SIGNAL!*) 

★IF!DEBUG.LT. I ) GO  TO  200 
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A DISPLH  "2-D  MAX" ,01 

A SETC(O)  El 

A DISPLH  "MODE", I I 

A DISPLH  "FREQ", I6,PINTI ,PINTI+63l 

A DISPLF  "K" , I 6, K, K+631 

A DISPLF  "FSTAP' ,!6,FSTAT,FSTAT+63l 

A DISPLF  "AZ»,l6,AZ,AZ«-63l 

A DISPLF  "VEL" , I 6, VEL, VEL+631 

A DISPLF  "CHANAV" , I 6,CHANAV, CHANAV+63t 

A DISPLF  "SIGNAL", 16, SIGNAL, SIGNAL+631 

PREALI  (*)®FMAX( ★, PINTI (★>> 

A DISPLF  "FMAX" , I 6, PREALI , PREALI ♦63l 

200*C()NTINUE 
MODE-ON 

★DO  250  1=1,64 
Tl  =NCHAN ( I ) 

★ IF(TI  .LT.I )G0  TO  245 
T1=C0UNT2( I) 

★IF (Tl .LT.I )G0  TO  245 
PINTI ( 1 )=TWTIME( I) 

A DISPLH  "7WTIME", I 6, PINTI , PINTI I 

★DO  240  J-l ,T1 
PREALI (1 )=AFRE0(I,J) 

PREALI (2 )=AK ( I , J ) 

PREALI ( 3 )“AFSTAT( I , J ) 

PREALI <4>«AAZ( I, J) 

PREALI (5)«AVEL(I,J) 

PREALI ( 6)=ACHANA( I , J ) 

PREALI ( 7 )=AFMAX( I , J ) 

PREALI (8>=AGRUP1  (I ,J) 

PREALI ( 9)=AGRUP2 ( I , J ) 

PREALI ( 10)=ASIGNA( I,J) 

A SKIP  , E2DI I 

A B2DI IDATA 

A "FREQ  K FSTAT  AZ"l 

A E2DI  * * I 

A DISPLS  , 1 6,B2D1 , E2DI -I  I 

A DISPLF,! 6, PREALI  ,PREALI+3l 

A SKIP  ,E2D2l 

A B2D2 iDATA 

A "VEL  CHANAV  FMAX"! 

A E2D2 ««  I 

A DISPLS  , I6,B2D2,E2D2-I I 

A DISPLF  , I6,PREALI+4,PREALI+6| 

A SKIP  ,E2D3l 

A B2D3 «DATA 

A "GROUP  I GR0UP2  SIGNAL/NOI SE" t 

A E2D3 « « I 

A DISPLS  , I6,B2D3,E2D3-I I 

A DISPLF  , I 6, PREALI +7, PREALI +9 1 

240  ^CONTINUE 


J 
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245  ★CONTINUE 
250  *C0NTINIJE 

PINTI  (★>-MAX(COUNT3(*>> 

II  =P INTI ( I ) 

★DO  400  1=1 ,TI 

MODE=MODE. AND.  ( I .LE .C0UNT3! ★ ) ) 

MBI T2=M0DE 

PINTI  ( *)  = L()C3D(  *,  I ) 

FMAX(*t PINTI (*>  >=FMAX(*f PINTI (★) ) /FLOAT! NCHAN! ★) *TWI N ) **2 
AZ(*)=ATAN(FKX!*, PINTI (★)) /FKY! ★.PINTI  (★))) 

MBITI=(FKY (★, PINTI (★)) .LT.O.O) 

MODE=MODH. AND. MBITI 
AZ(*)=AZ(*)+3. 14159265 
M0DE=MBIT2 
AZ ( ★ ) = AZ( *)*57 .3 
C RADIANS  TO  DEGREES. 

★ IF ( v AZ (*) .LT.O.O ) ) AZ( ★) =AZ( ★ ) + 360. 0 
' *■ *)=SQRT(FKX( ★, PINTI (★) ) **2+FKY (★, PINTI  (★))★*?) 

VEL!  *)=FL()AT(  PINTI  (*)-l  ) ★DELTAF/K!  ★) 

CHANAV ( * )=0.0 

OFFSEr(*)=( PINTI (★)-LOFREQ)*NCHAN( ★) 

MB ITI=MODE 
★DO  300  T2=l tMNCHAN 
MODE=M()DE.AND.  (T2.LE  .NCHAN ( ★ )) 

★CALL  REALE ( FFT( *, OFFSET! ★) +T2 ) , PREALI ( ★) ) 

★CA  LL  I MG( FFT( * , OFFSET! ★) +T2) , PREAL2 ( ★ ) ) 

CHANAV (* )=CHANAV( *)+PREALI ( ★ ) **2+PREAL2 ( ★) +*2 
300*C()NTINUE 
MODE=MBITI 

CHANAV ( * )=CHAN AV ( ★ ) /FLOAT! NCHAN! ★ ) *TWI N*TW I N ) 
SIGNAL(*)=FMAX(*. PINTI (★))/( CHANAV! ★ )-FMAX( ★, PINTI  (★) )) 
FSTAT!  ★)  = ( FLOAT (NCHAN!  ★)  )-l  .0)*SIGNALU) 

GROUP! ( ★ )=SQRT( ( FKX! ★, PI NTI ( ★ ) ) -FKX (★, PINTI  ( ★) -I  ))**2+ 

I (FKY!*, PINTI  (★)) -FKY (★, PINTI  (★)->!  ))★★?) 

GROUP  I (*)=DELTAF/GR0UPI  (★) 

GR0UP2 ( ★) =SQRT ((FKX!*, PI  NT I!*) ) -FKX (*, PINTI (*)  + | ))**2+ 

I (FKY (★, PINTI (*)) -FKY (*,P INTI (*)+l ))**2) 

GR0UP2 ( ★ )=DELTAF/GR0UP2(  ★) 

AFREQ! ★, I )=FLOAT( PINTI (★)) 

AK ( ★ . I )=K< ★) 

AFSTAT!*, I )=FSTAT(*) 

AAZ! ★, I )=AZ( ★) 

A VEL (* , I )=VEL( ★) 

ACHANA!*, I )=CHANAV(*) 

AFMAX! *,  I ) =FMA  X ( *,PI NTI ( ★ ) ) 

AGRUPI (*, I )=GROUP I (★) 

AGRUP2  ( ★ 1 1 )=GR()UP2  ( ★ ) 

AS IGNA ( *, I )=SIGNAL( ★ ) 

A DISPLH  "3-D  MAX" ,01 

A SETC(O)  El 

A DISPLH  "MODE", I » 


C7< 


¥ 


I 
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A DISPLH  "FREQ", 16.PINTI ,PINTI+63» 

A DISPLF  "K" , 1 6,K,K+63t 

A DISPLF  MFSTAT"  , 1 6, FSTAT, FSTAT+63* 

A DISPLF  "AZ" , I 6, AZ,AZ+63l 

A DISPLF  "VEL" , I 6, VEL, VEL+631 

A DISPLF  "CHANAV"  , I 6,CHANAV,CHANAV+63l 

A DISPLF  "SIGNAL", 1 6, SIGNAL, SIGN AL+631 

PREALI  (★)=FMAX(*,PINTI  (★)) 

A DISPLF  "FMAX",  I 6, PREALI .PREALI +631 

A DISPLF  "SIGNAL", 16, SIGNAL, SIGNAL+631 

400+CONTINUE 
M()DE=ON 

★DO  450  1=1,64 
TI=NCHAN C I ) 

★ IFCTI.LT.  I )G()  TO  445 
TI=C0UNT3(I) 

★IFCTI.LT. I )G0  TO  445 
PINT1  ( 1 )=TWTIMECI) 

A DISPLH  "TWTIME",  I6.PINTI  ,PINTI  t 

★DO  440  J=! ,TI 

PREALI Cl  )=AFREQ(  I ,J) 

PREALI (2 )=AK ( I , J) 

PREALI  C3)=AFSTAT(I,J) 

PREALI  ( 4 ) = AAZC I , J) 

PREALi  ( 5 )=AVEL( I , J ) 

PREALI (6)=ACHANA(  I ,J) 

PREALI  ( 7 )=AFMAX( I , J) 

PREALI  (8 )=AGRUPI ( I , J ) 

PREALI  C 9 )=AGRUP2(  I , J ) 

PREALI ( 10)=ASIGNA(I,J) 

A DISPLS  , I 6.B2D1 . E2D I —I  I 

A DISPLF,!  6, PREALI .PREALI+31 

A DISPLS  , I 6,B2D2,E2D2-1  I 

A DISPLF  , 1 6, PREALI + 4, PREALI  +61 

A DISPLS  , 1 6,B2D3,E2D3-I  I 

A DISPLF  , 16,PREALI+7,PREAL1+9l 

440  ^CONTINUE 
445  ★CONTINUE 
450  ★CONTINUE 
★RETURN 
★END 
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♦SUBROUTINE  PUTBYT 

n & SSeTofT¥,1 

♦ PE 


1 

2 

*CU 

1 

2 
3 

♦ CU 

1 

2 
3 


INTEGER 


LOGICAL 


INTEGER  CNTRL( *,6 ) , OtJTBUF  ( * , 64, 6) ,PI NT  I (*)  , INBIJFC  ♦.  I 28) 

TIME(*) ,OLDTIM(*),  * * 

OUPTWA , 6,  .SAVBCT,’,?K,,UP' **' " l”' .SCANS, 

* v ADB  • t??  ^ • 

So.BCT.ADBy^ 

r'RYTq^wnDnq^RTy,H^P7?;L,IiPTW,LSAVAD,LADRn[I(6)*L(,um' 

• LT2tLT3,LT4#LT5 ,LT6,L()IJPYT.LIT,LPRTIA 

LADDRS,LWORD,  LINBYT,LBYTCN(6) , LADBWR,L()RGCO,LPAGE 
,LUEBUG,LBCT,LADB,LENDAD  ’ 

♦external  rdprm,getbyt,cmvtim 

♦COMM ON /MA I N/CNTRL , OIJTBUF « I NBUF , P INTI  ,TI ME  f f)LDTI M SA VBCT, S AVPTW 
I TSTEPS,  SCANS,  OIJPTWA,  OUPAGE,  oTIMEA  ' 

♦ EOUI  VALFNCE(  I , ADB3IJF Cl), LADBBU ( I ) ) , ( 9,  ARRA Y , LA RRAY) . 

I ( IO,INPTB,LINPTR) , 

! (II , I NPTW , LINPTW) , { I2,SAVADB,LSAVAD) 

M3,  ADBOUT(  I ) ,LADBOU(  I )) 

! \\  V/Ttm  ' ! 5 ,*l?°  ’ ^TS ' LBYTS  } ’ ( 2 1 * W(,RDS  • LWORDS ) . 

’rLI'i  • 23,T2,LT2),(24,T3,LT3)t(25,T4,LT4), 

( 26 ,T5 ,LT5 ),(27,T6, LT6) , ( 28 , OUBYT, LOUBYT) ,(29. IT  LIT) 

J 5 * R JAL’  ) ! 5, 3 1 ♦ ADDRS  t LA  DDRS ) , ( 32,  WORD , LWORD ) , ' 

( 33,INE?YT,LINBYT)  ,(34,BYTCNT(  I ),LBYTCN(I  )). 

( 40, ADBWRD.LADBWR) 

a , Dkr.n  • ( 43 , ORGCOR , LORGCO ) , 

^•Ai2E;LPA<?E)*(45»DEBUG»LDEBlJG)*(46.BcT,LBCT) , 
(47,ADB,L/  jB)  ,(48,ENDADB,LENDAD) 

* 1 AREA  ()|jpuj6(  20)  ^ VnPUT^5 o?  J * ()UP^T3(^0>  «OUPUT4(20)  , 0UPUT5  ( 20)  , 

♦I  F (DEBUG. LT. I )G0  TO  5 
OISPLH  " PIJTBYT", 0? 

♦I  E (DEBUG. LT.  3) GO  T()  5 
DISPLH  "PUTBYT* " ,2! 

♦CONTINUE 

FIRST  WE  ROTATE  ADB()UT( ARRAY ) 

BYTCNT(ARRAY)  BY  I. 

L ADS OU  (ARRAY  ) = L.ADB()U ( ARRAY ) .RTL.  I 6 
LADBOU(ARRAY  )=LADB()U(  ARRAY  ) .OR  .LOUBYT 
kDUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS 
BYTCNT( ARRAY)=BYTCNT(ARRAY)+I  * 

♦IF (DEBUC.LT. 3)G0  TO  15 
DISPLH  "PUTBYT I " ,2  I 
♦CONTINUE 

♦IF  (BYTCNT(ARRAY) .LT.4JG0  TO  100 


LEFT  16  AND  STICK  IN  OUBYT  AND  BUMP 


OOO  o o o >*oo 
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ADB  WORD  IS  FULL. 

OUPTW-OUPTWA  ( ARRAY) 

DELAY  WE  RUN  INT°  TH°  CFD  RESTRICTI()N  ()N  DIMENSIONING.  EXCUSE  THE 

LT6= LOUPTW.SHR .6 
LT5* OFF. TURN  ON.. LAST. 6 

/.DUMMY  ASK  STATEMENT  TO  FORCE  DEALLOCATION  OF  REGISTERS. 

LT5= LT5 . AND . LOUPTW 

0UTBUF(T5+I  , T6-H  , ARRAY )=ADBOUT( ARRAY ) 

BYTCNT( ARRAY )=0 
ADB0UT(ARRAY)=0 
T6=0UPTW+I 
0UPTWA(ARRAY)=T6 
*IF(OUPTW+l  .LT.4096)G0  TO  100 
OUTBUF( ARRAY)  IS  FULL. 

♦ IF (DEBUG. LT. I )G0  TO  25 
DISPLH  "COREFULL '',01 

*IF(DEBUG.LT.3)G0  TO  25 
DISPLH  « ",  2 1 
25  *C0NTINUE 

T6=0UPAGE(  ARRAY) 

*1 F ( ARHAY.EO.  I ) WRITE ( 64,0UTBUF(  1,1,1), OUPUTI (T6),3) 

♦ IF (ARRAY. EQ.2)WRITE( 64,0UTBUF(  1 , 1 ,2)  .OUPUT2CT6) ,3) 

*IF ( ARRAY. EO. 3) WRITE ( 64,OUTBUF( 1,1,3),  0UPUT3( T6) ,3) 

♦ IF (ARRAY .E0.4)WRITE( 64,0UTBUF(  1,1,4), 0UPUT4 ( T6 ) , 3 ) 
*IF(ARRAY.EQ.5)WRITE(64,0UTBUF(I ,1 ,5) ,OUPUT5(T6) ,3) 

♦I F ( ARRAY. EQ. 6 ) WRITE ( 64, OUTBUF ( 1,1,6), 01JPUT6CT6 ) , 3) 

♦WAIT  64 

T6=0UPAGE(  ARRAY) 

T6=T6+3 

OUPAGE ( A RRA Y )=T6 

WE  WROTE  OUT  EVERYTHING  BUT  THE  LAST  PAGE  OF  THE  BUFFER.  SAVE 
THAT  AT  THE  TOP  OF  THE  BUFFER  FOR  FUTURE  USE. 

♦DO  40  T6«l , I 6 

40  OUTBUF (*,T6, ARRAY )=OUTBUF(+,T6+48, ARRAY) 

OUPTWA 'ARRAY)* 1024 

THE  RcST  OF  THE  BUFFER  REALLY  DOESNT  HAVE  TO  BE  ZFROED  OUT, 

BUT  IT  EXECUTES  SUPER  FAST  AND  MIGHT  SAVE  SOME  CONFUSION 
SOME  DAY. 

♦DO  50  T6* 17,64 
50  OUTBUF ( *,T6, ARRAY )*0 
00  ♦CONTINUE 

*IF(DEB!JG.LT.  1)  GO  TO  105 
DISPLH  "EPUTBYT * " ,0| 

*IF(DEBUG.LT. 3) GO  TO  105 
DISPLH  " ",  2 l 
05  ♦CONTINUE 

♦ RETURN 

♦ END 


70< 


j 
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★SUBROUTINE  RDPRM 

★PE  INTEGER  CNTRL( * , 6 ) , (>UTBUFU,64,  6) , P I NTI (★) , INBUFC ★, 1 28) . 

I TIMEU)  .OLDTIMU), 

1 SAVBCT, SAVPTW t OUPAGE (A) t TSTEPS( 6) .SCANS , 

2 OUPTWA(6),  OTIMEA(6)  .ORGADB,  INRUFI  (8192) 

*CU  I NTEGER  ADBBUF (8) , A RRAY , I NPTB , I NPTW , SAV fDB , ADBOUTI 6) , OUPTW, 

1 BYTS. WORDS, TI.T2.T3.T4. TC, T6,  IT, PRTI AL, ADDRS, 

2 WORD,  BYTCNT(6),ADBWRD,INBYT,OllRYT,()RGCOR,PAGE, 

3 DEBUG, BCT, ADR, ENDADB 

★CIJ  LOGICAL  LADBBU(  8)  ,LARRAY,LI NPTB.L  INPTW,LSAVa3,LADB()U(  6)  ,L()UPTW 

1 L BYTS , LWORDS , LT I , LT2 , LT3 , LT4 , LT5 , LT6, LOUB YT, L IT, LPRTI A 

2 LADDRS , LWORD.L I NBYT ,LBYTCN(  6 ) , L ADBWR, I.ORGCO, LPAGE 

3 , LDEBUG.LRCT, LADB.LENDAD 
★EXTERNAL  GETBYT.PUTBYT.CNVTIM 

★ COMMON/MAI N/CNTRL.OUTBUF , INBMF.PINTI  ,TI ME, OLDTIM, SAVBCT, SAVPTW, 

I TSTEPS, SCANS, OUPTWA, OUPAGE, OTIMEA, ORGADB 

★ EQIJI  VALENCEt  I ,ADBB!JF(  I ) , LADRBUI I ) ) , ( 9,  ARRAY, LARRAY ) , 

I ( I 0, 1 NPTB , L I NPTB) , 

I ( M ,INPTW,LINPTW),(I2,SAVADB,LSAVAD), 

1 ( I3,ADB0UT(6) ,LADB0U(6)) 

2 , (I  9, OUPTW  .LOUPTW ),( 20, BYTS, LBYTS),C2 1 .WORDS, LWORDS), 

3 <22,TI,LTI ),<23,T2,LT2),(24,T3,LT3),(25,T4,LT4), 

4 (26,T5 ,LT5 ),(27,T6,LT6),(28, OUBYT .LOUBYT) , (29, IT, LIT) , 

5 ( 30, PRTI AL, LPRTI A), (31  , ADDRS, LADDRS) , (32.WORD.LWORD) , 

6 ( 33, 1 NBYT, LI NBYT) , ( 34,BYTCNT( I ) ,LBYTCN( I ) ) , 

6 (40,ADBWRD,LADBWR) 

7 , (43,ORGCOR,LORGCO) 

8 ( 44 , PAGE , LPAGE) , ( 45  .DEBUG,  LDEBUG) , ( 46 , F.CT,  LBCT) , 

9 (47, ADB.LADB) ,(48, ENDADB, LENDAD) 

★EQUIVALENCE  ( INBUF  ( I ,1  ),INB(JF1  (I  )) 

★DISK  AREA  OUPUTI  (20)  ,OUPUT2(20)  ,0UPUT3(20)  ,()UPUT4(20)  ,()UPUT5(20)  . 

I ()UPUT6(20)  , INPUT(50) 

DEBUG=0 

★IF  (DEBUG.LT. I ) GO  TO  10 
DISPLH ''RDPRM'', 0| 

UISPLH  "DEBUG'', 2 1 
10  *CONTINUE 
★RETURN 
★END 
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★SUBROUTINE  REALE( IN, OUT) 
★PE  REAL  IN( ★) , OUT ( ★) 

★CU  INTEGER  DEBUG 
★EQUIVALENCE  (3, DEBUG) 

★IF (DEBUG.LT. 5 ) GO  TO  10 
DISPLH  "REALE", 01 
★ I F ( DEB UG . LT . 5 ) GO  TO  10 


% SAVE  IT. 

a ELIMINATE  THE  SIGN  SO  WE  CAN  Do  THF 
« EXPONENT. 

* ISOLATE  THE  EXPONENT. 

« SUBTRACT  THE  32  BIT  OFFSET. 

% ADD  THE  64  BIT  OFFSET. 

% PUT  IT  IN  THE  64  RIT  EXP.  FIELD. 

* nowEfor*the  SIGN. 

* SIGN  BIT  IS  ISOLATED. 

% NOW  WE  HAVE  EXPONENT  AND  SIGN. 

Si  SAVE  IT. 

% NOW  FOR  THE  MANTISSA. 

* ISOLATE  THE  MANTISSA. 

% PUT  IT  IN  64  BIT  EXP  FIELD. 

% DONE . 


★IF (DEBUG.LT. 5) GO  TO  110 
DISPLH  "EREALE" ,0| 

★IF (DEBUG.LT. 5) GO  To  110 
LDL(O)  $D50| 

LDA  0UT(0)» 

DISPLH  "OUT" ,32 1 
10  ^CONTINUE 
★RETURN 
★END 


LDL(O)  $D49| 

LDA  IN (0)| 
DISPLH  "IN", 32 l 
10  ★CONTINUE 
LDL(O)  $D49| 

LDA  IN (0)1 
LDR  $A» 

RAB  =01 

SHAR  =561 
SBM  >40H6l 
ADM  =4000«16| 
SHAL  =481 
LDS  $A I 
LDA  SHI 
SHAR  =63 1 
SHAL  =63 1 
OR  $S» 

LDS  $A  I 
LDA  SRI 
SHAL  =401 
SHAR  =161 
OR  $S| 

LDL(O)  SD50I 
STA  OUT(O) I 
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★FUNCTION  ROWSIJM(R) 

*PF  WEAL  R.ROWSUM 
★CU  INTEGER  DEBUG 
★EQUIVALENCE  (45, DEBUG) 

★CONTINUE 
A Ll)L(O)  SD49I 

A LDA  W(0)l 

A LIT(O)  1,6,11 

A LIT< I)  «ll 

A WONSUMLOOP* 

A WTL  SA ,0( I )l 

A ADWN  SRI 

A CSHL( I ) ll 

A TXEFP(O)  , ROWSUMI.OOPi 

A STA  RFUNVAH 

★IFIDEBUG.LT. I ) GO  TO  20 
A DISPLH  "EROWSUM" ,0| 

★Ir(DFBUG.LT.2)G0  TO  20 
1 DISPLH  "ROWSUM", I 6,RFUNVAL,RFUNVAL+63t 

A LDL(O)  $D49| 

A LDA  R( 0) I 

A DISPLH  «R",32l 

2n  ★CONTINUE 
★I  ETUWN 
★ END 
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★SUBROUTINE  RUN FFT 

^SUBROUTINE  CALLS  THE  UNIVERSITY  OF  ILLINOIS  (JIM  STEVFNS) 
EFTQR(,^Jj.NE  AFTER  SETTING  UP  THE  PROPPER  PARAMETERS.  THE  DATA 
T(i?LFFrEP  IS  IN  RUFF2‘  IT  STARTS  AT  EITHER  PUFF2+0  OR  R1IFF2 
+ 35200  DEPENDING  ON  THE  VALUE  OF  "NEW".  THE  NUMBER  OF  TIMEWI ND()W‘i 
IS  " NGDCH"  AND  THE  TIME  WINDOW  SIZE  IS  IN  "TWSZ". 

★PE  INTEGER  NBUFF  I (★,  64 ) t FI  NSCNU)  ,COMP(  ★)  ,T()TSCN  ( ★)  . PINTI  ( ★)  . 

ipF  REAL  SAVE(I!I2<*)*TIME<*>*0TIMEU,,THTIME<*),PEN<*,fCNTRLU*rt) 

★PE  REAL  GLCHFTU)  ,VARFT(*)  ,BUFF2U,  550,2) , BUFF3( *,640)  ,PREALI  (★)  , - 
I PREAL2( ★) ,ALLMSO(*) ,TVARFT(*) 

★PE  INTEGER  LOFREQ ,HIFREQ , I BUFF  I ( 4096) , IBUFF3( ★, 640) , APUFF2( 70400 ) - 

1 CHGOOD(80) ,SITEGD(80) ,SITES(80) 

★PE  REAL  CHMSQ( 80) ,RBUFFI (4096), ROWSUM, RBUFF2 ( 70400) 

★CU  REAL  RADBBU(8) 

★CU  INTEGER  ADBBUF(8) ,COREPT,  BYTE, ADBWRD, ARRAY .DEBUG. TWSZ, 

' OVLAP, NCHAN.NSITE .NROWS.DI  FFR.DI  FFW,NEW,()LD,GAP,T5CANS  - 

2 INDEXI , INDEX2, INDEX3, INDEX4.TI  ,T2,T3.T4,T5,T6,CH, 

3«..  . OFFSET,  INBYT, NGDCH, TWSZR,NGDST,NGDR,F,BF3PE.NGT. PAGE. T7 

★CU  LOGICAL  LADBBU( 8)  ,LC()REP,  LASTI 6,LBYTE,LADBWR,LARRAY,LDEBUG, 

L LT^SZ^LOVLAP.LNCHAN.LNSITE.LNROWS.LDIFFR.LDIFFW.LNEW,  - 

\ Kb£itG,A^ISCAN  .LTI,LT2,LT3,LT4,LT5,LTrt,LCH,LOFFSE  - 

^YT,DM1,  L NBYT,LF,LNGDCH,LTWSZR,LNGDST,LN0DR,LNGT,LT7 

★EXTERNAL  GTDATA , C I 6T64 , C64T32 , ROWSUM ,C32T64 

★COMMON /MAIN2/NBUFF I ,BUFF2,BUFF3,FINSCN, COMP, TOTSCN, PINTI  , PINT2 

1 TIME, OTIME.TWTIME, PEN, CNTRL.GLCHFT.VARFT.PREALI  PREAI  2 

2 ITw  A LLMSO, TVARFT , LOEREQ , H I FREQ , S ITEGD , S I TES , CH MSO  * ' 

★EQUI VALENCE  ( NBUFF 1(1,1), RBUFF I ( I ) , I BUFF  1(1  ) ) , ( BUFF2 (1,1,1), 

1 ABUFF2(  I ) ,RBUFF2(  I ) ) , ( CHG(X)D(  I ) ,SITFGD(  I )), 

2 ( BUFF 3 ( I , I ) , I BUFF 3(1,1  ) ) * 

★EQUI VALENCE(  I , RADBBU( I ) ) 

★EQUIVALENCE  ( I ,ADBBUF(  I ) ,LADBBU(  I ) ) , ( 9,C()REPT,LC0REP> , ( I O.BYTE, 

' l»A[JlBWRD,LADBWR)t(  12, ARRAY, LARRAY),(  13, TWSZ,  - 

\ f;°VLAP,LOVLAP),(l5,NCHAN,LNCHAN),M6tNSITE,  - 

t LDIffW),(20,NEW,LNEW),(2I  .’oLD.’loLD)  : (22,GAP,LGAP) , (23;  - 

5 TSCANS ,LTSCAN ) , (24 , INDEXI ) , ( 25 , INDEX2) , ( 26, INDEX3) , (27, - 

6 INDEX4),(28,Tl,LT!),(29fT2,LT2),(30,T3,LT3),(3I.T4  - 

Z fL^*J.32»P*LT5)  • ( 33, T6,LTft) , (34,CH,LCH) , (35, OFFSET, 

9 ^^qF5E^  • p3®  * .L^BYT) , (37,F ,LF) , ( 38, NGDCH, I.NGDCH) , - 

n 39,TWSZR,LTWSZR),(40,NGDST,LNGDST),(4| ,NGDR,LNGDR) f - 

? ! 1 ! "3,NGT,LNCT> , ( 44, LASTI  6) , (45, DEBUG, LOEBIJG)  ,- 

• i 46, PAGE) , (47, T7, LT7) 

★DISK  AREA  INDM2( 20) ,0UTDM2( 40) ,CONPRM( I ) 

★IF(DEBUG.LT.  I )G()  TO  I 0 
DISPLH  " RUNFFT" , Of 
*IF(DEBUG.LT.2)G0  TO  10 
DISPLH  ,2l 
DISPLH  "ERUNFFT" ,0l 
10  ★CONTINUE 

★TRANSFERS)  SAVE ( I )«RADPBU(  I ) 
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SAVE  THE  ADB  LOCATIONS  DESTROYED  BY  FFT  ROUTINE. 
A SLIT(O)  SIZE* 

A LDL!I)  TNSZl 

A STORE! 0)  $CI» 

A SI.IT(O)  NUMBER* 

A l.OL(l)  NGDCHl 

A STORE! 0)  SCI* 

A SLIT(O)  ARGLIST+2* 

A SLIT! I > BUFF2* 

A LDL! 2)  OLD* 

A ALIT12)  -I* 

A 7.ERXK2)  ,1* 

A ALIT(I)  35200* 

A STORE! 0)  SCI* 

A EXTERNAL  FFT* 

♦IF IDE3UG.LT. 2)00  TO  20 
A SLIT! I > ARGLIST+2* 

A LOAD! I ) SCI* 

A CSHR! I ) 6* 

A LOA  0(1)1 

A SLITI3)  ARGI.IST* 

A LOAD! 3)  SCO* 

A LOAD! 0)  SCO* 

A ALITI3)  I* 

A LOAD! 3)  SCI* 

A LOAD! I ) SCI* 

A ALITI3)  ll 

A LOAD! 3)  SC2* 

A DISPLH  *B4  FFT*", 33* 

20  *CONTINUE 
A CLC! 3) * 

A SLIT! 3)  FFT* 

A SLIT!2)  ARGLIST* 

A EXCHL! 3)  S ICR* 

A CACRB  10* 

♦IFIDEBUG.LT. 2)00  TO  30 
A SLIT! I ) ARGLIST42* 

LOAD! I ) SCI* 

A CSHR! I ) 6* 

a LDA  0! I )* 

A DISPLh  "AFTER**, 33* 

30  "CONTINUE 
A SKIP  .EARGLIST* 

A ARGLIST* 

A DATA  SIZE.NUMBER.O* 

SIZE • ADS  I* 

A NUMBER  I NDS  I* 

A EARGLIST** 

♦TRANSFER! 8)  RADBBU!  I )*SAVE( I ) 

I 10  *CONTINUE 
♦RETURN 
•END 
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